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In  this  research,  we  design  and  develop  a  virtual  cell  approach  for  the  transmis- 
sion of  IP  datagrams  in  mobile  computer  communications.  A  virtual  cell  consists  of 
a  group  of  physical  cells  whose  base  stations  are  implemented  by  remote  bridges  and 
interconnected  via  high-speed  datagram  packet-switched  networks.  Host  mobility  is 
supported  at  the  data  link  layer  using  the  distributed  hierarchical  location  informa- 
tion of  mobile  hosts.  It  eliminates  the  necessity  of  IP-level  mobile  host  protocols  that 
may  interfere  with  the  conventional  IP  protocol  in  a  practical  sense  and  achieves  a 
logically  flexible  coverage  area  according  to  mobility  and  communication  patterns 
among  physical  cells. 

Given  mobility  and  communication  patterns  among  physical  cells,  the  problem  of 
deploying  virtual  cells  is  transformed  to  the  optimization  problem  of  finding  a  cover 
of  disjoint  clusters  of  physical  cells.  The  objective  is  to  minimize  the  total  communi- 
cation cost  for  the  entire  system  where  intercluster  communication  is  more  expensive 
than  intracluster  communication.  Our  problem  differs  from  general  graph  partition- 
ing problems  in  that  it  must  meet  the  underlying  topology  constraints,  such  as  the  | 
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linear  arrangement  of  physical  cells  in  highway  cellular  systems  and  the  hexagonal 
mesh  arrangement  of  physical  cells  in  cellular  systems. 

For  highway  cellular  systems,  we  design  an  efficient  optimal  partitioning  algorithm 
of  O(mn^)  by  dynamic  programming,  where  m  is  the  number  of  clusters  in  a  partition 
and  n  is  the  number  of  base  stations.  For  hexagonal  cellular  systems,  we  develop 
several  heuristics  for  multiway  partitioning,  based  on  the  techniques  of  moving  or 
interchanging  boundary  nodes  between  adjacent  clusters.  These  heuristics  produce 
optimal  partitions  with  respect  to  the  initial  partitions  obtained  randomly  or  by 
centering.  The  heuristics  are  compared  and  shown  to  behave  quite  well  through 
experimental  testing  and  analysis. 

Once  an  optimal  partition  of  disjoint  clusters  is  obtained,  we  can  deploy  the  virtual 
cell  system  according  to  the  topology  of  the  optimal  partition  such  that  a  cluster 
corresponds  to  a  virtual  cell.  To  analyze  the  performance  of  a  virtual  cell  system,  we 
adopt  an  open  multiple  class  queueing  network  model.  In  addition  to  mobility  and 
communication  patterns  among  physical  cells,  the  topology  of  the  virtual  cell  system 
is  used  to  determine  service  transition  probabilities  of  the  queueing  network  model. 
By  solving  the  traffic  equations  of  the  queueing  network  model,  we  obtain  various 
performance  measures  such  as  the  network  response  time  for  each  type  of  message 
and  the  utilization  of  the  base  station  networks  and  the  backbone  network  of  the 
virtual  cell  system. 
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CHAPTER  1 
INTRODUCTION 

1.1  Objectives 

The  virtual  cell  system  supports  host  mobility  at  the  data  link  layer  using  high- 
speed base  station  networks  for  the  transmission  of  IP  datagrams  in  mobile  computer 
communications.  Integrating  host  mobility  into  the  IP  layer  in  an  internet  has  been 
a  common  approach  for  the  transmission  of  IP  datagrams  in  TCP/IP  environments. 
No  existing  mobile  computer  communications  networks  support  host  mobility  at  the 
data  link  layer  and  use  high-speed  base  station  networks. 

Supporting  host  mobility  at  the  IP  layer  and  using  an  internet  for  base  station 
networking  may  cause  several  problems  in  terms  of  the  compatibility  with  the  conven- 
tional IP  protocol.  Other  issues  that  need  to  be  addressed  are  the  efficient  location 
tracking  strategies  for  mobile  hosts  and  optimal  network  management  and  control. 
The  objectives  of  this  research  are 

•  to  design  the  communication  architecture  and  protocol  of  the  virtual  cell  system 
to  support  IP-level  mobile  communications, 

•  to  optimally  deploy  the  virtual  cell  system  according  to  mobility  and  commu- 
nication patterns  among  physical  cells  in  both  highway  cellular  systems  and 
hexagonal  cellular  systems,  and 

•  to  evaluate  the  performance  of  the  virtual  cell  system  in  hierarchical  networking 
environments. 
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1.2  Motivations 

A  fixed  host  in  TCP/IP  environments  is  always  assigned  an  IP  address.  The 
IP  address  is  not  only  used  as  a  unique  identifier  by  higher  layer  protocols  but  also 
represents  the  current  location  of  the  host.  An  inherent  problem  from  the  trans- 
mission of  IP  datagrams  in  mobile  computer  communications  is  that  the  IP  address 
of  a  migrating  mobile  host  (MH)  is  only  valid  as  its  identification  information  but 
not  able  to  represent  its  current  location  information.  To  solve  this  problem,  various 
researches  and  developments  have  focused  on  how  to  integrate  the  functionality  of 
host  mobility  into  the  IP  layer,  preserving  the  compatibility  with  the  conventional 
IP  protocol  [1,  2,  3,  4,  5,  6]. 

Although  existing  mobile  host  protocols  vary  on  how  to  represent  and  maintain 
location  information  for  efficient  tracking  of  MHs,  the  techniques  of  using  the  IP 
options  and  IP  encapsulation  have  been  considered.  Typical  examples  of  the  first 
technique  are  Virtual  Internet  Protocol  (VIP)  [2,  3]  and  the  mobile  host  protocol 
using  the  IP  Loose  Source  Routing  option  [4].  IP-within-IP  (IPIP)  [1]  and  Internet 
Packet  Forwarding  Protocol  [6]  use  the  second  technique.  Even  though  the  details 
of  those  mobile  host  protocols  are  different,  we  consider  one  representative  mobile 
host  protocol  for  each  technique,  VIP  and  IPIP,  to  illustrate  common  features  and 
problems  of  integrating  host  mobility  into  the  IP  layer. 

In  VIP,  the  network  layer  is  divided  into  two  layers;  the  VIP  layer  resides  on 
top  of  the  normal  IP  layer.  The  VIP  packet  header  is  implemented  as  an  option 
of  the  IP  packet  header.  An  MH  keeps  its  permanent  IP  address  used  at  the  VIP 
layer  for  its  identification  and  acquires  a  transient  IP  address  used  at  the  IP  layer 
for  its  current  location  information  when  it  migrates.  If  the  MH  is  in  its  native 
network,  both  addresses  are  the  same.  After  acquiring  a  transient  IP  address,  the 
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migrating  MH  sends  its  native  network  a  notification  packet  whose  header  contains 
the  permanent  and  transient  IP  addresses.  As  the  packet  propagates  to  the  native 
network,  every  network  entity  including  MHs,  mobile  support  gateways  (MSGs),  and 
even  intermediate  gateways  on  the  path  snoops  the  header  information  and  stores 
address  conversion  information  to  a  cache.  In  the  same  way,  the  header  information  of 
all  data  packets  in  transit  is  also  used  by  the  network  entities  to  maintain  their  caches. 
If  the  source  has  the  cache  entry  for  the  destination,  the  source  executes  address 
conversion  to  insert  the  correct  location  address  in  the  header.  The  existing  routing 
protocols  of  the  IP  layer  can  then  correctly  deliver  this  packet  to  the  destination. 
Otherwise,  the  source  assumes  that  the  destination  resides  in  its  native  network  and 
sends  the  packet  accordingly.  As  the  packet  traverses  to  the  native  network  of  the 
destination,  if  an  intermediate  gateway  has  the  cache  entry  for  the  destination,  the 
gateway  executes  address  conversion  and  forwards  the  packet  to  the  current  location 
of  the  destination. 

Unlike  VIP,  intermediate  gateways  in  IP  IP  are  not  involved  in  the  support  of 
host  mobility  but  merely  in  transport  service.  The  source  MSG  encapsulates  a  net- 
work control  packet  for  mobility  management  or  an  IP  datagram  from  an  MH  into 
another  IP  datagram  whose  source  and  destination  addresses  specify  the  communi- 
cating MSGs  and  transmits  it  over  an  internet.  The  existing  routing  protocols  of 
the  IP  layer  then  correctly  deliver  it  to  the  destination  MSG.  Thus,  MSGs  consider 
the  internet  as  a  full  mesh  of  logical  point-to-point  links  to  interconnect  them.  A 
mobile  network  consists  of  a  number  of  MSGs,  each  of  which  maintains  the  location 
information  of  its  own  MHs.  Every  MH  in  the  mobile  network  is  assigned  a  unique  IP 
address,  but  the  network  part  of  the  IP  address  is  the  same.  When  the  MH  migrates, 
a  forwarding  pointer  is  set  from  the  previous  MSG  to  the  new  MSG  for  location 
tracking.  If  the  MSG  has  no  location  information  for  a  specific  MH,  it  broadcasts 


4 


the  other  MSGs  in  the  same  mobile  network  an  inquiry  message  asking  who  has  the 
MH.  However,  IPIP  also  requires  a  transient  IP  address  when  the  MH  visits  a  foreign 
mobile  network. 

Regardless  of  which  technique  is  used,  the  integration  of  host  mobility  into  the 
IP  layer  reveals  several  implications.  First,  underlying  networks  differ  widely  in  their 
network  size,  bandwidth,  protocol,  and  packet  size  so  that  they  or  some  of  them 
may  not  meet  performance  needs  for  the  support  of  host  mobility  such  as  rapid 
migration  and  tracking  of  MHs.  Moreover,  because  they  are  usually  under  different 
administrative  controls,  efficient  network  management  and  optimization  for  the  global 
mobile  network  may  not  be  easy  tasks. 

Second,  intermediate  gateways  may  cause  some  performance  problems  no  matter 
if  they  are  involved  in  the  support  of  host  mobility  or  not.  If  they  are  involved,  as  in 
VIP,  they  must  be  modified  or  replaced  to  understand  VIP  so  that  the  benefit  from 
using  an  existing  internet  as  a  transport  network  is  diminished.  Furthermore,  because 
intermediate  gateways  have  to  snoop  every  packet  in  transit  to  maintain  location 
information  and  examine  every  data  packet  in  transit  in  order  to  perform  address 
conversion,  the  protocol  processing  load  at  intermediate  gateways  may  severely  affect 
overall  performance.  Even  in  those  cases  in  which  they  are  not  involved,  as  in  IPIP, 
because  they  usually  implement  packet  switching  operations  in  software,  protocol 
processing  time  coupled  with  possible  network  delay  between  multiple  hops  of  IP 
gateways  may  greatly  affect  message  delivery  latency. 

Third,  each  physical  cell  administered  by  an  MSG  is  in  principle  assigned  an  IP 
network  address  because  every  MH  and  intermediate  gateway  is  a  network  entity 
in  TCP/IP  environments.  It  means  that  an  MH  migrating  to  a  different  physical 
cell  requires  a  different  network  address  to  represent  its  current  location.  Some 
undesirable  features  of  IP-level  mobile  host  protocols  essentially  come  from  this  fact. 
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VIP  requires  a  large  amount  of  the  transient  IP  address  space,  which  becomes  a 
very  scarce  network  resource  because  internets  are  rapidly  growing.  IPIP  uses  one 
permanent  IP  address  for  each  MH  but  relies  on  the  broadcast  inquiry  among  MSGs 
when  the  location  information  is  not  available.  This  restricts  the  scalability  of  IPIP 
to  a  local  area. 

Fourth,  although  IP-level  mobile  host  protocols  based  on  options  keep  the  com- 
patibility with  the  conventional  IP  protocol,  in  a  practical  sense  they  may  interfere 
with  it  because  most  existing  fixed  hosts  and  intermediate  gateways  do  not  imple- 
ment the  IP  options.  The  implementations  to  support  the  IP  options  may  not  be 
feasible  in  the  near  future.  In  addition,  current  efforts  to  provide  IP  multicasting 
protocols  and  connection-oriented  protocols  require  IP-level  mobile  host  protocols  to 
be  compatible. 

As  high-speed,  connectionless,  packet-switched  networks  are  emerging  to  extend 
LAN-like  performance  across  a  wide  area,  we  believe  that  they  can  greatly  affect 
the  support  of  host  mobihty  in  mobile  computer  communications.  Examples  are 
ATM  networks  [7,  8]  and  Switched  Multimegabit  Data  Services  (SMDS)  networks 
[9,  10],  which  are  subnetworks  providing  a  MAC  service  across  a  wide  area  in  a 
large  interconnected  network.  To  take  advantage  of  these  networks,  it  is  desirable  to 
estabhsh  the  support  for  host  mobility  at  the  data  link  layer  between  base  stations 
(BSs). 

1.3    Design  of  the  Virtual  Cell  System 

In  this  research,  we  design  the  network  infrastructure  of  the  virtual  cell  system  for 
the  transmission  of  IP  datagrams  in  mobile  computer  communications.  The  virtual 
cell  system  shields  host  mobility  from  the  IP  layer  by  supporting  it  at  the  data  Hnk 
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layer  using  base  station  networking.  To  design  the  virtual  cell  system,  there  are 
several  issues  to  be  addressed: 

•  analyzing  the  requirements  of  base  station  networks, 

•  identifying  the  current  location  of  MHs  without  using  their  IP  addresses, 

•  constructing  a  distributed  location  information  structure  for  efficient  location 
tracking, 

•  reflecting  the  impact  of  host  mobility  on  the  distributed  location  information 
in  the  design  of  protocols, 

•  designing  the  communication  architecture  and  protocol,  and 

•  harmonizing  the  communication  architecture  and  protocol  with  the  conven- 
tional IP  protocol. 

Considering  the  requirements  of  base  station  networks  from  two  different  view- 
points, application  and  mobility  management,  the  virtual  cell  system  takes  advantage 
of  high-speed  datagram  packet-switched  networks  with  the  multicast  ability  for  base 
station  networking.  The  base  station  networks  interconnect  remote  bridge  BSs,  so  as 
to  preserve  the  interconnection  level  of  physical  cells  at  the  data  link  layer.  The  cur- 
rent location  of  MHs  is  identified  by  the  base  station  network  address.  It  represents 
which  BS  can  communicate  with  a  particular  MH.  For  efficient  location  tracking  of 
MHs,  the  virtual  cell  system  constructs  a  distributed  hierarchical  location  information 
structure.  Based  on  the  distributed  hierarchical  location  information,  the  virtual  cell 
protocol  for  handolF,  address  resolution,  and  data  transfer  is  designed.  The  handoff 
procedure  can  utifize  the  multicast  ability  of  the  base  station  network  to  achieve  the 
consistency  of  the  distributed  location  information.  Since  the  IP  network  address  of 
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a  migrating  MH  represents  at  least  the  near  location  information  of  the  MH  in  the 
virtual  cell  system,  the  distributed  location  information  coupled  with  the  existing 
routing  protocols  of  the  IP  layer  give  the  same  effect  of  maintaining  a  conceptually 
centralized  server  for  the  whole  mobile  network. 

1.4    Optimal  Deployment  of  the  Virtual  Cell  System 

A  virtual  cell  should  accommodate  a  limited  number  of  BSs  in  order  to  produce  an 
acceptable  performance  no  matter  how  the  network  parameters  such  as  the  topology 
and  capacity  of  the  virtual  cell  system  are  engineered.  In  addition,  the  cost  of  tracking 
MHs  within  a  virtual  cell  is  lower  than  between  virtual  cells.  Thus,  in  order  to 
decrease  the  total  communication  cost  of  tracking  MHs,  it  is  intuitively  desirable  that 
the  two  BSs  between  which  the  MH  migrates  or  communicates  very  often  should  be 
contained  in  the  same  virtual  cell,  while  the  two  BSs  between  which  the  MH  migrates 
or  communicates  infrequently  could  be  separated  into  different  virtual  cells. 

The  location  tracking  of  MHs  involves  two  primitive  operations  in  the  virtual 
cell:  the  move  operation  to  perform  handoff  and  the  find  operation  to  locate  the 
current  BS  of  an  MH  for  the  transmission  of  IP  datagrams.  Then  mobility  and  data 
communication  patterns  among  BSs  can  be  represented  by  the  frequencies  of  move 
and  find  operations  among  BSs.  Even  though  there  is  a  tradeoff  between  the  costs  of 
move  and  find  operations  within  a  virtual  cell,  the  optimal  deployment  of  the  virtual 
cell  system  is  concerned  with  optimal  partitioning  of  BSs  into  virtual  cells,  so  as  to 
minimize  the  total  communication  cost  for  a  sequence  of  move  and  find  operations 
in  the  entire  system. 

In  addition  to  mobihty  and  communication  patterns  among  BSs,  the  optimal 
deployment  of  the  virtual  cell  system  also  considers  the  physical  arrangement  of  BSs 
so  that  BSs  in  a  virtual  cell  are  contiguously  connected  in  their  underlying  topology. 
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not  scattered.  While  the  frequency  is  represented  by  a  complete  directed  graph 
among  BSs  for  a  possible  communication  between  MHs  through  different  BSs,  the 
physical  topology  is  represented  by  a  linear  graph  in  highway  cellular  systems  or  by 
a  hexagonal  mesh  graph  in  cellular  systems.  Hence,  the  optimal  deployment  of  the 
virtual  cell  system  is  concerned  with  the  two  types  of  graphs,  the  topology  graph  for 
partitioning  BSs  and  the  frequency  graph  for  optimizing  communication  cost. 

For  general  graphs,  the  problem  of  finding  a  cover  of  disjoint  sets  such  that  the 
sum  of  all  edge  weights,  whose  two  endpoints  are  in  two  diiferent  sets,  is  equal  to  or 
less  than  a  given  positive  integer  is  known  as  NP-complete  for  the  arbitrary  size  of  a 
set  [13].  The  k-cut  problem  of  finding  a  partition  of  vertices  into  k  nonempty  clusters 
such  that  the  total  edge  weight  between  clusters  is  minimum  is  also  NP-complete  for 
arbitrary  k.  The  A:-cut  problem  with  specified  vertices,  which  is  an  extension  of  the 
2-cut  problem  solvable  via  repeated  applications  of  a  max-flow  min-cut  algorithm, 
becomes  NP-hard  even  ioi  k  —  3  [14].  Given  a  planar  graph  with  an  s-outerplanar 
embedding,  an  optimal  solution  for  maximum  independent  set  can  be  obtained  in 
time  0(8' n)  by  a  dynamic  programming  technique,  where  n  is  the  number  of  nodes 
[15].  If  an  n- vertex  planar  graph  G  is  given  with  an  5-outerplane-separable  planar 
embedding  of  G,  then  the  optimal  partition  of  two  clusters  of  fixed  size  is  determined 
in  time  0{s^n^2^')  by  a  dynamic  programming  technique  [16].  This  implies  that 
without  considering  the  frequency  graph,  the  2-way  partition  in  a  hexagonal  mesh  of 
cellular  systems  is  exponential  where  s  =  0(n^/^). 

Given  the  frequencies  of  move  and  find  operations  among  n  BSs  in  highway  cellu- 
lar systems,  we  develop  an  efficient  optimal  partitioning  algorithm  of  0{mn^)  for  an 
arbitrary  number  of  clusters  m  by  dynamic  programming.  In  a  hexagonal  mesh  of 
physical  cells,  we  develop  several  heuristics  for  multiway  partitioning,  based  on  the 
techniques  of  moving  or  interchanging  the  boundary  nodes  between  adjacent  clusters. 
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These  heuristics  produce  optimal  partitions  with  respect  to  the  initial  partition  ob- 
tained randomly  or  by  centering.  The  heuristics  are  compared  and  shown  to  behave 
quite  well  through  experimental  tests  and  analysis. 

1.5    Performance  Analysis  of  the  Virtual  Cell  System 

Once  an  optimal  partition  of  m  disjoint  clusters  is  obtained,  the  virtual  cell  system 
can  be  deployed  so  that  each  cluster  corresponds  to  a  virtual  cell.  Virtual  cell  i,  where 
1  <  i  <  m,  is  constructed  by  interconnecting  the  BSs  of  the  ith.  cluster  and  a  location 
server  by  a  base  station  network.  These  m  virtual  cells  are  interconnected  by  the 
backbone  network. 

To  evaluate  the  performance  of  virtual  cells,  we  apply  an  open  multiple  class 
queueing  network  model.  There  are  three  types  of  messages  entering  or  leaving  the 
virtual  cell  system  via  BSs:  the  handoff  message,  the  data  message,  and  the  address 
resolution  message.  The  handoff  messages  are  generated  due  to  move  operations,  and 
the  data  and  address  resolution  messages  are  due  to  find  operations.  The  move  and 
find  frequencies  in  conjunction  with  the  topology  matrix  of  the  virtual  cell  system 
obtained  by  the  optimal  algorithms  are  used  to  determine  service  transition  proba- 
bilities for  each  type  of  message  in  the  queueing  network  model. 

With  various  system  parameters,  we  can  conduct  interesting  sensitivity  analyses 
to  determine  network  design  tradeoffs.  The  first  application  of  the  proposed  model 
is  to  determine  an  adequate  network  bandwidth  for  base  station  networks  such  that 
the  networks  would  not  become  a  bottleneck.  Also,  we  can  evaluate  the  network 
utilization  due  to  various  messages.  For  instance,  when  the  vehicles  begin  moving 
fast,  the  migration  rate  will  be  increased.  This  implies  more  handoff  messages  and 
more  forwarding  operations.   The  network  traffic  should  be  increased  accordingly. 
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The  performance  analysis  results  should  provide  a  good  evidence  to  demonstrate  the 
system  efficiency  under  different  mobility  assumptions. 


CHAPTER  2 
SURVEY  OF  RELEVANT  WORK 


In  this  chapter,  we  survey  three  representative  approaches  in  mobile  data  commu- 
nications: Cellular  Packet  Switch  (CPS)  [11,  12],  IP-within-IP  (IPIP),  and  Virtual 
Internet  Protocol  (VIP).  All  three  approaches  commonly  take  advantage  of  packet- 
switching  technology  to  keep  track  of  mobile  users.  However,  CPS  is  designed  for 
digitized  voice  and  data  services  in  personal  communications  and  handles  virtual  cir- 
cuit connection  at  the  data  link  layer  for  mobility  management.  IPIP  and  VIP  are 
designed  for  TCP/IP  environments  and  integrate  host  mobility  into  the  IP  layer. 

2.1    Cellular  Packet  Switch 

The  increasing  demand  in  wireless  communications  has  propelled  the  development 
of  new  and  improved  cellular  and  cordless  systems.  Even  with  these  systems,  how- 
ever, increased  voice  traffic  together  with  the  ever  growing  demand  for  data  services 
have  revealed  capacity  problems.  To  provide  the  necessary  capacity  in  a  densely  pop- 
ulated area,  more  efficient  use  of  the  limited  radio  spectrum  will  be  required,  and  a 
closely  spaced  grid  of  microcells  will  have  to  be  deployed.  While  small  cells  mitigate 
capacity  problems,  the  reduced  cell  size  makes  the  task  of  locating  and  controlling  the 
mobile  users  substantially  more  complex.  CPS  is  intended  to  support  high  density 
personal  communication  networks  using  a  new  network  architecture  based  on  packet 
switching.  A  packet  switched  architecture  itself  enables  distributed  control  of  mobil- 
ity management  and  radio  resource  management  functions  because  the  addresses  and 
other  information  in  packet  headers  make  it  possible  for  dispersed  network  elements 
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HDB:  Home  DataBase 

VDB:  Visitor  DataBase 

BIU:  Base  statioa  Interface  Unit 

VIU:  Visitor  database  Interface  Unit 

nU:  Trunk  Interface  Unit 

CIU:  Cellular  controller  Interface  Unit 

GIU:  Gateway  Interface  Unit 


Figure  2.1.  Cellular  Packet  Switch  Architecture 

to  respond  to  user  mobility  without  the  intervention  of  central  controllers.  In  addi- 
tion, it  is  efficient  in  integrating  the  diverse  types  of  information  that  will  be  offered 
to  future  networks. 

2.1.1    Network  Architecture 

As  shown  in  Figure  2.1,  CPS  takes  advantage  of  a  Metropolitan  Area  Network 
(MAN)  based  on  the  IEEE  802.6  standard  in  order  to  support  the  increased  traf- 
fic resulting  from  the  use  of  microcells.  The  MAN  connects  base  stations,  control 
units,  databases,  and  links  to  fixed  information  services,  such  as  the  Public  Switched 
Telephone  Network.  In  addition  to  routing,  the  interface  units  are  responsible  for 
information  transfer  and  protocol  conversion  through  different  media  in  the  network 
fabric  including  the  MAN,  the  fixed  network,  and  the  radio  transmission  systems.  In 
CPS,  call  processing  (setup  and  release)  is  handled  by  the  central  office  switch,  while 
mobility  management  and  radio  resource  management  are  handled  on  the  MAN.  In- 
formation on  the  MAN  is  organized  in  data  units  compatible  with  ATM,  and  data 
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units  are  routed  to  their  correct  destinations  according  to  addresses  in  packet  head- 
ers. The  base  station  interface  units  translate  information  between  the  MAN  format 
and  that  of  the  radio  link.  The  trunk  interface  units  convert  user  information  be- 
tween the  MAN  protocol  and  the  fixed  network  protocol,  which  may  be  ISDN  or 
broadband  ISDN,  for  example.  The  home  database  interface  unit  could  connect  the 
MAN  directly  to  a  database  that  serves  CPS.  It  also  could  provide  connectivity  to  a 
general  purpose  data  network  connected  to  home  databases.  Other  interface  units, 
such  as  the  cellular  controller  interface  unit  and  the  visitor  database  interface  unit, 
simply  provide  their  network  elements  with  access  to  the  MAN.  Because  the  coverage 
area  of  a  single  MAN  could  be  as  low  as  a  few  square  kilometers  in  densely  populated 
areas,  the  gateway  interface  unit  interconnects  adjacent  MANs.  It  permits  CPS  to 
perform  handoffs  between  base  stations  served  by  different  MANs. 

2.1.2    CPS  Protocols 

CPS  protocols  distribute  network  control  functions  among  the  interface  units 
using  virtual  circuit  packet  switching.  The  link  between  a  wireless  terminal  and  the 
fixed  network  is  a  virtual  circuit  identifier  (VCI),  which  is  a  number  carried  in  the 
header  of  each  packet  on  the  MAN.  By  referring  to  this  VCI,  MAN  interface  units 
route  each  packet  to  its  correct  destination.  Thus,  this  VCI  establishes  a  logical  link 
between  interface  units.  At  the  beginning  of  a  call,  the  cellular  controller  interface 
unit  assigns  a  call-specific  VCI  to  the  call.  The  base  station  interface  unit  that 
initially  handles  the  call  and  the  trunk  interface  unit  that  moves  the  call  to  and  from 
the  fixed  network  both  record  the  VCI.  When  the  call  moves  to  a  new  base  station, 
that  base  station  records  the  VCI.  It  also  sends  a  message  to  the  previous  base  station 
telling  that  base  station  to  delete  the  VCI  from  its  table. 


14 


The  following  example  shows  a  wireless  terminal  receiving  information  from  a 
trunk  interface  unit  when  it  moves  from  the  service  area  of  base  station  BSi  to  that 
of  base  station  BS2.  Initially,  the  terminal  is  communicating  with  and  has  a 
call  assigned  VCI  100  by  the  cellular  controller  interface  unit.  First  the  terminal 
determines  that  the  radio  link  to  is  inadequate  and  that  an  acceptable  link  is 
available  to  BS2-  It  spontaneously  sends  a  handolf  message  to  BS2-  This  message 
contains  VCI  100  and  the  identity  of  ^^i.  BS2  enters  VCI  100  in  its  address  table 
and  sends  a  message  to  BSi,  which  causes  to  delete  VCI  100  from  its  table.  The 
handoff  message  from  the  terminal  to  BS2  contains  the  sequence  number  of  the  last 
nonvoice  packet  received  accurately  through  BS^.  It  relays  this  information  to  the 
trunk  interface  unit,  which  transmits  to  BS2  nonvoice  packets  that  it  previously  sent 
to  BSi  while  the  handoff  was  in  progress.  Thus,  all  further  transmissions  from  the 
trunk  interface  unit  go  through  BS2. 

The  handoff  is  terminal  initiated,  which  means  that  the  terminal  makes  the  de- 
cision that  a  handoff  is  necessary.  It  is  also  terminal  controlled,  which  means  that 
the  terminal  decides  to  move  the  call  to  another  base  station.  In  the  overlapped  area 
between  two  different  radio  cells,  the  terminal  communicates  with  either  base  station, 
but  not  with  both  simultaneously.  The  network  stores,  for  possible  retransmission, 
critical  information  (data  packets)  in  a  source  buffer  at  the  trunk  interface  unit.  It 
could  also  store  this  information  in  a  destination  buffer  in  the  base  interface  unit. 

In  addition  to  this  handoff  procedure,  CPS  supports  other  protocols  that  perform 
a  variety  of  call  processing  and  mobility  management  functions,  including  terminal 
initiated  call  setup,  network  initiated  call  setup,  terminal  initiated  call  release,  net- 
work initiated  call  release,  location  updating  and  authentification,  inter-MAN  hand- 
off,  and  path  optimization.  At  the  wireless  terminal  and  in  the  fixed  network,  the 
call  setup  and  release  procedures  are  the  Q.931  protocols  of  ISDN.  Within  CPS  these 
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protocols  are  augmented  by  additional  functions  necessary  to  enter  or  delete  a  VCI 
from  the  address  tables  in  the  WIU,  the  initial  base  station  interface  unit,  and  the 
trunk  interface  unit.  For  location  tracking,  protocols  developed  for  existing  cellular 
systems  are  used.  Meier-Hellstern  et  al.  [12]  describes  wireless  services  based  on  the 
IEEE  802.6  standards  in  detail. 

2.2    Virtual  Internet  Protocol 

Teraoka,  et  al.  [2,  3]  identify  the  need  of  host  migration  transparency  to  provide 
the  same  computational  environment  with  mobile  users  in  a  large  interconnected 
network.  The  main  goal  is  to  support  host  portability,  which  implies  no  communica- 
tions during  the  migration  of  mobile  hosts.  To  achieve  host  migration  transparency, 
the  virtual  network  concept  that  uses  the  propagating  cache  method  has  been  pro- 
posed. As  an  example  of  the  virtual  network,  VIP  is  derived  from  the  conventional 
IP  protocol. 

2.2.1    Network  Architecture 

In  order  to  provide  a  transport  layer  with  host  migration  transparency,  two  net- 
work layer  identifiers  of  a  host  are  introduced:  one  is  migration  independent  and 
the  other  is  migration  dependent.  The  transport  layer  specifies  the  target  host  by 
the  migration  independent  identifier  so  that  the  transport  layer  is  not  aware  of  host 
migration.  To  support  this  migration  independent  identifier,  the  concept  of  virtual 
network  is  introduced.  The  virtual  network  is  a  logical  network.  Each  host  is  consid- 
ered to  be  permanently  connected  to  a  particular  virtual  network,  called  the  home 
network  of  the  host.  A  host  never  migrates  among  virtual  networks  even  if  it  migrates 
among  the  physical  networks. 

From  the  prospective  of  the  network  protocol  architecture,  a  conventional  net- 
work layer  is  divided  into  two  sublayers:  virtual  network  sublayer  (VN-sublayer)  and 
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physical  network  sublayer  (PN-sublayer).  A  virtual  network  address  (VN-address)  is 
assigned  to  a  host  in  the  VN-sublayer,  and  a  physical  network  address  (PN-address) 
is  assigned  to  a  host  in  the  PN-sublayer.  The  VN-address  of  a  host  never  changes  no 
matter  where  the  host  migrates,  while  the  PN-address  changes  if  the  host  migrates. 
Thus,  the  VN-address  of  a  host  indicates  its  home  network,  to  which  it  was  initially 
connected. 

2.2.2  Propagating  Cache  Method 

Since  the  transport  layer  specifies  the  target  host  by  its  VN-address,  the  main 
function  of  the  VN-sublayer  is  to  convert  a  VN-address  into  its  corresponding  PN- 
address.  To  solve  this  problem,  the  propagating  cache  method  is  proposed.  In  the 
propagating  cache  method,  each  host  and  gateway  has  a  cache  for  address  conversion, 
called  an  Address  Mapping  Table  (AMT).  If  the  source  host  has  an  AMT  entry  for 
the  destination  host,  the  source  host  executes  address  conversion  before  sending  out  a 
packet;  the  PN-sublayer  can  then  correctly  deliver  this  packet  to  the  destination  host. 
Otherwise,  the  source  host  assumes  that  the  destination  host  is  connected  to  its  home 
network  and  sends  the  packet  accordingly.  As  the  packet  traverses  the  network  to  the 
home  network  of  the  destination  host,  if  an  intermediate  gateway  has  the  AMT  entry 
for  the  destination  host,  the  gateway  executes  the  address  conversion  and  forwards 
the  packet  to  the  current  physical  location  of  the  destination  host.  Whenever  a  host 
or  gateway  receives  a  packet,  the  VN-sublayer  creates  or  updates  the  AMT  entry  for 
the  source  host  of  the  received  packet.  Thus,  the  AMT  information  propagates  across 
the  interconnected  networks  as  communication  progresses. 

2.2.3  Communication  Procedures 

VIP  is  derived  by  applying  the  virtual  network  concept  to  the  IP  layer.  The 
conventional  IP  layer  is  divided  into  two  sublayers:  the  VIP  sublayer  and  the  IP 
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sublayer.  A  host  has  an  IP  address  in  the  IP  sublayer  and  a  VIP  address  in  the  VIP 
sublayer.  Since  the  VIP  address  of  a  host  never  changes  even  if  the  host  migrates 
to  another  network,  the  TCP/UDP  layer  and  its  application  programs  can  uniquely 
specify  a  host  by  the  VIP  address.  The  IP  address  becomes  invisible  to  the  TCP/UDP 
layer.  Since  the  IP  address  is  location  dependent,  however,  it  changes  when  a  host 
migrates  to  another  network.  This  means  that  the  host  must  be  assigned  a  new  IP 
address.  The  function  to  allocate  an  IP  address  to  a  host  is  not  included  in  VIP  and 
implemented  with  daiemon  processes. 

In  packet  transmission,  the  TCP/UDP  layer  issues  a  transmission  request  to  the 
VIP  sublayer  with  the  VIP  address  of  the  destination  host.  The  VIP  sublayer  then 
generates  the  VIP  header  and  sets  each  field.  Next,  the  VIP  sublayer  searches  for 
the  AMT  entry  for  the  destination  host.  An  AMT  entry  consists  of  four  data  fields: 
VIP  address  used  as  key,  IP  address  of  the  host,  address  timestamp  to  specify  the 
version  number  of  the  VIP/IP  address  pair  of  this  entry,  and  idle  timer  used  for  aging 
the  AMT  entry.  If  the  entry  is  found,  the  VIP  sublayer  converts  the  VIP  address  of 
the  destination  host  into  the  corresponding  IP  address.  Otherwise,  the  VIP  sublayer 
assumes  that  the  VIP  address  and  the  IP  address  of  the  destination  host  are  the 
same.  This  means  that  the  VIP  sublayer  assumes  that  the  destination  host  exists 
in  its  home  network.  The  procedures  in  the  IP  sublayer  are  the  same  as  that  of  the 
conventional  IP  layer. 

In  packet  reception,  some  modifications  to  the  existing  IP  procedures  are  neces- 
sary. When  the  data  fink  layer  notifies  the  IP  sublayer  of  packet  reception,  the  IP 
sublayer  calls  a  function  in  the  VIP  sublayer  that  creates  or  updates  the  AMT  entry 
for  the  source  host  of  the  received  packet.  If  the  AMT  entry  for  the  source  host  does 
not  exist,  an  entry  is  created.  If  the  entry  already  exists  and  either  the  IP  address  of 
the  source  host  has  changed  or  the  value  of  the  timestamp  field  of  the  received  packet 
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is  newer  than  that  of  the  timestamp  field  in  the  AMT  entry  for  the  source  host,  then 
the  entry  is  updated. 

After  this  modification  to  the  AMT  entry,  the  IP  sublayer  compares  the  destina- 
tion IP  address  field  of  the  received  packet  with  its  own  IP  address.  If  the  two  IP 
addresses  match,  the  IP  sublayer  passes  the  packet  to  the  VIP  sublayer.  Otherwise, 
the  IP  sublayer  calls  another  function  in  the  VIP  sublayer  in  order  to  convert  the 
VIP  address  of  the  destination  host  into  the  corresponding  IP  address,  if  necessary. 
If  the  AMT  entry  for  the  destination  host  is  found  and  the  value  of  the  timestamp 
field  of  the  received  packet  is  older  than  that  of  the  timestamp  field  in  the  AMT 
entry  for  the  destination  host,  this  function  returns  the  values  of  the  IP  address  field 
and  the  timestamp  field.  If  the  address  conversion  occurs,  the  IP  sublayer  modifies 
the  destination  IP  address  field  in  the  IP  header  and  the  VIP  sublayer  modifies  the 
timestamp  field  in  the  VIP  header.  The  packet  is  then  forwarded  to  the  next  gateway 
or  to  the  destination  host. 

2.2.4    Migration  Procedures 

When  a  migrating  host  is  connected  to  a  subnetwork,  the  host  sends  a  connection 
notification  packet  to  its  home  network  and  waits  for  a  connection  acknowledgement 
packet.  The  connection  notification  packet  announces  that  the  source  host  has  just 
been  attached  to  a  subnetwork  and  is  then  relayed  by  gateways.  Since  this  connection 
notification  packet  includes  the  VIP  and  IP  addresses  of  the  source  host,  intermediate 
gateways  can  learn  the  relation  of  these  two  addresses  and  create  or  update  the  AMT 
entry  for  the  source  host.  The  home  gateway,  which  is  the  gateway  in  the  home 
network  of  the  migrating  host,  receives  the  connection  notification  packet,  creates 
the  AMT  entry  for  the  migrating  host,  and  returns  the  connection  acknowledgement 
packet  to  it.  Note  that  the  home  gateway  never  deletes  this  AMT  entry  by  timeout. 
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The  home  gateway  also  broadcasts  the  connection  notification  packet  within  the  home 
network  of  the  migrating  host  if  the  home  network  is  a  broadcast  type  network  such 
as  Ethernet.  In  addition,  the  home  gateway  creates  an  entry  in  the  ARP  table  for 
the  migrating  host  to  answer  the  ARP  request  for  that  host. 

Upon  packet  transmission  to  a  migrating  host,  if  the  source  host  does  not  have 
the  AMT  entry  for  the  migrating  host,  the  host  assumes  that  the  IP  address  of  the 
migrating  host  is  equal  to  its  VIP  address.  The  packet  then  heads  for  the  home 
network  of  the  migrating  host.  If  a  gateway  on  the  path  has  the  AMT  entry  for 
the  migrating  host,  the  gateway  converts  the  VIP  address  into  the  corresponding  IP 
address  and  forwards  the  packet  to  the  actual  location  of  the  migrating  host.  If  the 
migrating  host  returns  a  packet  to  the  source  host,  the  source  host  can  create  the 
AMT  entry  for  the  migrating  host  so  that  the  host  can  thereafter  directly  specify  the 
IP  address  of  the  migrating  host. 

When  a  migrating  host  is  about  to  disconnect  from  a  subnetwork,  it  sends  a 
disconnection  notification  packet  to  its  home  gateway.  When  the  home  gateway 
receives  that  packet,  it  broadcasts  a  control  packet  to  all  subnetworks  to  which  it 
is  connected  so  that  a  relevant  AMT  entry  is  deleted.  Gateways  in  the  subnetwork 
receive  the  control  packet.  If  a  gateway  has  the  AMT  entry  for  the  migrating  host, 
it  deletes  the  entry  and  also  broadcasts  the  control  packet  to  any  other  connected 
subnetworks.  If  a  gateway  does  not  have  an  AMT  entry  for  the  migrating  host,  the 
gateway  does  nothing.  Thus,  most  AMT  entries  of  the  migrating  host  are  deleted 
when  the  host  is  about  to  disconnect  from  a  subnetwork. 

2.3  IP-Within-IP 

IPIP  is  designed  to  provide  continuous  network  access  to  mobile  hosts  in  a  campus 
environment.  The  key  idea  is  to  assign  Mobile  Support  Stations  (MSSs)  to  every  radio 
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Figure  2.2.  A  Model  for  Typical  Campus  Environments 

cell  on  which  mobile  internetworking  protocols  are  supported,  and  attach  them  to  an 
interconnected  network.  Figure  2.2  shows  a  typical  campus  environment  used  in  this 
work. 

2.3.1    Network  Architecture 

All  MHs  have  the  same  internet  network  address  so  that  they  can  be  immediately 
distinguished  from  fixed  hosts  by  their  IP  addresses.  In  addition,  an  MH  has  a 
unique  IP  address  regardless  of  its  migration.  This  eliminates  the  need  of  a  transient 
IP  address,  the  change  of  the  IP  address  at  name  servers,  and  the  notification  of  it  to 
higher  protocols.  Thus,  MSSs  are  responsible  for  the  addressability  of  all  MHs  on  the 
network.  An  MSS  has  its  regular  IP  address  on  the  wired  network,  a  mobile  IP  address 
for  its  wireless  interface,  and  a  second  mobile  IP  address  for  its  wired  interface.  The 
radio  cell  is  the  geographical  area  around  the  MSS  and  is  defined  by  the  propagation 
characteristics  of  the  electromagnetic  waves  at  the  operating  frequency. 

As  a  requirement,  link-level  broadcasts  are  assumed  in  order  for  the  MSS's  beacon 
to  be  received  from  all  the  MHs  in  a  radio  cell.  The  ability  of  the  hardware  to  detect 
signal  strength  and  report  it  to  the  device  driver  is  also  required  to  realize  that 
the  MH  is  at  the  outer  reaches  of  a  radio  cell  and  that  it  should  try  to  switch  to 
a  different  one.  In  addition,  in  the  case  of  spread-spectrum  or  multi-channel  radio 
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communications,  the  ability  to  receive  on  multiple  channels  at  once  is  assumed  so  as 
to  detect  areas  where  radio  cells  overlap  and  decide  when  MHs  switch  to  a  new  MSS 
whose  signal  is  stronger  than  another  MSSs. 

2.3.2    Data  Link  Layer 

There  are  two  kinds  of  interactions  that  take  place  in  the  data-link  layer:  mapping 
logical  addresses  to  physical  addresses  and  dynamically  acquiring  IP  addresses  when 
moving  to  a  foreign  campus.  For  address  mapping,  we  have  the  following  problem 
in  mobile  networking.  When  an  MH  wants  to  communicate  with  another  MH,  and 
the  two  MHs  have  IP  addresses  in  the  same  subnet,  the  first  MH  will  send  an  ARP 
request  for  the  physical  address  of  the  second,  since  it  thinks  the  latter  is  on  the  same 
connected  network.  If  the  second  MH  is  in  the  same  radio  cell,  then  it  will  respond 
to  the  ARP  packet.  If  it  is  not,  the  MSS  will  execute  the  proxy  ARP  protocol  for 
that  other  MH.  When  that  happens,  the  first  MH  will  send  all  traffic  destined  for 
that  other  MH  through  its  MSS,  which  will  now  encapsulate  it  in  an  IPIP  datagram, 
and  deliver  it  to  the  remote  MSS  handling  the  second  MH's  traffic. 

When  an  MH  moves  into  a  foreign  campus,  it  requests  a  transient  IP  address, 
which  it  will  use  to  talk  to  the  local  MSS  and  also  to  encapsulate  IPIP  packets  in 
order  to  maintain  open  connections  to  its  home  campus.  The  MH  knows  it  moved 
into  a  foreign  campus  if  the  network  portion  of  the  MSS's  IP  address  is  different 
from  its  own.  This  can  be  achieved  by  extending  a  number  of  protocols  such  as  the 
Reverse  Address  Resolution  Protocol  and  the  Bootstrap  Protocol  that  provide  such 
physical-to-logical  address  mappings  as  part  of  their  functionality  or  by  using  existing 
and  under-development  protocols  such  as  the  Dynamic  Host  Configuration  Protocol. 
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2.3.3    Network  Layer 

The  IP  layer  of  the  MSS  must  be  modified  to  decide  how  to  route  a  datagram. 
An  entry  is  kept  in  the  kernel  for  each  MH  the  MSS  knows  about.  The  key  field 
in  each  entry  is  the  MH's  home  address.  Additional  fields  include,  for  local  MHs, 
the  IP  address  to  be  used  when  communicating  the  MH  if  it  is  not  the  same  as  the 
home  address;  for  remote  MHs,  the  IP  address  of  the  MSS  handling  the  MH's  traffic; 
a  timer  field,  which  is  reset  every  time  the  entry  is  accessed  and  is  used  to  expire 
entries  that  have  not  been  used  recently.  There  are  two  IP-level  protocols:  IP-within- 
IP  (IPIP)  used  for  the  encapsulation  and  Mobile  Internetworking  Control  Protocol 
(MICP)  used  to  locate  the  other  MSS. 

The  IPIP  protocol  is  used  to  tunnel  IP  datagrams  from  one  cell  to  another  or 
from  one  network  to  another,  essentially  using  IP  for  virtual  point-to-point  connec- 
tions. The  problem  is  how  to  deliver  an  IP  datagram  to  an  MH  whose  location 
cannot  be  deduced  from  its  IP  address.  MSSs  keep  track  of  MHs  and  tunnel  packets 
from  MSS  to  MSS  using  IPIP  encapsulation,  which  is  equivalent  to  IP  loose  source 
routing.  However,  they  differ  in  that  source  routing  depends  on  an  IP  option  being 
implemented  on  all  routers  along  a  path,  whereas  encapsulation  only  depends  on  a 
protocol  module  being  implemented  at  the  two  endpoints  of  the  tunnel. 

When  an  IP  datagram  is  to  be  sent  out  at  the  IP  layer,  the  routing  and  MH 
tables  are  consulted.  If  the  destination  MH  is  in  another  cell  with  a  known  MSS,  it  is 
encapsulated  in  an  IP  datagram  and  sent  to  the  remote  MSS.  Thus,  the  IP  destination 
and  source  addresses  of  the  datagram  can  be  considered  the  two  endpoints  of  the 
tunnel.  Upon  receipt  by  the  remote  MSS,  the  IP  code  hands  the  datagram  to  the 
IPIP  protocol  module.  The  latter  strips  the  IPIP  header  and,  after  decrementing  the 
time-to-live  field,  feeds  the  packet,  which  now  has  the  real  destination  and  source  IP 
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addresses,  back  into  the  IP  queue  so  that  it  can  be  dehvered.  In  addition,  the  code 
checks  whether  the  source  of  the  IP  IP  packet  was  another  MSS. 

The  MICP  is  used  by  the  MSSs  to  exchange  control  information  and  by  the 
MHs  to  signal  to  their  MSSs  that  they  have  changed  cells.  MICP  datagrams  are  also 
encapsulated  in  IP  datagrams  with  a  different  protocol  type  number.  There  are  three 
classes  of  MICP  traffic:  MH-MSS  handshakes,  MSS-MSS  exchanges  to  distribute  MH 
location  information,  and  MSS  configuration  exchanges. 

The  MSS  periodically  broadcasts  its  identity.  An  MH  may  receive  beacons  from 
more  than  one  MSS  at  a  time.  This  is  usually  the  case  where  adjacent  cells  of  radio 
networks  overlap.  If  the  hardware  gives  an  indication  of  the  signal  strength,  and  the 
MH  is  within  the  overlapping  area  of  two  cells,  the  MH  may  want  to  switch  to  the 
cell  that  has  the  strongest  signal.  Otherwise,  the  MH  may  wait  until  it  no  longer 
receives  the  beacon  from  the  previous  MSS  before  it  switches  to  the  new  one. 

The  beacon  is  an  MICP  packet  that  contains  the  IP  address  of  the  MSS's  primary 
interface  and  the  subnet  mask  of  the  current  cell.  It  is  broadcast  periodically  on  the 
local  broadcast  address  of  the  cell.  When  an  MH  receives  a  beacon  packet  that  is 
different  from  the  previous  one  it  received,  it  responds  with  a  greeting  MICP  packet. 
This  packet  contains  the  home  address  of  the  MH,  the  IP  address  of  the  MSS  of  the 
previous  cell  it  was  in. 

When  the  MSS  receives  a  greeting  packet,  it  responds  with  a  greeting  acknowl- 
edgement MICP  packet.  It  updates  the  relevant  system  structures  to  indicate  that 
the  MH  is  now  local  and  records  the  information  supplied  by  the  MH.  Upon  receipt 
of  the  greeting  acknowledgement  packet,  the  MH  changes  its  routing  tables  to  route 
all  packets  through  the  new  MSS  and  now  considers  it  its  new  MSS. 

Beacon  packets  should  be  broadcast  often  enough  that  missing  some  will  not  cause 
the  MH  to  assume  it  is  no  longer  in  the  cell.  If  the  greeting  packet  is  lost,  the  MSS  will 
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not  respond  to  it,  and  the  next  beacon  packet  will  cause  the  MH  to  send  a  new  greet- 
ing. If  the  acknowledgement  packet  is  lost,  the  MH  will  send  another  greeting  until 
it  gets  an  acknowledgement.  After  the  first  successful  beacon/greeting/acknowledge 
handshake,  duplicate  beacons  are  ignored.  Duplicate  greetings  cause  the  MSS  to 
overwrite  the  MH  information,  and  duplicate  acknowledgements  have  no  effect  since 
the  MH  can  tell  they  are  coming  from  its  current  cell. 

When  an  MSS  receives  the  greeting,  it  also  sends  a  forwarding  pointer  MICP 
packet  to  the  MSS  previously  handling  that  MH.  The  previous  MSS  uses  this  infor- 
mation to  handle  properly  any  traffic  for  that  MH  originating  in  the  segments  of  the 
network  it  is  on.  The  MSS  will  keep  retransmitting  the  forwarding  pointer  until  it 
gets  a  forwarding  acknowledgement  MICP  packet.  This  is  needed  in  order  to  make 
certain  that  the  previous  MSS  will  not  attempt  to  handle  traffic  for  an  MH  that  it 
no  longer  serves.  After  an  MH  moves  away  from  a  cell,  its  MSS  will  continue  receiv- 
ing IPIP-encapsulated  datagrams  for  it  from  other  MSSs.  It  knows  where  the  MH 
migrated,  so  it  informs  those  MSSs  of  the  fact  by  sending  them  a  redirection  MICP 
packet. 

When  an  MSS  receives  from  another  MSS  a  datagram  for  an  MH  that  it  no  longer 
serves,  but  for  which  it  has  in  its  cache  the  address  of  the  MSS  currently  serving  it,  it 
will  send  a  redirection  packet  to  the  sender.  To  avoid  having  to  rely  on  timeouts  from 
higher-level  protocols,  the  MSS  will  also  attempt  to  deHver  the  received  datagram 
to  where  it  thinks  it  should  have  gone.  If  the  MH  had  moved  again,  this  delivery 
will  result  in  a  new  redirect  being  sent.  Eventually,  all  MSSs  involved  in  a  particular 
MH's  traffic  will  either  converge  to  point  to  the  same  MSS,  or  their  information  will 
expire. 

Assuming  that  the  MH  does  not  change  cells  faster  than  these  packets  can  travel, 
the  network  should  be  able  to  track  the  MH  at  all  times.  The  only  danger  with  this 
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redirect  packet  being  lost  is  that  the  originating  MSS  will  not  be  informed  of  the 
MH's  movement.  This  does  not  cause  problems,  since  if  another  packet  is  received, 
another  redirect  will  be  sent.  Eventually,  one  will  reach  the  offending  MSS.  Hence, 
no  explicit  acknowledgement  of  the  redirect  packet  is  needed. 

When  an  MSS  is  asked  to  deliver  a  packet  to  an  MH  it  does  not  know  about,  it 
must  discover  which  MSS  is  responsible  for  it.  It  drops  the  packet  and  then  sends  an 
MICP  packet  to  all  the  MSSs  on  the  campus  to  find  who  is  responsible  for  it.  The 
MSS  currently  handling  the  MH  will  respond  with  an  MICP  packet.  This  on-demand 
discovery  of  peer  MSSs  allows  them  to  know  only  the  MSSs  involved  in  their  MHs' 
traflBc. 

2.4    Chapter  Summarv 

The  existing  mobile  data  networks  can  be  largely  categorized  into  two  different 
groups  depending  on  what  applications  are  intended:  TCP/IP  applications  and  pack- 
etized  voice  and  data  services.  In  the  mobile  data  networks  such  as  IPIP  and  VIP  for 
TCP/IP  applications,  a  common  approach  is  to  implement  the  IP-level  mobile  host 
protocols  and  use  the  internet  for  base  station  networking.  On  the  other  hand,  in 
the  mobile  data  networks  such  as  CPS  for  personal  communications,  the  importance 
of  base  station  networking  has  been  emphasized  and  MAN  technologies  are  used  for 
distributed  mobility  management  at  the  data  link  layer.  However,  it  has  not  consid- 
ered TCP/IP  environments,  which  accommodate  a  large  population  of  the  internet 
users  and  applications. 

The  next  chapter  is  about  the  virtual  cell  approach  to  solve  the  problems  from 
integrating  host  mobility  into  the  IP  layer  and  using  an  internet  for  base  station 
networking  in  TCP/IP  environments. 


CHAPTER  3 
THE  VIRTUAL  CELL  SYSTEM 


This  chapter  describes  a  virtual  cell  concept  for  the  transmission  of  IP  datagrams 
in  mobile  computer  communications.  A  virtual  cell  consists  of  a  group  of  physical 
cells  whose  base  stations  are  implemented  by  remote  bridges  and  interconnected  via 
high-speed  datagram  packet-switched  networks,  supporting  host  mobility  at  the  data 
link  layer.  Thus,  as  far  as  the  IP  layer  is  concerned,  it  appears  as  if  the  communication 
between  two  mobile  hosts  in  different  physical  cells  were  taking  place  directly  as  in 
the  same  physical  cell.  It  eliminates  the  necessity  of  IP-level  protocols  for  mobihty 
management,  which  may  interfere  with  the  conventional  IP  protocol  in  a  practical 
sense,  and  achieves  a  logically  flexible  coverage  area  according  to  mobility  and  data 
traffic  patterns  among  physical  cells.  Virtual  Cell  Protocol  (VCP),  which  consists 
of  handoff,  address  resolution,  and  data  transfer  modules,  is  designed  based  on  the 
distributed  hierarchical  location  information  of  mobile  hosts. 

3.1    Virtual  Cell  Concept 

Consider  the  transmission  of  IP  datagrams  between  two  MHs  within  a  physical 
cell.  Because  radio  links  provide  the  broadcast  ability,  the  source  can  deliver  IP 
datagrams  to  the  destination  using  the  normal  Address  Resolution  Protocol  (ARP). 
Thus,  the  broadcast  ability  of  radio  links  eliminates  the  necessity  of  locating  MHs, 
resulting  in  no  mobility  management  protocols  at  the  IP  layer  and  no  modifications  to 
the  normal  ARP  protocol.  To  apply  a  similar  rationale  to  MHs  crossing  the  physical 
cell  boundary,  we  propose  the  virtual  cell  concept,  as  depicted  in  Figure  3.1. 
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Figure  3.1.  Virtual  Cell  Concept 

A  virtual  cell  is  a  logically  extended  cell  of  physical  cells  whose  BSs  are  imple- 
mented by  remote  bridges  and  interconnected  via  a  base  station  network.  Because 
remote  bridge  BSs  make  it  possible  to  preserve  the  interconnection  level  of  physical 
cells  at  the  data  link  layer,  the  whole  virtual  cell  is  assigned  an  IP  network  address.  In 
order  to  make  it  possible  for  MHs  in  a  virtual  cell  to  directly  deliver  IP  datagrams  us- 
ing the  conventional  IP  and  ARP  protocols,  a  mobility  management  protocol,  called 
Virtual  Cell  Protocol  (VCP),  is  implemented  at  the  data  link  layer  of  BSs.  Each 
virtual  cell  is  also  allocated  an  ARP/Location  server  (ALS),  which  implements  VCP 
and  supports  the  normal  gateway  function  with  other  virtual  cells  or  fixed  networks. 
VCP  supports  handoff,  address  resolution,  and  location  tracking  for  datagram  deliv- 
ery, based  on  the  distributed  hierarchical  location  information  of  BSs  and  the  ALS, 
as  described  in  detail  in  Section  3.3. 
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Figure  3.2.  Data  Path  Between  Virtual  Cells 

It  should  be  noted  that  the  identification  and  location  information  of  an  MH 
is  represented  by  two  different  data  link  layer  addresses  in  the  virtual  cell,  Radio 
network  MAC  (RMAC)  address  and  Base  station  network  MAC  (BMAC)  address, 
respectively.  Note  that  the  BMAC  address  is  not  only  used  as  the  location  information 
of  the  MH  but  also  as  the  identification  information  of  a  virtual  cell.  For  example, 
the  SMDS  addresses  can  be  formatted  with  a  similar  structure  used  for  the  North 
American  Numbering  Plan  to  represent  the  geographic  semantics.  This  separation 
eliminates  the  necessity  of  acquiring  a  temporary  address  to  represent  the  MH's 
current  location  as  it  migrates.  In  addition,  unlike  IP-level  mobile  host  protocols, 
the  IP  network  address  also  represents  the  correct  location  information  of  the  MH 
while  it  moves  to  a  different  physical  cell  in  the  virtual  cell.  This  is  possible  because 
the  whole  virtual  cell  is  assigned  an  IP  network  address  and  host  mobility  is  shielded 
from  the  IP  layer. 

Even  when  the  MH  migrates  between  virtual  cells,  the  IP  network  address  is  still 
valid  as  the  near-location  information  in  the  virtual  cell  environment.  Figure  3.2 
shows  the  data  path  from  a  fixed  host  or  an  MH  in  a  different  network  to  MHi, 
which  migrates  from  its  native  virtual  cell  A  to  virtual  cell  B  to  virtual  cell  C. 
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Whenever  MHi  crosses  the  boundary  of  virtual  cells,  a  one-hop  forwarding  pointer 
is  maintained  from  ALSa  in  its  native  virtual  cell  to  the  ALS  of  its  current  virtual 
cell,  ALSb  or  ALSc,  by  the  handoff  procedure  of  VCP.  Note  that  the  forwarding 
pointer  is  not  maintained  at  the  IP  layer  but  at  the  VCP  layer.  When  a  fixed  host  or 
an  MH  in  a  different  network  wants  to  send  an  IP  datagram  to  MHi  in  virtual  cell 
B  or  C,  the  existing  routing  protocols  of  the  IP  layer  correctly  deliver  it  to  ALSa- 
Then  ALSa  redirects  it  at  the  VCP  layer  to  the  corresponding  ALS,  ALSb  or  ALSc, 
which  keeps  track  of  the  exact  location  of  MHi  in  its  virtual  cell. 

Thus,  every  MH  resides  in  its  native  virtual  cell  from  the  IP  layer's  viewpoint  but 
practically  may  reside  in  an  adjacent  virtual  cell  because  of  a  one- hop  forwarding 
pointer  from  the  native  virtual  cell  to  the  current  virtual  cell  maintained  at  the  VCP 
layer.  It  means  that  the  combination  of  the  existing  routing  protocols  of  the  IP  layer 
and  the  forwarding  pointer  between  virtual  cells  at  the  VCP  layer  gives  the  same 
effect:  that  we  have  a  fully  duplicated  location  information  among  virtual  cells  for 
the  global  mobile  network. 

To  configure  a  number  of  neighboring  physical  cells  into  virtual  cells,  the  underly- 
ing transport  network  should  provide  high-speed  datagram  packet-switched  services 
and  the  group  addressing  capability  at  the  data  link  layer.  The  flexible  configura- 
tion allows  network  managers  and  designers  to  customize  the  logical  coverage  area  of 
the  virtual  cell  to  their  requirements  according  to  host  mobiHty  and  communication 
patterns  among  physical  cells. 

Consider  an  environment  in  which  a  number  of  physical  cells  are  deployed  in  a 
metropolitan  area.  The  cell  size  in  the  urban  area  is  usually  smaller  than  that  in  the 
suburban  area  in  order  to  accommodate  a  larger  population  of  mobile  users.  Even 
though  user  mobiHty  is  inherently  unpredictable,  there  could  be  a  high  possibihty 
that  the  daily  routine  of  mobile  users  is  usually  confined  to  several  physical  cells  in 
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Figure  3.3.  Virtual  Cell  Architecture 

the  urban  area.  Because  of  a  relatively  small  cell  and  large  population  of  mobile  users 
in  the  urban  area,  it  becomes  very  important  to  achieve  the  fast  location  tracking  of 
mobile  users  and  the  small  cell  svi^itching  latency  between  physical  cells,  mitigating 
the  effect  of  a  large  volume  of  mobility  management  information.  Thus,  the  physical 
cells  covering  the  urban  area  could  be  combined  into  one  or  a  few^  number  of  virtual 
cell(s)  for  efficient  mobility  management  and  data  delivery. 

3.2    Virtual  Cell  Architecture 

Because  the  virtual  cell  is  a  logical  concept,  the  same  transport  network  may 
support  several  virtual  cells  simultaneously.  As  shown  in  Figure  3.3,  there  are  two 
roles  of  the  transport  network:  base  station  network  and  backbone  network.  The 
base  station  network  is  used  to  interconnect  a  number  of  BSs  and  an  ALS  to  build 
a  virtual  cell,  and  the  backbone  network  is  used  to  interconnect  among  virtual  cells 
and  fixed  networks.  Note  that  the  base  station  network  utilizes  both  point-to-point 
and  multicast  communications,  while  the  backbone  network  uses  only  point-to-point 
communication . 
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3.2.1    Base  Station  Network 

The  base  station  network  should  meet  some  requirements  from  two  different  view- 
points: 

•  Application's  viewpoint.  IP  operates  under  the  assumption  that  packet  arrivals 
are  independent  and  unpredictable.  If  IP  is  coupled  with  the  connection- 
oriented  base  station  network,  one  connection  establishment  and  release  may- 
be needed  for  each  individual  packet  over  a  fixed  link  between  a  pair  of  BSs. 
Although  a  bundle  of  packets  could  be  transmitted  over  one  connection  by 
some  multiplexing  techniques,  when  communicating  MHs  migrate  to  different 
physical  cells  with  the  connection,  the  problem  becomes  complex  and  even  in- 
tractable. In  addition,  because  location  tracking  has  inherently  connectionless 
properties,  the  base  station  network  should  support  the  datagram  delivery.  The 
current  ubiquitous  coverage  of  TCP/IP  networks  and  applications  also  requires 
that  the  base  station  network  be  able  to  cover  a  wide  area. 

•  Mobility  management's  viewpoint.  Combined  with  the  distributed  location  in- 
formation of  MHs  among  BSs,  host  mobility  is  a  main  cause  of  the  inconsistency. 
In  order  to  maintain  the  consistency  of  the  distributed  location  information 
efficiently,  the  base  station  network  should  support  the  multicast  ability.  Fur- 
thermore, because  the  network  control  information  for  mobility  management 
is  expected  to  be  greatly  increased  as  the  cell  size  becomes  smaller,  the  base 
station  network  should  have  high-bandwidth  not  so  as  to  affect  the  normal  data 
traffic  greatly. 

We  consider  SMDS  as  an  example  with  these  characteristics.  Both  individual 
and  group  addressing  capabihties  combined  with  a  set  of  addressing-related  service 
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Figure  3.4.  Base  Station  Architecture 

features  (e.g.,  source  address  validation,  source  and  destination  address  screening) 
enables  to  create  a  number  of  logical  networks  over  SMDS.  Every  BS  and  ALS  is 
attached  to  a  Subscriber  Network  Interface  (SNI)  and  individually  addressed.  At  the 
same  time,  a  group  address  identifies  all  BSs  in  a  virtual  cell,  ensuring  that  each  group 
address  identifies  uniquely  only  one  set  of  individual  addresses.  The  number  of  SNIs 
to  be  supported  by  a  switching  system  is  at  least  256  SNIs  and  the  future  number  is 
up  to  4096  SNIs  within  a  Local  Access  Transport  Area.  Considering  that  the  range  of 
a  physical  cell  is  usually  3-5  km  in  the  urban  area  or  10-15  km  in  the  suburban  area, 
the  capacity  of  a  very  few  number  of  switching  systems  may  be  enough  to  support 
base  station  networking  within  a  metropolitan  area. 

3.2.2    Base  Station 

The  communication  architecture  in  a  BS  is  shown  in  Figure  3.4.  The  internal 
port  connected  to  a  physical  cell  implements  an  RMAC  entity,  while  the  external 
port  connected  to  the  base  station  network  implements  a  BMAC  entity.  The  MAC 
relay  entity  translates  the  information  format  between  the  physical  cell  and  the  base 
station  network.  The  protocol  identifier  field  of  both  the  RMAC  and  BMAC  frame 
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headers  is  set  for  LLC  type  1  Unnumbered  Information  format.  For  example,  the 
Higher  Layer  Protocol  Identifier  field  of  the  SMDS  Interface  Protocol  L3_PDU  used 
as  BMAC  frames  is  set  to  1.  The  LLC  Service  Access  Point  of  the  RMAC  frame  is 
set  for  Subnetwork  Access  Protocol  (SNAP),  while  that  of  the  BMAC  frame  is  set 
for  VCP. 

The  VCP  header  has  three  fields:  VCP  type,  VCP  length,  and  RMAC  address. 
The  VCP  type  is  set  to  one  of  ADDRESS,  DATA  and  HANDOVER  to  specify  the 
address  resolution,  data  transfer  and  handoff  modules  of  VCP,  respectively.  Depend- 
ing on  the  VCP  type,  the  information  field  of  the  VCP  frame  has  a  different  format. 
The  DATA  type  indicates  that  the  information  field  has  an  IP  datagram,  and  the 
information  fields  for  the  HANDOVER  and  ADDRESS  types  are  given  in  the  next 
section.  The  VCP  length  indicates  the  length  of  the  VCP  header.  Note  that  the 
RMAC  address  is  only  used  with  the  DATA  type. 

Each  BS  maintains  a  partial  location  information  of  the  virtual  cell  in  its  Internal 
Membership  Table  (IMT)  and  External  Membership  Table  (EMT).  IMT  keeps  track 
of  all  MHs  currently  roaming  in  its  physical  cell.  The  IMT  entry  for  an  MH  has  two 
fields,  the  IP  and  RMAC  addresses,  each  of  which  is  used  as  an  identifier  but  has  a 
diflFerent  role.  The  IP  address  is  used  as  an  identifier  by  the  higher  network  layers 
as  usual,  but  the  network  part  of  the  IP  address  represents  the  MH's  native  virtual 
cell  to  which  it  initially  belongs.  This  IP  address  obtained  by  the  handoff  procedure 
is  only  needed  for  fast  address  resolution,  as  described  in  the  next  section.  On  the 
other  hand,  the  RMAC  address  is  the  other  identifier  used  by  VCP  for  the  support  of 
host  mobility  at  the  data  link  layer.  EMT  has  a  relatively  small  number  of  entries  for 
MHs  roaming  in  different  physical  cells  of  the  same  virtual  cell.  The  EMT  entry  for 
an  MH  has  three  fields,  the  IP,  RMAC  and  BMAC  addresses,  respectively,  of  which 
the  BMAC  address  represents  the  MH's  current  network  address. 
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If  the  source  MH  wants  to  send  an  IP  packet  to  the  destination  MH  whose  native 
virtual  cell  is  the  same,  the  source  performs  ARP  and  sends  the  IP  packet  using  the 
destination  RMAC  address.  On  the  other  hand,  if  the  source  wants  to  send  an  IP 
packet  to  the  destination  whose  native  virtual  cell  is  different,  the  source  does  not 
perform  ARP  and  sends  the  IP  packet  using  the  RMAC  address  of  its  current  BS, 
which  is  obtained  by  a  beacon  message.  When  the  BS  receives  an  SNAP/LLC/RMAC 
frame  from  the  internal  port,  it  checks  the  Protocol  Identification  (PID)  field  of 
the  SNAP  header.  If  the  PID  field  is  set  for  ARP,  the  ARP  packet  is  sent  to  the 
address  resolution  module  which  performs  Virtual  Cell  Address  Resolution  Protocol 
(VCARP).  VCARP  achieves  the  IP-to-RMAC  address  binding  for  the  destination 
and  also  distributes  the  location  information  of  the  source  and  destination.  If  the 
PID  field  is  set  for  IP,  the  IP  packet  is  sent  to  the  data  transfer  module  with  the 
destination  RMAC  address.  If  the  destination  RMAC  address  is  the  BS's  RMAC 
address,  the  data  transfer  module  sets  the  RMAC  address  field  of  the  VCP  header  to 
the  BS's  RMAC  address  and  relays  a  VCP/LLC/BMAC  frame  to  the  ALS  without 
searching  IMT  and  EMT.  Otherwise,  the  data  transfer  module  keeps  track  of  the 
destination  MH's  current  location  using  IMT  and  EMT.  If  the  corresponding  entry 
is  not  found,  the  data  transfer  module  transmits  to  the  ALS  a  VCP/LLC/BMAC 
frame  whose  VCP  header  contains  the  destination  RMAC  address. 

3.2.3    ARP/Location  Server 

Each  ALS  maintains  Global  Membership  Table  (GMT)  for  MHs  roaming  in  its 
virtual  cell,  as  shown  in  Figure  3.5.  The  entry  format  of  GMT  is  the  same  as 
that  of  EMT.  There  are  two  types  of  frames  from  the  backbone  network  for  data 
transfer.  One  is  the  VCP/LLC/BMAC  frame  redirected  from  other  virtual  cells 
that  maintain  one-hop  forwarding  pointers  at  the  VCP  layer,  and  the  other  is  the 
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Figure  3.5.  Communication  Architecture  in  an  ARP/Location  Server 

SNAP/LLC/BMAC  frame  from  fixed  networks  or  other  virtual  cells.  The  first  type 
of  frame  is  sent  to  the  data  transfer  module  which  identifies  the  destination  BS  by 
searching  GMT  with  the  RMAC  address  of  the  VCP  header,  and  transmits  it  to  the 
base  station  network.  The  second  type  of  frame  is,  however,  directly  sent  to  the 
routing  module  which  extracts  the  IP  packet.  If  the  network  part  of  the  destination 
IP  address  indicates  the  same  virtual  cell,  the  packet  is  sent  to  the  data  transfer  mod- 
ule of  VCP  where  the  IP-to-RMAC  address  binding  occurs  and  a  VCP/LLC/BMAC 
frame  is  transmitted  to  the  base  station  network  or  to  the  backbone  network,  de- 
pending on  whether  the  destination  MH  resides  in  this  native  virtual  cell  or  moved 
to  another  virtual  cell,  respectively.  If  the  IP  packet  is  in  transit,  it  is  transmitted  to 
the  next-hop  router. 

The  VCP/LLC/BMAC  frame  from  the  base  station  network  is  always  sent  to  the 
VCP  module.  If  the  frame  carries  a  VCARP  packet,  the  address  resolution  module 
performs  the  IP-to-RMAC  address  binding  with  GMT  and  sends  it  back  to  the  base 
station  network.  If  the  frame  carries  an  IP  packet  with  a  BS's  RMAC  address  of 
the  virtual  cell  in  the  VCP  header,  the  data  transfer  module  dose  not  search  GMT 
and  directly  pass  the  IP  packet  to  the  routing  module  where  the  next-hop  router  is 
determined.  Otherwise,  the  data  transfer  module  searches  GMT  with  the  destination 
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RMAC  address  and  then  transmits  a  VCP/LLC/BMAC  frame  back  to  the  base 
station  network  or  to  an  adjacent  virtual  cell  following  a  forwarding  pointer.  At  the 
same  time,  if  the  destination  MH  resides  in  this  virtual  cell,  the  ALS  sends  to  the 
source  BS  a  VCARP  reply  which  conveys  the  destination  BMAC  address  so  that  the 
subsequent  data  transfer  can  directly  go  to  the  destination  BS. 

3.3    Virtual  Cell  Protocol 

3.3.1    Distributed  Location  Information 

There  are  generally  three  different  ways  to  distribute  location  information:  cen- 
tralized, partitioned,  and  duplicated.  Depending  on  which  way  is  used  to  treat  lo- 
cation information,  there  is  a  tradeoff  between  location  registration  and  paging.  In 
a  centralized  system,  a  large  volume  of  location  updates  at  one  physical  site  may 
degrade  the  network  performance  significantly.  However,  the  consistency  of  location 
information  is  obtained  and  simple  paging  is  achieved.  A  typical  example  is  the 
first  generation  of  mobile  cellular  telephone  systems.  In  a  partitioned  system,  the 
different  partitions  of  location  information  are  held  at  different  sites.  An  example  is 
IPIP  where  every  BS  maintains  its  local  location  information  and  paging  is  basically 
required  when  two  remote  physical  cells  are  involved  in  communication  unless  a  for- 
warding pointer  is  found.  Thus,  this  approach  can  alleviate  the  problem  of  location 
registration  in  the  centralized  system,  but  may  need  frequent  paging.  In  a  duplicated 
system,  on  the  other  hand,  the  same  location  information  is  held  at  the  different 
sites.  VIP  is  an  example  where  the  location  information  of  an  MH  is  held  on  several 
intermediate  gateways.  Although  it  can  give  an  optimal  routing  path  in  the  normal 
case,  when  a  communicating  MH  migrates,  the  inconsistent  location  information  may 
be  spread  over  several  intermediate  gateways. 
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Figure  3.6.  Distributed  Location  Information  in  the  Virtual  Cell  System 

In  order  to  support  mobility  management  in  a  distributed  manner,  the  virtual  cell 
takes  advantage  of  the  distributed  hierarchical  location  information  which  involves  a 
combination  of  partitioning,  duplication,  and  centralization,  as  shown  in  Figure  3.6. 
GMT  is  partitioned  between  virtual  cells.  On  the  other  hand,  IMT  is  partitioned  with 
other  IMTs  and  duplicated  with  GMT  in  a  virtual  cell.  EMT  partially  duplicates  a 
part  of  GMT  for  MHs  in  different  physical  cells  in  the  same  virtual  cell  so  that  address 
resolution  and  location  tracking  for  data  transfer  are  first  tried  to  accomplish  among 
BSs  independent  of  the  ALS.  Thus,  GMT  is  only  referred  in  case  that  they  cannot 
be  solved  among  BSs. 

It  should  be  noted  that  there  is  another  conceptual  hierarchy  above  the  GMT 
level.  When  MEx  migrates  from  its  native  virtual  cell  A  to  virtual  cell  5,  forwarding 
pointers  are  maintained  from  ALS  a  to  ALSb  to  the  current  BS  at  the  VCP  layer. 
From  the  prospective  of  the  IP  layer,  however,  ME^  is  regarded  as  if  it  were  in 
its  native  virtual  cell.  Thus,  when  a  remote  fixed  host  or  an  MH  in  a  third  network 
wants  to  send  a  packet  to  MHx ,  the  existing  routing  protocols  of  the  IP  layer  correctly 
deliver  the  packet  to  ALS  a-  Then  ALS  a  traces  MH^  using  the  forwarding  pointers. 
Therefore,  the  IP  network  address  of  a  migrating  MH  represents  at  least  the  near- 
location  information  of  the  MH,  and  when  coupled  with  the  existing  routing  protocols 


38 


■\^BS4j  MH4 
■X  ^BS3^  MH3 
(  BS2  )  MH2 


— »■  X  Message  transfer  failure 


 *■  Long  message  latency 


MH:  Mobile  Host 


BS:  Base  Station 


■»  Movement  of  a  MH 


Short  message  latency 


Figure  3.7.  An  Example  of  the  Inconsistency  Problem 


of  the  IP  layer  it  can  give  the  same  effect  that  we  maintain  a  conceptually  centralized 
server  for  the  whole  mobile  network.  If  MHi  migrates  within  a  virtual  cell,  the  IP 
network  address  of  MHi  gives  the  exact  location  information. 

3.3.2    Impact  of  Mobility 

The  location  registration  due  to  the  handoff  procedure  is  a  main  source  of  incon- 
sistent location  information.  To  achieve  the  consistency  of  distributed  information  is 
not  trivial  at  all  and  many  researchers  have  been  working  on  this  problem  in  fixed 
network  environments  [17,  18,  19].  In  the  virtual  cell  environment,  however,  host 
mobility  is  integrated  with  fixed  networks,  which  makes  the  problem  even  more  com- 
plex. To  solve  this  problem,  the  handoff  procedure  should  utilize  the  multicasting 
ability  of  the  base  station  network. 

Assuming  that  the  distributed  location  information  is  initially  consistent,  con- 
sider that  MHl  moves  from  BS4  to  BSi  which  periodically  broadcasts  its  identity, 
as  shown  in  the  Figure  3.7.  When  MHi  receives  a  beacon  packet,  it  sends  a  greeting 
message  containing  its  identity  to  jB^i.  BSi  detects  the  entry  of  MHi  to  its  local  cell, 
deletes  its  corresponding  EMT  entry  if  exists,  and  inserts  a  new  IMT  entry.  Then, 
BSi  multicasts  a  location  registration  message  to  every  BS  so  that  it  maintains  the 
consistent  location  information  for  MH^.  However,  some  BSs  may  fail  to  receive  the 
message  due  to  communication  errors  and  buffer  overflows.  The  long  propagation 
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delay  also  has  the  similar  effect  as  the  message  loss  at  a  point  in  time  due  to  mobil- 
ity. To  explore  a  practical  solution,  we  should  begin  with  the  assumption  that  the 
base  station  network  supports  the  atomic  multicasting  by  using  an  existing  protocol 
such  eis  the  Trans  protocol  [19].  The  basic  idea  of  the  Trans  protocol  is  that  acknowl- 
edgements for  multicast  messages  are  piggybacked  on  messages  that  axe  themselves 
multicast,  using  a  combination  of  positive  and  negative  acknowledgement  strategies. 
This  piggyback  feature  is  suitable  for  the  handoff  procedure  because  in  the  steady 
state,  if  some  MHs  move  out  then  there  will  be  a  high  possibility  that  other  MHs 
move  in. 

Note  that  although  the  atomic  multicast  is  supported,  the  message  loss  can  still 
occur  if  BSi  multicasts  the  location  registration  message.  Assume  that  BS3  failed  to 
receive  it  and  BS4  received  it  with  a  long  message  latency.  Then,  the  following  four 
cases  can  happen  when  other  MHs  transmit  a  message  to  MHi  during  the  process  of 
the  atomic  multicast: 

1.  If  BS4  has  the  old  location  information  for  MHi,  the  message  transmitted  from 
MHi  to  MHi  will  be  lost. 

2.  If  BS4  has  the  new  location  information  for  MHi,  the  message  transmitted 
from  MH4  to  MHi  will  be  received. 

3.  Because  BS3  has  the  old  location  information  for  MHi,  the  message  transmit- 
ted from  MH3  to  MHi  will  be  lost  or  received,  depending  on  what  value  BS4 
has.  If  it  has  the  old  value,  the  message  will  be  lost.  Otherwise,  the  message 
will  be  redirected  to  BSi  by  BS4. 

4.  Because  BS2  has  the  new  location  information  for  MHi ,  the  message  transmit- 
ted from  MH2  to  MHi  will  be  received. 
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From  the  above  observation,  we  can  see  that  at  least  the  previous  base  station 
BS4  must  have  the  new  location  information  for  MHi  as  soon  as  possible  in  order 
to  avoid  a  possible  message  loss  or  to  mitigate  the  effect  of  long  message  latency. 

3.3.3    HandofF  Procedure 

The  new  BS  informs  the  previous  BS  of  an  MH's  migration  using  point-to-point 
communication  immediately  after  detecting  the  MH's  identity  so  that  all  messages 
received  at  the  previous  BS  during  the  handofF  procedure  can  be  correctly  redirected 
to  the  new  BS.  Next,  the  previous  BS  is  responsible  for  maintaining  the  consistency 
of  the  MH's  location  information  at  all  BSs  using  an  atomic  multicasting.  Note 
that  the  ALS  is  excluded  from  the  multicasting  group.  If  the  ALS  is  involved,  every 
movement  of  the  MH  will  need  an  access  to  the  ALS  and  there  is  no  difference  from  the 
centraHzed  system.  However,  it  may  cause  another  inconsistency  problem  between 
the  ALS  and  BSs  because  the  ALS  may  have  old  location  information  for  some  MHs. 
Thus,  the  previous  BS  periodically  backups  the  location  information  updated  during 
a  predefined  time  interval  to  the  ALS  using  point-to-point  communication.  When  an 
MH  continues  to  move  to  different  physical  cells,  there  can  be  a  redirection  chain  with 
multiple  hops  from  the  first  BS  to  the  current  BS.  However,  the  chain  is  eliminated 
when  the  newest  location  information  is  received  through  an  atomic  multicasting. 

Figure  3.8  shows  a  communication  model  for  the  handoff  procedure.  For  conve- 
nient description,  we  define  the  following  notations  where  the  VCP  header  informa- 
tion is  omitted: 

•  ^  B,...,N  :  {frame}  means  that  A  multicasts  {frame}  to  B,...,N  over 
the  base  station  network,  where  {frame}  consists  of  {source  address,  destination 
address  |  frame  data}. 
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Figure  3.8.  A  Communication  Model  for  HandofF  Procedure 

A  — »  B  :  {frame}  means  that  A  sends  {frame}  to  B  via  a  point-to-point  link. 

A-^  B  :  {frame}  means  that  A  broadcasts  {frame}  whose  destination  is  B 
via  radio  links. 


•  IPjA)  denotes  the  IP  address  of  an  MH  A. 

•  RMAC{A)  denotes  the  radio  MAC  address  of  A,  where  A  may  be  a  BS  or  an 
MH. 

•  BMACjA)  denotes  the  address  of  A  in  the  base  station  network  (BMAC), 
where  A  may  be  a  BS  or  an  ALS. 

Within  a  virtual  cell.  Consider  that  MH^  moves  from  552^  to  BSnx-  In  addition 
to  its  own  addresses,  IP{MHa)  and  RMAC{MH^),  MHa  initially  keeps  the  current 
base  station  addresses,  RMAC{BS2^)  and  BMAC{BS2x). 


1.  BSnx MHa  •■  {RMAC{BSnx),  broadcast  address  |  BMAC{BSnx)} 
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MHa  receives  a  beacon  including  BMAC{BSnx)  from  BSnx,  decides  to  switch 
from  BS2X  to  BSnx,  and  updates  its  base  station  address  from  RMAC^BS^x) 
and  BMAC{BS2x)  to  RMAC{BSr,.)  and  BMAC{BSnx),  respectively. 

2.  MH,     BSnx  :  {RMAC{MHa),  RMAC{BSnx)  \  IP{MHa),  BMAC{BS2x)} 

BSnx  receiving  a  greeting  message  from  MHa  processes  its  IMT  and  EMT  using 
the  source  address  RMAC{MHa)  as  a  key.  If  the  EMT  entry  for  MHa  is  found, 
the  entry  is  deleted  and  then  the  IMT  entry  for  MHa  is  added.  IP{MHa)  in 
the  EMT  entry  is  used  for  fast  address  resolution,  as  described  in  Section  3.3.4. 

3.  BSnx       BS2X  :  {BMAC{BSnx),  BMACiBS^,)  \ 
IP{MHa),  RMAC{MHa),  BMAC{BSnx)} 

BS2X  receiving  the  notification  of  MHa^s  migration  sends  an  acknowledgement 
to  BSnx  and  then  processes  IMT  and  EMT.  The  IMT  entry  for  MHa  is  deleted 
and  the  new  EMT  entry  for  MHa  is  added. 

4.  BS2x  =J>  BS^x, . . . ,  BSnx  ■  {BMAC{BS2,),  multicast  address  | 
IPiMHa),  RMAC{MHa),  BMAC{BSn,)} 

Every  BS  receiving  the  notification  of  MHa's  migration  except  BSnx  processes 
its  EMT.  If  the  EMT  entry  for  MHa  already  exist,  the  entry  is  updated. 

Between  virtual  cells.  Consider  that  MHa  moves  from  BSnx  to  ^^i,^.  In  addition 
to  its  own  addresses,  IP{MHa)  and  RMAC{MHa),  MHa  keeps  the  current  base 
station  addresses,  RM  AC  (BSnx)  and  BMAC{BSna:)- 

1.  BSiy-^  MHa  •■  {RMAC{BSiy),  broadcast  address  |  BMAC{BSiy)} 

2.  MHa     BS^y  :  {RMACiMHa),  RMAC{BS,y)  \  IP{MHa),  BMAC{BSnx)} 
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BS\y  receiving  a  greeting  message  from  MHa  detects  from  the  previous  base 
station  address  BMAC{BSnx)  that  MHa  was  in  another  virtual  cell.  BSiy 
creates  the  IMT  entry  for  MHa  ■ 

3.  BSxy  — >  ALSy  :  {BMAC{BSiy),  BMAC{ALSy)  \ 
IP{MHa),  RMAC{MHa),  BMAC{BS^y),BMAC{BSn.)) 

ALSy  receiving  a  location  registration  message  creates  the  GMT  entry  for  MHa 
using  IP{MHa),  RMAC{MHa),  and  BMAC{BSiy).  Then,  ALSy  knows  that 
BSnx  is  a  base  station  in  virtual  cell  X,  using  BMAC{BSnx)- 

4.  ALSy  — ^  ALSx  :  {BMAC{ALSy),  BMAC{ALSx)  \ 
IP{MHa),  RMAC{MHa),  BMAC{ALSY),BMAC{BSn.)} 

ALSx  receiving  a  location  registration  message  updates  the  GMT  entry  for 
MHa  using  IP{MHa),  RMAC{MHa),  and  BMAC{ALSy).  This  entry  serves 
as  a  forwarding  pointer  at  the  VCP  layer  if  ALSx  is  the  ALS  of  MHa's  native 
virtual  cell.  If  it  is  not,  ALSx  informs  the  ALS  of  MH^s  native  virtual  cell  of 
MHa's  migration  to  virtual  cell  Y  so  that  a  forwarding  pointer  is  maintained 
from  MHa^s  native  virtual  cell  to  virtual  cell  Y. 

5.  ALSx  — *  BSnx  :  {BMAC{ALSx),  BMAC{BSnx)  \ 
IPiMHa),  RMAC{MHa),  BMAC{ALSx)} 

BSnx  receiving  the  notification  of  MHaS  migration  sends  an  acknowledgement 
to  ^^ij,  in  the  reverse  direction  and  then  processes  IMT  and  EMT.  The  existing 
IMT  entry  is  deleted  and  an  EMT  entry  using  IP  {MHa),  RMAC{MHa),  and 
BMAC{ALSx)  is  added. 

6.  BSnx  =^  BSix,BS2x, . .  .,BS(n-i)x  ■  {BMAC{BSnx),  multicast  address  | 
IP{MHa),  RMAC{MHa),  BMAC{ALSx)} 
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Every  BS  receiving  the  notification  of  MHa's  migration  except  BSnx  manipu- 
lates its  EMT.  If  the  EMT  entry  exists,  it  is  updated. 

3.3.4    Address  Resolution 

VCARP  is  designed  to  get  the  physical  address  of  the  MH  in  a  remote  physical 
cell  of  the  virtual  cell  and  should  satisfy  at  least  the  following  three  requirements. 
First,  in  order  to  achieve  the  compatibility  with  ARP,  VCARP  must  be  shielded  from 
MHs  as  if  they  were  in  a  physical  cell.  Second,  because  the  VCARP  packet  latency 
can  directly  affect  the  scalability  of  a  virtual  cell,  a  distributed  address  resolution 
mechanism  should  be  considered  rather  than  broadcasting  or  a  centralized  solution. 
Third,  when  an  MH  moves  to  an  adjacent  physical  cell  immediately  after  sending  an 
ARP  request  and  further  moves  to  another  physical  cell,  the  ARP  reply  may  be  lost 
and  then  the  ARP  request  may  be  repeatedly  generated.  Hence,  VCARP  must  also 
support  host  mobility. 

Figure  3.9  shows  the  ARP/VCARP  packet  format.  The  ARP  packet  follows  ex- 
actly the  same  format  as  the  existing  ARP  standard  and  has  longer  fields  SENDER 
HA  and  TARGET  HA  in  order  to  make  it  possible  to  accommodate  the  use  of  hi- 
erarchical 64-bit  E.164  network  addresses  in  radio  networks.  The  VCARP  packet 
format  has  an  additional  64-bit  field  BASE  HA  which  is  used  to  convey  the  location 
information  which  is  expected  to  be  used  for  data  transfer  in  the  very  near  future. 
The  source  MH  performs  ARPing  only  when  it  resides  in  its  native  virtual  cell  and 
the  destination  MH  belongs  to  the  same  native  virtual  cell.  Otherwise,  it  directly 
transmits  IP  packets  without  ARPing. 

Within  a  physical  cell.  The  procedure  is  the  same  as  in  the  normal  ARP.  The 
same  ARP  request/reply  packet  formats  are  used  where  an  8-bit  hardware  address 
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Figure  3.9.  ARP/VCARP  Packet  Format 

length  field  allows  to  accommodate  arbitrary  radio  network  addresses.  The  IP  pro- 
tocol of  the  source  MH  checks  the  destination  IP  address  with  its  address  resolution 
table.  If  the  address  binding  is  successful,  the  source  MH  uses  the  RMAC  address 
to  transfers  datagrams.  Otherwise,  the  source  MH  broadcasts  an  ARP  request.  As 
soon  as  the  destination  MH  receives  the  ARP  request,  it  responds  with  an  ARP  reply 
that  includes  the  requested  physical  address  of  the  destination.  At  the  same  time, 
the  BS  that  received  the  ARP  request  searches  IMT  and  knows  that  the  destination 
MH  is  in  the  local  physical  cell  so  it  discards  the  ARP  request. 

Within  a  virtual  cell.  Figure  3.10(a)  shows  the  flow  of  messages  when  MHi 
broadcasts  an  ARP  request  to  get  the  physical  address  of  MH2  and  it  is  resolved 
at  BSi.  BSi  that  received  the  ARP  request  searches  its  EMT.  The  EMT  entry  for 
MH2  already  contains  its  physical  address  and  BSi  directly  sends  an  ARP  reply 
within  the  physical  cell.  At  the  same  time,  BS\  sends  a  VCARP  reply  to  BS2. 
Then,  BS2  transforms  the  format  of  the  VCARP  reply  to  that  of  the  ARP  reply  by 
stripping  out  the  field  BASE  HA  and  broadcasts  the  transformed  ARP  reply  in  its 
local  cell.  This  VCARP  reply  also  conveys  the  location  information  of  MHi  to  BS2 
through  the  field  BASE  HA  which  contains  the  network  address  of  ^^i.  Then,  BS2 
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creates  an  EMT  entry  so  that  MH2  can  directly  transmit  data  to  MH^  without  any 
additional  address  resolution  in  the  near  future  communication.  It  is  based  on  the 
observation  that  computer  communication  is  usually  bidirectional;  if  MHi  has  some 
reason  to  talk  to  MH2,  then  MH2  will  probably  have  some  reason  to  talk  to  MHy. 

Figure  3.10(b)  shows  the  flow  of  messages  when  the  address  resolution  is  accom- 
plished at  the  ALS.  BSi  sends  a  VCARP  request  to  the  ALS  because  it  has  no  entry 
for  MH2.  Then,  the  ALS  sends  a  VCARP  reply  to  BS^  after  resolving  M/fz's  ad- 
dress binding  using  GMT.  At  the  same  time,  the  ALS  also  sends  BS2  a  VCARP 
reply  with  the  same  reason  as  in  the  above  case.  The  location  information  of  MH\ 
and  MH2  is  also  conveyed  via  both  VCARP  replys  to  establish  a  logical  bidirectional 
link  between  them.  Now  we  describe  the  support  of  host  mobihty  in  VCARP.  Con- 
sider the  case  that  MHx  requested  an  address  resolution  to  send  datagrams  to  MH2 
and  then  moved  to  a  third  BS^  before  the  VCARP  reply  is  arrived.  By  the  handoff 
procedure,  BS^  knows  that  MHi  currently  belongs  to  BSn-  As  soon  as  BSi  receives 
the  VCARP  reply,  it  redirects  the  packet  to  BSn- 

3.3.5    Data  Transfer 

The  transmission  of  packets  between  virtual  cells  and  fixed  networks  have  been 
described  in  the  previous  sections.  This  section  deals  with  the  packet  transmission 
within  a  physical  cell  and  within  a  virtual  cell.  Let  BSij  denote  the  i-th  BS  of  virtual 
cell  j  in  Figure  3.6,  where  i  =  1, . . . ,  n  when  j  =  A  and  i  =  1, . . . ,  m  when  j  =  B. 
We  assume  that  MHi,  belongs  to  BSij  and  initially  belongs  to  BS-^a-  We  also 
assume  that  the  address  binding  is  achieved  by  VCARP. 

Within  a  physical  cell.  Consider  the  case  that  MH^a  wants  to  send  a  packet  to 
MH^.  By  the  address  resolution  protocol,  MH^a  knows  the  RMAC  address  of  MH^ 
and  directly  broadcasts  the  packet  to  it.  BSia  knows  that  MH^  resides  in  its  physical 
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cell  by  checking  IMT  using  the  RMAC  address  of  MH^:  as  a  key,  and  discards  the 
received  packet. 

Within  a  virtual  cell.  Consider  the  case  that  MH^  moves  from  ^^i^  to  55*2^ 
and  each  of  MHia,  MH2A,  and  MHnA  wants  to  send  a  packet  to  MH^^.  By  the 
handolf  procedure,  a  forwarding  pointer  is  set  from  -S^i^  to  BS2A  at  the  EMT  entry 
of  BSiA  and  the  packet  from  MHxa  is  correctly  delivered  to  MHx  by  the  forwarding 
pointer.  The  packet  from  MH2A  is  directly  delivered  to  MH^  and  BS2A  discards  it 
because  the  IMT  entry  for  MH:^  is  found.  The  delivery  of  the  packet  from  MH^a 
has  three  cases.  If  the  EMT  entry  for  MH^  is  found  before  the  completion  of  ^^i^'s 
multicasting  operation,  the  packet  will  be  delivered  to  BS2A  via  BSia-  If  the  EMT 
entry  is  found  after  the  completion  of  its  multicasting  operation,  the  packet  will  be 
directly  dehvered  to  BS2a-  If  the  EMT  entry  is  not  found,  the  first  packet  from  MHua 
will  be  delivered  to  BS2A  through  the  ALS.  At  the  same  time,  the  ALS  generates  a 
VCARP  reply  in  order  to  have  BS'nA  learn  that  MH^  is  currently  belongs  to  BS2A- 
It  makes  it  possible  to  directly  transfer  the  subsequent  packets  following  the  first  one 
to  BS2A  without  going  through  the  ALS  repeatedly. 

3.4    Comparison  and  Discussion 

Figure  3.11  summarizes  the  characteristics  of  the  three  representative  mobile  data 
networks  and  compares  them  with  those  of  the  virtual  cell  system.  When  considering 
the  ubiquitous  coverage  of  TCP/IP  networks  and  applications,  the  virtual  cell  has 
several  advantages  against  IP-level  mobile  host  protocols.  Because  the  virtual  cell 
uses  a  underlying  datagram  network  to  interconnect  BSs  and  to  form  flexible  coverage, 
the  difficulties  arising  from  diverse  underlying  networks  can  be  avoided  in  terms  of 
the  performance  of  mobility  management  function  and  the  complexity  of  network 
management  function.  Moreover,  because  the  virtual  cell  is  a  logical  cell  and  the 
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Figure  3.11.  Comparison  of  the  Virtual  Cell  System  With  Other  Networks 


bridge  function  is  usually  much  faster  than  the  gateway  function,  it  is  possible  to 
increase  performance  significantly  if  we  can  properly  engineer  the  deployment  of 
virtual  cells  so  that  the  intra  virtual  cell  traffic  is  as  much  larger  than  the  inter 
virtual  cell  traffic  as  possible.  Preserving  the  interconnection  level  of  BSs  at  the  data 
link  layer  also  makes  it  possible  to  shield  host  mobility  from  the  IP  layer,  and  using 
the  base  station  network  address  as  a  contact  point  of  an  MH  makes  it  possible  to 
avoid  acquiring  a  temporary  address  by  some  kind  of  local  utilities  whenever  the  MH 
migrates.  Thus,  the  problems  of  interfering  with  the  conventional  IP  layer  and  of 
reserving  a  large  amount  of  IP  address  space  for  location  information  can  be  solved. 

There  are  two  major  issues  related  to  the  performance  of  VCP.  Within  a  virtual 
cell,  a  BS  first  tries  to  perform  address  resolution  and  data  transfer  locally  without 
the  intervention  of  the  ALS,  based  on  its  two  membership  tables,  IMT  and  EMT. 
Thus,  the  hit  ratio  of  the  EMT  search  directly  affect  the  degree  of  referring  to  the 
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ALS  when  remote  MHs  in  different  physical  cells  are  involved  in  communication. 
This  issue  directly  addresses  the  size  and  update  scheme  of  EMT,  and  it  is  closely 
related  to  the  IP  traffic  and  host  mobility  pattern.  A  number  of  research  results  on 
the  internet  protocol  traffic  analysis  in  fixed  network  environment  reveals  that  certain 
hosts  communicate  more  with  one  another  than  with  other  hosts  [39,  21,  22].  Even  in 
the  virtual  cell  environment,  the  locality  of  destination  hosts  can  have  an  important 
implication  for  each  of  the  address  resolution  and  data  transfer  functions.  A  small 
ARP  cache  at  each  MH  can  greatly  reduce  VCARP  traffic  within  a  virtual  cell.  The 
address  resolution  should  not  severely  burden  the  size  of  EMT.  If  we  do  not  consider 
host  mobility,  the  locality  property  also  means  that  data  transfer  does  not  severely 
burden  the  size  of  EMT.  When  considering  host  mobility,  the  update  scheme  of  EMT 
could  be  rather  important  if  we  can  properly  deploy  virtual  cells  based  on  mobility 
and  traffic  patterns  among  physical  cells. 

The  other  issue  is  the  amount  of  multicast  traffic  among  BSs,  which  is  generated 
only  by  the  handoff  function  of  the  virtual  cell.  Because  the  multicasting  ability  is 
supported  by  switches  in  the  base  station  network,  each  BS  is  not  related  to  the 
load  of  generating  the  multicast  traffic,  and  the  bandwidth  usage  of  the  base  station 
network  is  affected  by  it.  The  relatively  short  length  of  the  handoff  message  should 
not  severely  affect  the  high-bandwidth  base  station  network.  This  will  be  addressed 
in  Chapter  6. 

3.5    Chapter  Summary 

We  have  presented  a  new  network  infrastructure  for  the  transmission  of  IP  data- 
grams in  mobile  computer  communications.  With  the  virtual  cell  concept,  physical 
cells  are  grouped  into  larger  logical  cells  where  host  mobility  is  shielded  from  the  IP 
layer.  It  eliminates  the  need  of  IP-level  protocols  for  host  mobility  within  a  virtual 
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cell  and  provides  the  flexible  coverage  area  that  can  be  properly  engineered  according 
to  mobility  and  communication  patterns  among  physical  cells.  In  order  to  achieve 
this  concept,  we  have  described  the  virtual  cell  protocol  that  consists  of  handofF, 
address  resolution,  and  data  transfer  modules,  based  on  the  distributed  hierarchical 
location  information  of  mobile  hosts. 

Given  mobility  and  communication  patterns  among  physical  cells,  the  next  chap- 
ter deals  with  the  problem  of  deploying  virtual  cells  in  highway  cellular  systems. 
The  problem  is  transformed  to  an  optimization  problem  of  finding  a  cover  of  disjoint 
clusters  of  physical  cells,  where  each  cluster  corresponds  to  a  virtual  cell. 


CHAPTER  4 

OPTIMAL  PARTITIONING  IN  HIGHWAY  CELLULAR  SYSTEMS 

Given  a  linear  array  of  n  base  stations  which  generate  multiple  types  of  traffic 
among  themselves,  we  consider  the  problem  of  finding  a  cover  of  disjoint  clusters  of 
neighboring  base  stations.  The  objective  is  to  minimize  the  total  communication  cost 
for  the  entire  system  where  the  cost  of  intra-cluster  communication  is  lower  than  that 
of  inter-cluster  communication  for  each  type  of  traffic.  The  optimization  problem  is 
transformed  into  the  dual  based  on  the  relative  cost  that  is  the  communication  cost 
if  a  pair  of  nodes  are  in  different  clusters  of  a  partition  and  reflects  multiple  types 
of  traffic  with  multiple  dynamic  cost  functions.  Using  the  relative  cost  matrix,  an 
efficient  algorithm  of  O(mn^),  where  m  is  the  number  of  clusters  in  a  partition,  is 
designed  by  dynamic  programming. 

4.1    Problem  Formalization 

4.1.1    Problem  Statement 

The  base  stations  in  a  highway  cellular  system  axe  ordered  linearly  as  stations 
1,2, ...  ,n.  We  are  interested  in  finding  clusters  of  neighboring  stations,  such  that  a 
location  server  can  be  allocated  to  each  cluster.  In  other  words,  the  partition  of  base 
stations  must  follow  the  underlying  topology  constraint,  i.e.,  a  linear  arrangement  in 
highway  cellular  systems. 

The  communication  among  base  stations  is  considered  as  a  full  mesh  of  point-to- 
point  logical  network.  This  network  represents  a  possible  communication  between 
mobile  hosts  through  different  base  stations.  The  network  is  described  by  a  complete 
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directed  graph  G  =  {V,E),  where  |y|  =  n.  The  vertices  of  the  graph  represent 
the  base  stations  of  the  network,  and  their  indices  are  numbered  from  1  to  n  in 
sequence  for  n  concatenated  cells.  The  edges  represent  directional  communication 
links  between  vertices,  and  each  of  which  is  assigned  a  move  frequency  by  a  function 
fjn  -.V  xV  ^  and  a  find  frequency  by  //  :  V  x  V  — >  Denote  w]^  and  w"^  as 
the  weight  of  a  move  operation  within  a  cluster  and  between  clusters,  respectively, 
and  w^j  and  lu^  as  that  of  a  find  operation  within  a  cluster  and  between  clusters, 
respectively.  Then  we  define  a  cost  function  :  x  V  — >  ,  where  i  =  m  or  f  and 
J  =  1  or  2,  such  that  for  all  (u,  u)  G  F  x  V, 

4{u,v)    =    fi{u,v)wi.  (4.1) 

The  set  of  equations  given  by  (4.1)  implies  that  each  edge  is  dynamically  assigned  one 
of  two  different  costs  for  each  type  of  operation,  depending  on  whether  it  is  involved 
in  a  cluster  or  between  clusters.  Given  the  number  of  clusters  m  for  1  <  m  <  n,  a 
sequence  {1, . . . ,  n}  of  n  vertices  can  be  linearly  partitioned  into  11  =  {Pi, . . . ,  Pm} 
such  that  PiDPj  =  0  and  U,P,-  =  V,  where  i,j  —  1, . . . ,  /7i  and  i  ^  j.  The  intra-cluster 
communication  cost  of  IT  is  then 

C0Sti^,ra{U)  =  E    E  +  ^))  (4-2) 

1=1  u,veP, 

and  the  inter-cluster  communication  cost  of  11  is 

m-l 

CosUnUr{n)=Yl        i:       {icl{u,v)  +  cj{u,v))  +  {ci{v,u)  +  c}{v,u))).  (4.3) 

•=1  uGP,,v€Pj,:<i 

From  (4.2)  and  (4.3),  the  total  communication  cost  of  11  is  therefore 

CosttotalC^)  =  Costintrai^)  +  Cost„UeriU).  (4.4) 

Figure  4.1  illustrates  a  network  model  of  size  6  with  two  frequency  matrices, 
and  //.  When  u;^  =  u;}  =  1,  u;^  =  3,  and  wj  -  2,  the  cost  matrices  of  Figure  4.2  are 


54 


VI   V2  V3   V4  V5  V* 


2    3    4    5  6 


2    3    4    5  6 


0 

2 

4 

6 

3 

2 

4 

0 

7 

2 

2 

1 

3 

5 

0 

7 

5 

3 

2 

6 

3 

0 

1 

5 

6 

2 

2 

4 

0 

3 

4 

3 

5 

3 

6 

0 

The  phyrictl  topology    The  logical  catnmunicaticgi  network   Tlie  move  frequency  matrix        The  find  frequency  matrix 

Figure  4.1.  An  Example  of  the  Network  Model 
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(a)  The  intra-cluster  cost  matrices 


(b)  The  inter-cluster  cost  matrices 


Figure  4.2.  Intra-  and  Inter-cluster  Cost  Matrices 

driven  from  (4.1).  If  a  partition  U  =  {Pi,P2,^3},  where  Pi  =  {1,2},  =  {3,4}  and 
P3  =  {5,6},  Co5f,„<ro(n)  and  Costinteri^)  IS  the  sum  of  all  elements  in  the  shaded 
areas  of  Figure  4.2  (a)  and  (b),  respectively.  Thus,  Costtotaii^)  =  Costintrai^)  + 
C0Stinter{U)  =U  +  493  =  537. 

Given  n  and  m,  our  task  is  to  find  an  optimal  partition  11  which  minimizes 
Costtotai  subject  to  1  <  |Pi|  <  6  <  n  for  every  i  =  1, . . . ,  m,  where  |P,|  is  the  number 
of  vertices  in  the  ith  cluster.  Since  the  cluster  size  constraint  does  not  add  the  time 
and  space  complexity  and  affect  our  algorithm,  we  first  present  our  solution  for  the 
general  case  without  the  constraint,  and  then  we  show  that  the  solution  is  simply 
modified  to  consider  the  constraint.  Before  presenting  our  solution,  let  us  point  out 
that  exhaustively  checking  all  possible  partitions  does  not  yield  an  efficient  algorithm. 
For  example,  if  an  ordered  sequence  {1, 2, 3, 4}  of  4  vertices  is  given,  we  can  partition 
it  in  8  distinct  ways: 


{{1,  2,  3,  4}} 


if  m  =  1; 
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Table  4.1.  Number  of  Partitions 


The  number  of  clusters 

n 

m=l 

m=:2 

ni=3 

m=4  m 

=5 

m=6 

m=7 

m=8 

m=9 

1 

1 

2 

1 

1 

3 

1 

2 

1 

4 

1 

3 

3 

1 

5 

1 

4 

6 

4 

1 

6 

1 

5 

10 

10 

5 

1 

7 

1 

6 

15 

20 

15 

6 

1 

8 

1 

7 

21 

35 

35 

21 

7 

1 

9 

1 

8 

28 

56 

70 

56 

28 

8 

1 

{{1,2,  3},  {4}},  {{1,2},  {3,  4}},  and  {{1},  {2,  3,  4}}  if  m  =  2; 

{{1,  2},  {3},  {4}},  {{1},  {2,  3},  {4}},  and  {{1},  {2},  {3,  4}}  if  m  =  3; 
and  {{1},  {2},  {3},  {4}}  if  m  =  4. 

Table  4.1  gives  the  number  of  ways  to  partition  a  sequence  {1,2,...,  n}  of  n  vertices 
into  m  clusters  for  1  <  m  <  n,  preserving  order.  Since  the  numbers  in  Table  4.1  form 
Pascal's  triangle,  if  given  n  and  m,  the  number  of  partitions  is  G{n,  m)  =  For 
example,  if  n  =  9  and  m  =  5,  G{9, 5)  —  (4)  -  70. 

Because  row  n  of  Table  1  sums  to  2"~^  by  the  binomial  theorem,  the  total  number 
of  partitions  of  n  ordered  vertices  is  2""^  for  n  >  1.  Thus,  the  time  complexity  of 
exhaustive  search  is  0(2"n^)  and  its  space  complexity  is  0(n'^)  for  every  m.  Given 
n  and  m,  it  has  0{G{n,m)n^)  time  complexity  and  O(n^)  space  complexity.  The 
brute-force  method  of  exhaustive  search  is  therefore  a  poor  strategy  for  determining 
an  optimal  partition.  The  results  of  our  algorithm  by  dynamic  programming  based 
on  a  relative  cost  matrix  show  0{n^)  time  complexity  and  0{n^)  space  complexity 
for  every  m,  and  O(mn^)  time  complexity  and  0{n)  space  complexity  for  each  m. 


56 


4.1.2    Dual  Problem 

The  dynamic  property  of  multiple  cost  functions  for  multiple  types  of  traffic  not 
only  makes  the  optimization  problem  complex  in  terms  of  computation  and  storage 
but  also  makes  it  difficult  to  design  an  efficient  algorithm.  It  is  therefore  desirable 
to  express  the  optimization  problem  with  multiple  cost  functions  as  an  equivalent 
problem  with  one  cost  function  without  affecting  the  objective.  Thus,  we  combine 
the  set  of  equations  given  by  (4.1)  into  a  relative  cost  function  g  :  V  xV  such 
that  for  all  (u,  v)  eV  xV, 

g{u,  v)  =  fm{u,  v){wl  -  wl)  +  fj{u,  v){wj  -  w]),  (4.5) 

where  {wf^  —  w}^)  is  the  relative  weight  of  a  move  operation  and  (wj  —  '^})  is  the 
relative  weight  of  a  find  operation.  g{u,v)  then  represents  the  total  relative  cost  of 
fm{u,v)  move  and  ff{u,v)  find  operations  from  vertex  u  to  vertex  v  if  they  are  in 
different  clusters.  Using  the  equation  (4.5),  the  dual  to  the  optimization  problem  can 
now  be  written  as 

m-l 

t=i  uePi,veP3,i<j 
s.t.      l<\Pi\<b<n,     i  =  l,...,m, 

where  |P,|  is  the  number  of  vertices  in  the  ith  cluster.  Let  us  now  define 

m  — 1 

RcosUr,t„{U)   =    J2       E      {g{u.v)  +  g{v,u)),  (4.6) 

!=i  uePi,veP],i<j 

m— 1 

t/pper,„,,,(n)    =    Yl        E       icl{u,v)  +  c]{u,v)),  (4.7) 

T=i  uePi,vePj,i<j 

m  — 1 

Lowenr^train)     =      Yl  E  (4 (u,  u)  +  c} (u,  u)) .  (4.8) 

1=1  uePi,vePj,i<j 

The  correctness  of  the  transformation  is  then  proved  by  the  lemma  and  theorem  be- 
low. 
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Lemma  4.I  Let  Costconstanti^)  =  CosUntrai^)  +  Upperintrai^)  +  Lowerintra{^)- 
Then  CosUotali^)  =  Costconstant{^)  +  i?COSf,„ter (H) . 

Proof  of  Lemma  4.1     From  (4.1),  (4.2),  (4.3)  and  (4.4), 

m 

m-1 

«=i  uePi,veP,,i<j 

m-l 

5:        ^       (/„.(z,,u),.,U//(^,«W)-  (4.9) 

Adding  and  subtracting  the  same  value  to  and  from  the  right-hand  side  of  (4.9)  still 
hold  the  equality.  Thus,  by  adding  Upperintrai^)  (4.7)  and  Zoit;er,„<ro(n)  (4.8)  to 
the  first  term  of  the  right-hand  side  of  (4.9),  and  by  subtracting  C/pper,„tra(n)  from 
the  second  term  and  Lowerintra{^)  from  the  third  term,  we  obtain 

Costtotali^)    =  Cost 

constant  (n)  + 

m-l 

i-\  u€Pi,vePj,i<j 

m-l 

E        E       iUv,n){rvl  -  rol)  +  fj{v,u){w}  -  w])).{i.lO) 

i=l  uePi,vePj,i<j 

Hence,  applying  (4.5)  and  (4.6)  to  (4.10)  in  sequence, 

m  — 1 

Costtotali^)     =     Costconstanti'n)  +  "Y  {9{u,  v)  +  g{v ,  u)) 

i=i  ueP,,vePj,i<j 

=     CostconstantO^)  +  B.COStinteri^)- 

□ 

Figure  4.3(a)  shows  the  relative  cost  matrix  driven  from  the  move  and  find  fre- 
quency matrices  in  Figure  4.1  using  (4.5)  if  wf,^  —  wl^  —  2  and  —  w^j  =  1.  For 
the  partition  11  =  {{1, 2},  {3,4},  {5,  6}},  therefore,  the  dual  problem  is  concerned 
with  the  sum  of  all  elements  in  the  shaded  areas,  say,  Rcostinter{^)  =  300.  Since 
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(a)  relative  cost  matrix 


(b)  transformation  process 


(c)  cumulative  cost  matrix 


Figure  4.3.  Constructing  the  Cumulative  Cost  Matrix 

Costconstanti^)  =  Costintra{^)  +  U  pperintraC^)  +  Lowerintra{^)  =  44  +  96  +  97 
237,  Costtotali^)  =  CostconstantCn)  +  RcOsUnter{U)  =  237  +  300  =  537. 


Theorem  4-1  If  H  is  an  optimal  partition  which  minimizes  RcosUnteT,  then  IT  is  also 
an  optimal  partition  which  minimizes  Cost^dah 

Proof  of  Theorem  4-1  By  Lemma 4.1,  it  follows  that  C osttotai{n.)  =  C ostconstant{^)+ 
Rcost inter (ii)-  Note  that  Costconatanti^)  is  Constant  independent  of  how  to  partition 
because  it  consists  of  three  intra-cluster  communication  costs  defined  in  the  equations 
(4.2),  (4.7),  and  (4.8).  Thus,  it  holds  that  11  minimizes  Rcostinter  iff  11  minimizes 

C  OStioial-  □ 

4.1.3    Cumulative  Cost  Matrix 

We  have  the  relative  cost  matrix  A  of  size  n  such  that  A{u,v)  =  g{u,v)  for  all 
{u,v)  £  V  X  V,  as  shown  in  Figure  4.3(a).  The  dual  problem  can  then  be  restated 
as  the  problem  of  partitioning  A  into  submatrices  Ajj  for  i,j  =  1, . . . ,  m  such  that 
the  main  diagonal  submatrices  A,-,  are  square  and 

m  — 1  7n 

RcOStiieri^)  =  E    E  (A*,  +  a;.)  (4.11) 
i-1  j=i+l 

is  minimized,  where  A'^  is  the  sum  of  all  elements  in  A  which  represents  the  relative 
cost  from  the  ith  cluster  to  the  jth  cluster.  In  the  example  of  Figure  4.3(a)  where 
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Pi  =  {1,2},  P2  =  {3,4},  and  P3  =  {5,6},  i2co<,„(n)  =  A;,  +  A*,,  +  A;,+  A*, 
+  A^i  +  AI2  =  51  +  52  +  48  +  50  +  45  +  54  =  300  =  RcosUr^teri^)- 

Because  A*^  and  A^,  vary  for  each  partition,  it  is  desirable  to  obtain  their  values 
directly  without  duplicated  computations.  Thus,  we  transform  the  relative  cost  ma- 
trix A  to  a  cumulative  cost  matrix  B  so  that  for  an  arbitrary  partition  it  is  possible 
to  obtain  the  value  of  the  expression  T,f=i+i{A*j  +  A^,)  of  (4.11)  by  simply  accessing 
to  an  element  in  B  which  is  defined  as: 

B(5,t)   =   Y.EMkJ)    if  •s<i;  (4.12) 

fc=s  /=( 
t-1  n 

B{t,s)   =   J2J2Mi^k)    a  s<t.  (4.13) 

k=s  l=t 

The  B  matrix  can  be  computed  by  accumulating  the  elements  above  the  main 
diagonal  of  A  from  right  to  left  for  each  row  and  then  from  bottom  to  top  for  each 
column,  while  the  elements  below  the  main  diagonal  of  A  are  accumulated  from 
bottom  to  top  for  eeich  column  and  then  from  right  to  left  for  each  row,  as  shown  in 
Figure  4.3(b). 

Lemma  4-2  Let  the  size  of  P,  be  such  that  fi,  —  n.  Given  an  integer  i  for 
1  <  i  <  m  —  1, 

a)  B{s,  t)  =  E7=,+i       if  s<t,  where  s  =  di  +  --  •  +       + 1  and  i  =  </i  -f-  ■  •  ■  +  dj_i  +  1. 

b)  B{t,s)  =  E^i+i  A*,-  if  5  <  f,  where  5  =  ffi  +  -  •  •  +  f/,_i  +  l  and  t  =  di  +  - ■  ■  +  dj_i  +  l. 
Proof  of  Lemma  4-2     The  equation  (4.12)  can  be  divided  into  (m  —  i  +  1)  terms 

k=»  i=t 

t-l  <+cij+i-l  t+dj  +  i+(ij+2-l  n 

=  E(  E  A(^,/)+   y:  A{kj)+...+    y:  Mk,i)). 
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Since  each  term  is  the  sum  of  all  elements  in  A,j,  where  j  =  i  +  . . .  ,m,  the  left-hand 
side  of  a)  can  be  expressed  as 

B{s,t)   =   A*(,+i)  +  A*(,+,)  +  •  •  •  +  A*„ 

=  E  A*. 

i=t+i 

Proving  in  the  same  way  as  a),  it  is  obvious  that  b)  holds.  □ 
Given  the  tth  cluster  Pi,  Lemma  4.2  shows  that  for  all  j  =  i -f  1, . . . ,  m,  B{s,t) 
represents  the  sum  of  the  relative  cost  from  P,  to  Pj,  and  B(^,  s)  represents  the  sum 
of  the  relative  cost  from  Pj  to  P,.  Thus,  B(.s,^)  +  B(/,5)  gives  the  total  relative 
cost  among  them.  Since  the  positions  of  B(s,/)  and  B(t,5)  are  symmetric  over  the 
main  diagonal,  we  combine  them  by  adding  B(i,s)  to  B(5,^)  without  affecting  the 
optimization  goal.  Then,  B(5,i)  =  b^t  >  0  ii  s  <  t  and  B(5,^)  =  0  if  s  >  t.  Hence, 
the  equation  (4.11)  can  be  rewritten  by 

m  — 1 

RcostZ,,,{U)  =  Y:  Kt,  (4.14) 

t=l 

where  s  —  Ylk=i  ^k-i  +  1,  c?o  =  0  and  t  =  Ylk=i  +  1  for  every  i.  In  the  example 
of  Figure  4.3(c),  Rcostf^t„{\l)  =  613  +  635  =  198  +  102  =  300  =  Rcostit„{U).  The 
613  entry  is  the  cost  between  {1,2}  and  {3, 4, 5,  6},  and  the  635  entry  is  that  between 
{3,4}  and  {5,6}.  Thus,  the  sum  of  613  +  ^'35  is  the  total  cost  among  {1,2},  {3,4}, 
and  {5,6}. 

4.2  Algorithm 

4.2.1    Optimal  Structure  and  Overlapping  Subproblems 

Given  the  cumulative  cost  matrix  B  =  (hij),  1  <  i  <  j  <  n,  we  wish  to  find  a  se- 
quence of  (m— 1)  elements  whose  sum  is  minimum  such  that  the  first  element  is  chosen 
among  bij,  and  if  6,j  is  selected  as  the  rth  element,  the  next  possible  element  is  one  of 
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bjkS  foT  j  +  1  <  k  <  n  — (m  — 1— r)  +  l.  We  denote  the  problem  to  yield  this  optimal  se- 
quence as  LPART(1,  n,  m).  Let  <  bij,  bjk, . .  .,bst>  be  an  optimal  sequence  of  (m- 1) 
elements  for  the  problem  LPART(1,  n,  m).  Then,  <  bjk, . . . ,  t^t  >  must  be  an  optimal 
sequence  of  the  subproblem  LPART(j,  n,  m  —  1).  If  it  is  not,  there  is  another  sequence 
<  bji, . . .  ,bxy  >  oi  {m  —  2)  elements  whose  sum  is  less  than  that  of  <  bjk,  ■  ■  ■  ,bst  >■ 
Hence,  <  feij,  bji, . . . ,  i^y  >  is  a  sequence  of  (m  —  1)  elements  with  less  value,  which  is 
a  contradiction.  Thus,  the  LPART  problem  has  an  optimal-substructure  property.  In 
addition  to  the  optimal-substructure  property,  the  LPART  problem  has  overlapping 
subproblems.  LPART(j,  n,p)  {or  1  <  p  <  m  is  referenced  G{j  —  1,  m  — p)  times  during 
the  computations  of  LPART(ji',  n,p -|- 1)  for  j'  =  1,2, ...  ,j  —  1  if  p  <  m  —  1.  For  ex- 
ample, consider  that  n  —  9  and  m  =  5.  LPART(5,9,3)  is  referenced  0(4,2)  =  (i)  =  3 
times  during  the  computations  of  LPART(4,9,4),  LPART(3,9,4),  and  LPART(2,9,4) 
which  correspond  to  the  way  of  partitioning  vertices  1  through  4  into  {{1,2,3},  {4}}, 
{{1, 2},  {3, 4}},  and  {{1},  {2, 3, 4}},  respectively.  Thus,  many  other  subproblems 
share  subsubproblems.  Since  an  optimal  solution  has  optimal  structure  and  over- 
lapping subproblems,  dynamic  programming  should  be  applied  to  the  optimization 
problem. 

4.2.2    Recursive  Definition 

We  now  define  the  value  of  an  optimal  solution  recursively.  Let  ^j(p)  be  the 
optimal  value  of  the  problem  when  a  sequence  {j,j  +  l,...,n}  of  vertices  is  linearly 
partitioned  into  p  clusters.  Then,  we  wish  to  find  (f)\{m).  The  optimal  cost  <^j(p)  can 
be  defined  as  follows.  If  p  =  1,  vertices  j  through  n  are  grouped  into  one  cluster,  so  no 
inter-cluster  cost  is  necessary.  If  p  >  2,  each  bjk  for  j  -|- 1  <  k  <  n  —  p4-2  is  considered 
as  the  first  element,  representing  the  cost  between  the  first  cluster  {j,j-\-\,...,k  —  \] 
and  the  remaining  (p— 1)  clusters  consisting  of  vertices  k  through  n.  Thus,  the  optimal 
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cost  of  partitioning  {j,j  +  l,...,n}  into  p  clusters  is  the  minimum  value  among  the 
sums  of  bjk  and  (f>k{p  —  I)  for  all  possible  k.  The  dynamic  programming  formulation 
can  be  written  as 

for  l<j<n  —  From  now  on,  we  assume  that  each  of    and  A:  has  the  same 

range  as  in  this  section  if  it  is  not  specified  explicitly. 

4.2.3    Computing  Optimal  Costs 

Based  on  the  recurrence  relation  (4.15),  we  present  an  algorithm  to  compute  the 
optimal  cost  by  using  a  bottom-up  approach.  The  following  pseudocode  inputs  the 
cumulative  cost  matrix  B  =  for  I  <  i  <  j  <  n  and  the  number  of  clusters 
m.  During  the  execution  of  the  procedure  for  each  p,  1  <  p  <  m,  the  pseudocode 
uses  the  main  diagonal  element  bjj  for  storing  the  optimal  cost  <f>j{p)  and  the  lower 
triangle  element  6(„_p^2)j  for  recording  which  index  of  k  achieved  the  optimal  cost 

Mp)- 

LPART(1,  n,m) 

1  for  j  i—  1  to  n 

2  bjj  ^  0 

3  for  p  <—  2  to  m 

4  for  j  i—  n  —  p  +  1  downto  1 

5  begin 

6  bjj  <—  oo 

7  for  k*— j  +  lton  —  p  +  2 

8  begin 

9  sum  <-  bjk  +  bkk 

10  if  sum  <  bjj 

11  begin 

12  bjj  <—  sum 

13  6(„_p+2)j  <—  k 

14  end 

15  end 

16  end 
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Figure  4.4.  Computation  of  the  LPART(l,6,p)  Problem  for  1  <  p  <  6 


The  algorithm  first  computes  bjj  0  for  j  =  1, . . . ,  n  in  lines  1-2,  each  represents 
the  minimum  cost  of  partitioning  {j, . . . ,  n}  into  one  cluster,  <f>j{l)-  It  then  uses  the 
recurrence  relation  (4.15)  to  compute  the  minimum  cost  of  partitioning  {j, . . . ,  n}  into 
two  clusters,  <^j(2),  during  the  first  execution  of  the  loop  in  lines  3-16,  and  so  forth. 
At  each  execution  of  the  loop,  the  optimal  cost  is  computed  for  1  <  ^  <  n  —  p  -f- 1, 
assuming  that  the  size  of  a  cluster  is  at  least  1.  For  a  choice  of  p  and  j,  the  bjj  cost 
computed  in  lines  5-16  depends  on  the  sum  of  the  bjk  cost  and  the  bkk  cost  which 
Wcis  already  computed  from  the  previous  execution  of  the  loop  in  lines  3-16.  Once  k 
is  determined,  it  is  recorded  at  the  b(^n-p+2)-i  entry. 

Figure  4.4  illustrates  this  procedure  for  1  <  p  <  6  when  n  =  6.  The  shaded  entries 
in  the  table  for  a  particular  p  are  computed  from  the  upper  triangle  and  the  main 
diagonal  of  the  table  for  (p  —  1).  Note  that  the  computation  is  practically  achieved 
using  only  one  table  and  each  separate  table  is  used  to  show  the  execution  step  cis 
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p  increases.  The  upper  triangle  containing  tlie  cumulative  costs  among  vertices  does 
not  vary  and  the  lower  triangle  is  initially  set  to  0.  When  p  =  1,  every  main  diagonal 
element  is  set  to  0.  For  the  case  oi  p  ^  1,  consider  how  to  compute  the  in  entry 
when  p  =  3,  which  is  equivalent  to  </>i(3). 

<^i(3)   =   min{6i2  +  ^2(2),foi3  +  '^3(2),?h4  +  <^4(2),6i5  +  '^5(2)} 

From  the  table  for  p  =  2, 

<j)i{3)    =   min{6i2  +  &22,  ^i3  +  ^33,  ^14  +  ^44,  ^>i5  +  ^55} 
=  min{118  +  98, 198  +  71, 220  +  49, 199  +  27} 
=  216 

when  k  =  2.  Thus,  bu  =  216  and  651  =  2. 

4.2.4    Constructing  an  Optimal  Partition 

Since  LPART(1, n,m)  only  gives  the  optimal  cost  (l)i{m),  we  need  to  construct  an 
optimal  set  of  m  clusters  from  the  recorded  k  information  in  the  lower  triangle.  Let 
<  ^1,  ^2,  •  •  • )  km-i  >  be  a  sequence  of  values  of  k  for  an  optimal  set  of  m  clusters.  The 
6(„_m+2)i  entry  gives  the  ki  value,  and  the  ?>(„_(„t_i)+2)A:i  entry  gives  the  k2  value,  and 
the  6(n-(m-2)+2)Jt2  entry  gives  the  ks  value,  and  so  forth.  Once  <  ki,k2, . . . ,  k^-i  > 
is  obtained,  we  can  easily  divide  the  vertex  set  into  m  clusters,  {l,...,A;i  —  1}, 
{ki, . . .  ,k2  —  1},  . . . ,  {km-i, . . . , n},  respectively.  This  is  computed  by  the  following 
pseudocode  LPART-CLUSTER(1,  n,  m). 

LPART-CLUSTER(1,  n,  m) 

1  LPART(1,  n,m) 

2  CLUSTER  ^fH 

3  ko  <—  1 

4  i      n  —  m  +  2 

5  j^l 
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6  for  /  <—  1  to  m  —  1 

7  begin 

8  k  i-  bij 

9  CLUSTER  <-  CLUSTERS  {ki_i,...,k,  -  1} 

10  i 

11  j  ^  k 

12  end 

13  CLUSTER  ^  CLUSTER  U  {k,  ...,n} 

In  the  example  of  Figure  4.4,  the  call  LPART-CLUSTER(1,6,3)  first  refers  to  the 
entry  651  which  contains  ki  =  2.  Then,  the  subsequent  value  of  k  is  obtained;  k2  =  3. 
Hence,  an  optimal  partition  of  3  clusters  is  {{1},  {2},  {3,4,5,6}}. 

Once  we  have  computed  the  cumulative  cost  matrix  in  O(n^),  the  procedure 
LPART  computes  the  optimal  cost  in  0{mn^).  The  other  steps  in  the  loop  of  lines 
6-12  require  0{m).  Hence,  the  total  time  complexity  is  O(mn^). 

4.2.5    Constraint  on  the  Cluster  Size 

Let  us  now  consider  the  constraint  imposed  on  the  size  of  a  cluster.  Denote  s  as 
the  maximum  size  of  a  cluster  allowed  in  the  process  of  partitioning  such  that  ms  >  n. 
The  following  pseudocode  MLPART(1,  n,  m,  5)  is  obtained  by  simply  adjusting  the 
k  values  in  line  9  of  the  LPART(1,  n,m)  procedure.  The  lowest  k  value  is  changed 
from  (i  + 1)  to  max(j  +  l,n  —  (p— 1)5-|-1)  and  the  highest  is  changed  from  {n  —  p  +  2) 
to  min(n  —  p  +  2,j  +  s)  in  order  to  generate  the  valid  k  values  for  the  constraint. 

MLPART(l,n,m,5) 

1  if  ms  <  n 

2  return 

3  for  j  <—  1  to  n 

4  bjj  <-  0 

5  for  p  <—  2  to  m 

6  for  j  <—  n  —  p  +  1  downto  1 

7  begin 

8  bjj  <—  00 

9  low  <—  max{j  +  1,  n  —  (p  —  1)5  +  1) 
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Figure  4.5.  Computation  of  the  MLPART(l,6,p,3)  Problem  for  1  <  p  <  6 
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high  <—  min(n  —  p  +  2,  j -f  s) 
if  low  <  high 

for  k  *—  low  to  high 
begin 

sum  <—  hjk  +  hkk 
if  sum  <  bjj 
begin 


sum 


end 


end 


end 


Figure  4.5  illustrates  the  MLPART(l,6,p,3)  problem  for  1  <  p  <  6.  The  indices 
of  the  lightly  shaded  area  for  a  particular  p  represent  the  valid  k  values  which  are 
used  for  the  next  step  of  (p  +  1),  satisfying  that  the  size  of  any  cluster  in  a  resulting 
optimal  partition  is  equal  to  or  less  than  3.  For  example,  consider  the  first  row  of 
the  table  for  p  =  1  which  is  used  to  compute  ^i(2)  recorded  in  the  bu  element  of  the 
table  for  p  =  2.  Unlike  in  the  example  of  Figure  4.4,  only  the  614  element  is  used  to 
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compute  (f>i{2)  because  the  consideration  of  other  elements  results  in  a  cluster  with 
the  size  greater  than  3.  It  means  that  the  low  variable  is  set  to  n  —  (p  —  1)5  +  1  = 
6-(2-l)x3  +  l=  4  and  the  high  variable  is  set  to  j  +  5  =  1  +  3  =  4.  Thus, 
?!>i(2)  =  614  +  644  =  220  +  0  =  220.  Hence,  the  optimal  partition  is  {{1, 2, 3},  {4, 5, 6}} 
if  p  =  2.  Note  that  for  each  p  >  4,  the  optimal  cost  (f>i{p)  and  partition  11  are  the 
same  as  the  example  of  Figure  4.4  because  the  constraint  does  not  limit  the  range  of 
the  k  values.  It  always  holds  when  n  —  (p  —  1)  <  s. 

Given  the  maximum  size  of  a  cluster  s  and  an  integer  K  which  is  the  maximum 
communication  cost  allowed  for  the  entire  system,  the  MLPART  procedure  can  also 
be  used  to  generate  all  possible  optimal  partitions  for  a  sequence  of  n  nodes.  In  the 
example  of  Figure  4.5,  consider  the  case  when  n  =  6,  s  =  3,  and  K  —  500.  Because 
the  MLPART  procedure  optimizes  the  relative  cost,  we  first  compute  the  relative  cost 
K'  =  K  —  Cost  constant  =  500  —  237  =  263.  Then,  the  MLPART  procedure  computes 
(t>i{p)  for  2  <  p  <  n  until  (^i(p)  >  K' .  Because  (^i(2)  =  220  and  <^i(3)  =  257  are 
less  than  K\  the  possible  optimal  partitions  are  {{1, 2,  3},  {4, 5, 6}}  if  p  =  2  and 
{{l},{2,3},{4,5,6}}ifp  =  3. 

4.3    Chapter  Summaiy 

In  highway  cellular  systems,  a  linear  array  of  n  heterogeneous  base  stations  may 
generate  multiple  types  of  traffic  among  themselves.  We  have  considered  the  problem 
of  finding  a  set  of  disjoint  clusters  to  cover  n  base  stations  such  that  it  minimizes 
the  total  communication  cost  for  the  entire  system.  The  total  communication  cost 
consists  of  the  costs  incurred  by  two  different  types  of  traffic,  mobility  management 
data  and  user  data.  For  each  type  of  traffic,  the  cost  of  intra-cluster  communication 
is  usually  lower  than  that  of  inter-cluster  communication.  The  optimal  partitioning 
problem  is  solvable  in  polynomial  time  of  0{mn^)  by  dynamic  programming  for 
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an  arbitrary  number  of  clusters  and  size  of  a  cluster,  where  m  is  the  number  of 
clusters  in  a  partition.  In  addition,  the  algorithm  finds  all  valid  partitions  in  the 
same  polynomial  time  if  given  a  constraint  on  the  cluster  size  and  the  total  allowable 
communication  cost  for  the  entire  system. 

The  optimal  deployment  of  virtual  cells  in  a  hexagonal  mesh  arrangement  of  n 
base  stations  is  described  in  the  following  chapter. 


CHAPTER  5 

MULTIWAY  PARTITIONING  IN  HEXAGONAL  CELLULAR  SYSTEMS 

Given  a  hexagonal  mesh  of  base  stations  in  cellular  systems  we  consider  the  prob- 
lem of  finding  a  cover  of  disjoint  clusters  of  base  stations  which  generate  multiple 
types  of  traffic  among  themselves.  The  problem  differs  from  general  graph  parti- 
tioning problems  in  that  it  considers  not  only  communication  costs  but  also  the 
underlying  topology  among  base  stations.  The  objective  is  to  minimize  the  total 
communication  cost  for  the  entire  system  where  inter-cluster  communication  is  more 
expensive  than  intra-cluster  communication  for  each  type  of  traffic.  The  problem  is 
transformed  into  the  dual  based  on  a  topology  matrix  and  a  relative  cost  matrix.  We 
develop  several  heuristics  for  the  dual.  These  heuristics  produce  optimal  partitions 
with  respect  to  the  initial  partition,  based  on  the  techniques  of  moving  or  interchang- 
ing the  boundary  nodes  between  adjacent  clusters.  The  heuristics  are  compared  and 
shown  to  behave  quite  well  through  experimental  tests  and  analysis. 

5.1    Problem  Formalization 

5.1.1    Labeling  Scheme 

The  physical  deployment  of  n  base  stations  is  represented  by  a  planar  graph  of 
a  hexagonal  mesh  of  n  base  stations,  where  the  vertices  represent  base  stations  and 
the  edges  represent  the  adjacency  of  base  stations.  The  vertices  on  the  exterior  face 
of  the  graph  are  level  1  vertices  and  the  vertices  on  the  exterior  face  of  the  subgraph 
induced  by  removing  level  1  vertices  are  level  2  vertices,  and  so  on.  A  planar  graph 
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(a)  The  (column,  row)  indexing  (b)  Tlie  two-dimcaisional  labeling 


Figure  5.1.  Labeling  of  the  Physical  Topology  for  H3 

is  s-outer planar  if  it  has  no  vertices  of  level  greater  than  s.  Denote  Hg  as  an  s- 
outerplanar  graph  of  a  hexagonal  mesh.  It  is  then  known  that  the  number  of  vertices 
in  an  is  n  =  3s^  —  35  -f  1  and  the  number  of  columns  in  each  of  three  directions 
IS  d  =  2s  —  I.  Before  formalizing  the  problem,  it  is  necessary  to  introduce  a  labeling 
scheme  to  describe  the  problem  mathematically. 

Starting  from  the  left-most  column  of  an  Hg,  each  column  is  indexed  from  0 
through  J  —  1  in  sequence.  Then  the  bottom  vertices  of  every  column  constitute  row 
0,  the  next  vertices  of  every  column  row  1,  and  so  forth.  Once  the  column  index  c 
and  the  row  index  r  of  a  vertex  is  determined,  the  vertex  is  labeled  such  that 
i  =  2c  and  j  =  2r.  In  this  two-dimensional  coordinate  system  of  an  Hg,  a  point  (i,  j) 
for  0  <  I,  j  <  2[d  —  1)  can  represent  either  a  vertex  Vij  if  i  and  j  are  even  or  an  edge 
e,j  otherwise.  Figure  5.1  illustrates  the  labeling  of  the  physical  topology  for  an  Hz 
of  n  =  19  nodes. 

If  given  a  vertex  at  most  six  edges,  each  leading  to  an  adjacent  vertex,  are 
directly  identified  by  the  labeling  scheme  as  follows: 

•  ^tO-i-i),  e(,+i)(j+i),  e(,+i)_,-,  e,(j_i),  e(i_i)(,_i),  and  e(i_i)j  if  i  <  J  -  1; 

•  e(,.,.i)(j_i),  e,(j_i),  e(i_i)(j_i),  and  e(,_i)j  \ii  =  d  -  l\ 

•  e»(i+i)>  e(,+i)(_,_i),  e,(j_i),  e(i_i)j,  and  e(,_i)(j+i)  \i  i  >  d  -  I. 
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On  the  other  hand,  if  given  an  edge  e,j,  the  two  vertices  connected  by  the  edge 
are  directly  identified  by  the  labeling  scheme  as  follows: 

•  ^(.-1);  and  ^(,+1);  if  i  is  odd  and  j  is  even; 

•  Vi^j-i)  and  if  i  is  even  and  j  is  odd; 

•  and  if  both  i  and  j  are  odd  and  i  <  d; 

•  and  t;(,_i)(j+i)  if  both  i  and  j  are  odd  and  i  >  d. 

In  the  example  of  Figure  5.1(b),  a  vertex  V24  leads  to  the  six  edges  613,614,623,625, 
634,  and  635,  which  are  connected  to  the  neighboring  vertices,  Uo2)i'o45f22)'^26)V44,  and 
^46,  respectively. 

5.1.2    Topology  Matrix 

To  handle  the  underlying  topology  constraint  on  clustering,  we  construct  a  topol- 
ogy matrix  T  =  for  0  <  i,j  <  2{d  —  1),  where  tij  corresponds  to  either  a 
vertex  u.j  if  i  and  j  are  even  or  an  edge  e,j  otherwise  in  the  labeling  scheme.  If 
Vij  currently  belongs  to  cluster  P^,  the  element  tij  is  set  to  x.  If  the  two  vertices 
connected  by  e,j  belong  to  clusters  and  Py,  the  element  Uj  is  set  to  xy  (ov  x  <  y. 
Figure  5.2(b)  shows  the  topology  matrix  T  derived  from  an  initial  partition  of  3 
clusters  in  Figure  5.2(a). 

From  the  topology  matrix  T  ,  the  indices  of  the  boundary  edges  between  a 
pair  of  adjacent  clusters  P^  and  Py  can  be  represented  by  a  list  of  edge  elements 
Lxy  =  {e.jl^.j  =  xy,x  <  y}.  Using  the  list  L^y,  the  boundary  vertices  between 
them  can  be  directly  identified  by  the  labeling  scheme  itself.  Denote  V^y  and  V^y  as 
the  indices  of  the  boundary  vertices  for  P^  and  Py,  respectively.  In  the  example  of 
Figure  5.2,  L12  =  {603,613,623,633,643},  Vj^  =  {'{^04,^^24,^^44},  and  V^^  -  {^02,^22,^42}- 
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(a)AF«ititi(»of3clusten  (b)  The  topology  imtlixT 


Figure  5.2.  Topology  Matrix  for  a  Partition  of  H3 

5.1.3    Relative  Cost  Matrix 

The  communication  among  base  stations  is  considered  a.s  a  full  mesh  of  point-to- 
point  logical  network  to  represent  a  possible  communication  between  mobile  hosts 
through  different  base  stations.  The  communication  network  is  described  by  a  com- 
plete directed  graph  G  =  {V,E),  where  \V\  =  n.  The  vertices  of  the  graph  represent 
base  stations  and  the  edges  represent  directional  communication  links  between  base 
stations.  Each  edge  is  assigned  a  move  frequency  by  a  function  /m  :  x  V  — * 
and  a  find  frequency  hy  ff  :  V  x  V  ^  i?"*".  Denote  lu}^  and  w'^  as  the  weight  of  a 
move  operation  within  a  cluster  and  between  clusters,  respectively,  and  luj  and 
as  that  of  a  find  operation  within  a  cluster  and  between  clusters,  respectively.  Then 
we  define  a  relative  cost  function  c  :  V  x  V  — >  T?"*"  for  all        Vki)  E  V  x  V  as 

c{Vij,  Vkl)  =  fm{Vij,  Vkl){wl^  -  +  ff{v,j,Vkl){w'^f  -  w]),  (5.1) 

where  (lu^  —  ly^)  is  the  relative  weight  of  a  move  operation  and  {w'j  —  is  the 
relative  weight  of  a  find  operation.  The  cost  c{vij,vi;i)  represents  the  total  relative 
cost  of  fm{vij,Vki)  move  and  ff{vij,Vki)  find  operations  from  Vij  to  Vki  if  they  are  in 
different  clusters.  Denote  c{vij,Vki)  as  c^j-ki-  The  relative  cost  matrix  C  =  {cij.^kl) 
for  0  <  i,  j,  /  <  2{d  —  1)  is  derived  by  the  equation  (5.1)  from  the  move  and  find 
frequencies  among  n  vertices. 
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5.1.4    Dual  Problem 

Let  n  =  {Pi, . . .  ,Pm}  be  a  partition  of  m  clusters  of  contiguous  vertices  such 
that  n  Pj,  =  0  and  U^Pr  =  V  for  x,  y  =  1, . . . ,  m  and  x  ^  y.  The  intra-cluster 
communication  cost  for  all  {vij,Vki)  E:  V  x  V  is 

Costintra  =  ^         {fm  {vij  ,  +  //  {vij ,  Vkl)w)) 

and  the  relative  inter-cluster  communication  cost  of  11  is 

m— 1 

RcOStintcTi^)     =     X)  m  (^•'J;'^'  +  Ckl-ij). 

X=l  VijePx,VktePy,l-<y 

The  total  communication  cost  of  11  is  then 

Because  Costintra  is  constant  independent  of  how  to  partition,  it  holds  that  11 
minimizes  Costtotai  iff  H  minimizes  Rcostinter  ■  Let  us  define  the  relative  intra-cluster 
communication  cost  of  11  as 

m 

RcOSti„tra{n)  =  X       X      [Cij-kl  +  Cfc/jij). 
x-1  ii,j,i)jfci€Pi 

Because  the  sum  of  Rcostintra  and  Rcostinter  constant,  our  task  of  finding  an 
optimal  partition  11  which  minimizes  Rcostinter  is  equivalent  to  that  of  finding  an 
optimal  partition  11  which  maximizes  Rcostintra- 

Hence,  the  dual  problem  is:  given  the  topology  matrix  T  =  and  the  rel- 
ative cost  matrix  C  =  (c,j;jt;)  for  a  system  of  n  hexagonal  cells,  find  a  cover  of  m 
disjoint  clusters  of  contiguous  base  stations  H  =  {Pi,...,Pm},  so  as  to  maximize 

RcOStir,tra{^)- 

5.2  Heuristics 

We  consider  heuristics  for  the  problem:  starting  with  an  arbitrary  partition  H  of  m 
clusters,  we  try  to  increase  the  initial  relative  intra-communication  cost  Rcostintrai'^) 
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by  repeated  applications  of  a  two-way  optimization  procedure  to  pairs  of  adjacent 
clusters.  In  the  two-way  optimization  procedure,  we  try  to  increase  the  sum  of  the 
initial  relative  intra-communication  cost  of  each  cluster  by  moving  or  interchanging 
boundary  nodes  until  no  further  improvement  is  possible.  The  resulting  pair  of  adja- 
cent clusters  are  then  pairwise  optimal  with  respect  to  the  initial  partition.  Because 
the  two-way  optimality  for  a  pair  of  adjacent  clusters  may  affect  that  for  another 
pair  of  adjacent  clusters,  more  than  one  pass  through  all  pairs  of  adjacent  clusters 
may  be  required.  This  process  can  be  repeated  with  another  initial  partition  11  of  m 
clusters,  and  so  on,  so  as  to  obtain  as  many  locally  maximum  partitions  as  we  desire. 
Then,  one  of  the  resulting  partitions  has  a  fairly  high  probability  of  being  a  globally 
maximum  partition. 

5.2.1    Interchanging  or  Moving  Boundary  Nodes 

Consider  a  pair  of  adjacent  clusters  and  P,,  in  the  two-way  optimization  pro- 
cedure. The  interchange  of  a  pair  of  boundary  nodes,  u,j  G  V^y  and  Vki  €  Vxy,  is  said 
to  be  feasible  if  the  resulting  clusters  preserve  the  underlying  topology  constraint.  In 
other  words,  all  nodes  in  P^  {Py)  adjacent  to  Vij  {vki)  must  be  connected  in  their 
physical  topology  before  the  interchange  and  Vij  (vkt)  must  be  connected  to  at  least 
one  of  the  nodes  in  Py  (Px)  after  the  interchange. 

To  determine  feasible  pairs  of  interchanging  nodes,  we  use  the  topology  matrix 
T  where  an  element  tij  corresponds  to  a  node  Vi^  or  an  edge  e,j.  Given  a  node 
Vij  €  Vxy,  at  most  six  edges  adjacent  to  the  node  can  be  directly  identified  by  the 
labeling  scheme  itself.  Denote  Adj{vij)  as  an  ordered  list  of  those  edges  which  are 
sequentially  arranged  in  a  circular  fashion.  Then  it  is  said  that  the  node  v,j  can  be 
moved  into  Py  for  the  interchange  when  the  following  two  conditions  are  satisfied: 
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Figure  5.3.  Interchanging  Boundary  Node  Elements 

1.  There  is  only  one  continuous  subsequence  of  edges  which  are  equal  to  xx  in 
Adj{vij). 

2.  There  is  at  least  one  edge  in  Adj[vij)^  which  is  equal  to  xy  and  does  not  lead 
to  the  node  Vki  €  Vxy. 

Assume  that  initially  all  nodes  of  each  cluster  are  connected  in  their  physical 
topology.  Condition  1  implies  that  the  remaining  nodes  of  after  removing 
still  preserve  the  connectivity  in  their  physical  topology.  Condition  2  implies  that 
the  removed  node  u.j  is  also  connected  to  its  adjacent  cluster  Py  whose  nodes  are 
already  connected  in  their  physical  topology.  In  the  same  way,  the  node  Vki  must  also 
satisfy  the  above  two  conditions  so  that  it  can  be  moved  into  P^  for  the  interchange. 
Therefore,  the  interchange  of  a  pair  of  boundary  nodes  Vij  6  Vxy  and  Vki  €  Vxy  is 
feasible  iff  both  and  vm  satisfy  Condition  1  before  the  interchange  and  Condition 
2  after  the  interchange. 

In  the  example  of  Figure  5.2(b),  we  consider  interchanging  V22  €  P2  and  U44  6 
Pi.  For  node  ^22,  since  the  only  one  continuous  subsequence  {612,611,621,632}  in 
Adj{y22)  =  {623,612,611,621,632,633}  is  equal  to  22,  Condition  1  holds.  In  addition, 
since  623  =  12  and  it  does  not  lead  to  V44,  Condition  2  also  holds.  Thus,  V22  can  be 
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moved  to  Pi  for  the  interchange.  At  the  same  time,  for  node  V44,  since  the  only  one 
continuous  subsequence  {645,634}  in  Adj{v44)  =  {634,633,643,653,654,645}  is  equal  to 
11,  Condition  1  holds.  In  addition,  since  643  =  12  and  it  does  not  lead  to  V22, 
Condition  2  also  holds.  Thus,  V44  can  be  moved  to  P2  for  the  interchange.  Hence,  V22 
and  V44  are  a  feasible  pair  of  interchanging  nodes  between  Pi  and  P2.  Figure  5.3(a) 
shows  the  resulting  topology  matrix  T  after  the  interchange.  On  the  other  hand,  the 
nodes  V42  €  P2  and  ^44  G  Pi  cannot  be  interchanged,  as  depicted  in  Figure  5.3(b). 
The  node  V42  is  isolated  after  the  interchange  due  to  the  fact  that  no  edge  except  643 
adjacent  to  V42  is  equal  to  12  in  Figure  5.2(b). 

It  should  be  noted  that  the  move  of  a  boundary  node  u.j  €  Vxy  into  Py  is  rather 
simple  than  the  interchange  because  we  need  to  check  that  only  all  nodes  in  Px 
adjacent  to  f,j  are  connected  in  their  physical  topology  before  the  move. 

5.2.2    Computing  Gains 

When  interchanging  boundary  nodes.  Assume  that  a  pair  of  boundary  nodes, 
Vij  €  Vxy  and  Vki  G  V^y,  is  feasible  for  the  interchange.  Define  the  internal  cost  of 
node  V{j  with  respect  to  P^;  for  the  interchange  to  be 

and  the  external  cost  of  node  Vij  with  respect  to  Py  for  the  interchange  to  be 

Ee{vij)  =  ^  {Cij-k'l'  +  Ck'l'-ij). 

^k'l'^Py^kii^'^k'ii 

Similarly,  we  define  /e(ujt/)  and  Ee{vki)  for  the  boundary  node  Vki.  If  u,j  and  Vki  are 
interchanged,  then  the  gain  g  of  the  increase  in  cost  is  given  by 

g  =  {E,{vij)  +  E,{vk,))  -  (7e(v,,)  +  h{vki)). 
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When  moving  boundary  nodes.  Assume  that  a  node  u,j  €  V^y  is  feasible  for  the 
move.  Define  the  internal  cost  of  node  u,j  with  respect  to      for  the  move  to  be 

Im{Vij)  =         J2        i^ij'i'j'  + 
and  the  external  cost  of  node  Vij  with  respect  to  Py  for  the  move  to  be 

JEmivij)  =      ^    {Cij-k'V  +  Ck'l'-ij)- 

If  Vij  G  Px  is  moved  into  Py,  then  the  gain    of  the  increase  in  cost  is  given  by 

g  =  Emi'Vkl)  -  Im{Vxj). 

5.2.3    Heuristic. 1 

Given  a  pair  of  adjacent  clusters  Px  and  Py,  Heuristic.l  interchanges  only  one 
feasible  pair  of  boundary  nodes  with  a  positive  maximum  gain  g  between  the  two 
clusters.  This  process  is  repeated  with  an  updated  boundary  until  no  feasible  pair 
produces  a  positive  gain.  Let  k  be  the  number  of  feasible  pairs  interchanged.  Then 
the  the  total  gain  with  respect  to  the  sum  of  the  initial  costs  for  the  two  clusters 
is  Gxy  =  I2v=i  9i-  Note  that  a  pair  of  nodes  interchanged  in  the  previous  step  is 
not  interchanged  again  at  the  next  step  because  the  positive  maximum  gain  in  the 
previous  step  becomes  the  minimum  gain  with  the  same  negative  value  in  the  next 
step. 

Heuristic.  l(n) 

1    for  every  pair  of  adjacent  clusters  P^  and  Py  in  H 


2  begin 

3  Determine  the  boundary  node  lists  Vxy  and  V^y 

4  forever 

5  begin 

6  Determine  a  set  of  feasible  pairs 

7  if  there  is  no  feasible  pair 

8  break 
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9  Compute  gains  for  all  feasible  pair 

10  if  there  is  no  feasible  pair  with  a  positive  gain 

11  break 

12  Choose  a  feasible  pair  {vij,Vki)  with  the  maximum  gain 

13  Interchange  v,j  and  Vki 

14  Update  the  boundary  node  lists  V^y  and  V^y 

15  end 


16  end 

5.2.4  Heuristic.2 

Once  a  set  of  feasible  pairs  of  boundary  nodes  is  determined  from  the  boundary 
node  lists  V^y  and  V^y,  Heuristic.2  interchanges  all  feasible  pairs  with  positive  gains 
before  updating  the  boundary  node  lists  V^y  and  K;7-  Initially,  all  feasible  pairs  are 
unmarked.  A  feasible  pair  with  the  maximum  positive  gain  is  first  interchanged  and 
then  all  remaining  unmarked  feasible  pairs  which  involve  the  boundary  nodes  of  the 
interchanged  pair  are  marked  so  that  in  the  next  step  they  cannot  be  considered 
again.  After  marking,  the  gains  for  all  unmarked  feasible  pairs  are  computed  again 
and  then  an  unmarked  feasible  pair  with  the  largest  positive  gain  is  next  interchanged. 
Heuristic.2(n) 

1    for  every  pair  of  adjacent  clusters  Poo  and  Py  in  IT 


2  begin 

3  Determine  the  boundary  node  lists  V^y  and  V^y 

4  forever 

5  begin 

6  Determine  and  unmark  a  set  of  feasible  pairs 

7  if  there  is  no  feasible  pair 

8  break 

9  Compute  gains  for  all  feasible  pairs 

10  while  there  are  unmarked  feasible  pairs  with  positive  gains 

11  begin 

12  Choose  a  pair  {vij,Vki)  with  the  largest  gain 

13  Interchange  Vij  and  v^i 

14  Mark  all  unmarked  feasible  pairs  involving  Vij  or  Vki 

15  Compute  gains  for  all  unmarked  feasible  pairs 

16  end 
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17  Update  the  boundary  node  lists  V^y  and  Vxy 

18  end 

19  end 

5.2.5  Heuristic.3 

Heuristic. 3  is  ba^ed  on  the  observation  that  infeasible  pairs  of  boundary  nodes 
might  become  feaisible  pairs  after  interchanging  feasible  pairs.  For  example,  consider 
a  pair  of  adjacent  clusters  Pi  and  P2  in  the  topology  matrix  of  H3,  where  Pi  = 
{^00,^02,^04}  and  p2  =  {^^20,  "22,^24,1^26}-  Then  V^^  =  Pi  and  =  P2-  A  pair  of 
boundary  nodes  (uo2,  V22)  is  infeasible  for  the  interchange  because  vq2  and  V22  violate 
Condition  1.  However,  if  a  feasible  pair  (^04,  ^20)  is  interchanged,  the  infeasible  pair 
("02,^22)  becomes  a  feasible  pair  for  the  next  interchange. 

The  algorithm  of  Heuristic.3  is  obtained  by  slightly  modifying  that  of  Heuristic. 2. 
Before  performing  line  15  in  Heuristic. 2,  infeasible  pairs  of  boundary  nodes  which 
become  feasible  pairs  due  to  the  interchange  of  nodes  Vij  and  Vki  in  line  13  are  added 
to  the  set  of  unmarked  feasible  pairs. 

5.2.6  Heuristic. 4 

While  the  previous  three  heuristics  interchange  boundary  nodes  between  a  pair 
of  adjacent  clusters,  Heuristic.4  moves  a  boundary  node  in  a  cluster  into  the  other 
cluster  between  a  pair  of  adjacent  clusters.  Since  the  move  of  a  boundary  node 
between  a  pair  of  adjacent  clusters  changes  their  cluster  sizes,  the  constraint  on  the 
cluster  size  is  given  by  the  minimum  and  maximum  cluster  sizes. 

Given  a  pair  of  adjacent  clusters,  Heuristic.4  determines  a  feasible  boundary  node 
with  a  maximum  positive  gain  for  each  cluster  and  compares  them  to  determine  a 
boundary  node  to  be  moved.  Once  the  boundary  node  to  be  moved  is  determined,  it 
is  removed  from  its  cluster  and  added  to  the  other  cluster  as  long  as  its  move  does 
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not  violate  the  constraint  on  the  cluster  size.  This  process  is  repeated  until  there  is 

no  feasible  boundary  node  for  the  move  or  one  of  the  clusters  has  the  minimum  or 

maximum  cluster  size. 
Heuristic. 4(n) 

1    for  every  pair  of  adjacent  clusters  Px  and  Py  in  11 


2  begin 

3  Determine  the  boundary  node  lists  V^y  and  Vxy 

4  while  the  constraint  on  the  cluster  size  is  satisfied 

5  begin 

6  Determine  a  set  of  feasible  nodes 

7  if  there  is  no  feasible  node 

8  break 

9  Compute  gains  for  all  feasible  nodes 

10  if  there  is  no  feasible  node  with  a  positive  gain 

11  break 

12  Choose  a  feasible  node  Vij  with  the  maximum  gain 

13  Move  Vij  into  the  other  cluster 

14  Update  the  boundary  node  lists  V^y  and  Vxy 

15  end 

16  end 


5.2.7    Heuristic. 5.  Heuristic.6,  and  Heuristic.7 

Heuristic. 5,  Heuristic.6,  and  Heuristic.7  are  based  on  a  repeated  application  of 
two-phase  optimization.  In  the  first  phase,  the  heuristics  use  Heuristic.  1,  Heuristic.2, 
and  Heuristic. 3  for  interchanging  the  boundary  nodes  between  adjacent  clusters,  re- 
spectively. Their  second  phase  uses  Heuristic. 4  for  moving  the  boundary  nodes  be- 
tween adjacent  clusters  to  improve  the  gains  obtained  by  their  first  phase  if  possible. 
The  improvement  by  their  second  phase  implies  the  change  of  the  boundary  nodes 
and  the  possibility  that  their  first  phase  further  improves  their  gains.  Thus,  the  two- 
phase  optimization  is  repeatedly  applied  until  no  improvement  is  possible  by  their 
second  phase.  The  basic  idea  of  the  heuristics  is  that  changing  the  cluster  size  by 
their  second  phase  might  overcome  the  limited  improvement  due  to  preserving  the 
cluster  sizes  of  the  initial  paxtition  by  their  first  phase. 
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5.2.8    Heuristics.  Heuristic.9.  and  Heuristic.  10 

Heuristic. 8,  Heuristic.9,  and  Heuristic.  10  are  also  based  on  a  repeated  application 
of  two-phase  optimization.  However,  in  the  first  phase,  the  heuristics  use  Heuristic.4 
for  moving  the  boundary  nodes  between  adjacent  clusters.  Their  second  pha^e  uses 
respectively  Heuristic.  1,  Heuristic. 2,  and  Heuristic. 3  for  interchanging  the  boundary 
nodes  between  adjacent  clusters  to  improve  the  gains  obtained  by  their  first  phase  if 
possible.  The  heuristics  are  derived  from  the  fact  that  the  limited  gain  of  increase 
in  cost  due  to  the  constraint  on  the  cluster  size  in  their  first  phase  might  be  further 
improved  by  interchanging  the  boundary  nodes  without  violating  the  constraint  on 
the  cluster  size. 

5.3    Experimental  Testing  and  Analysis 

We  obtain  the  initial  partition  for  experimental  testing  of  the  heuristics  by  two 
different  methods:  random  and  centering.  The  random  partition  is  based  on  the 
topology  matrix  because  it  is  only  concerned  with  the  geographical  arrangement 
of  base  stations  and  the  cluster  size  constraint.  On  the  other  hand,  the  centering 
partition  is  based  on  both  the  topology  and  relative  cost  matrices  because  it  further 
considers  the  traffic  pattern  among  base  stations. 

The  centering  partition  of  m  clusters  is  achieved  by  a  two-phase  algorithm.  In 
the  first  phase,  the  m  center  nodes  on  which  traffics  are  concentrated  are  selected  by 
using  the  relative  cost  matrix.  Each  center  node  forms  an  initial  intermediate  cluster. 
For  every  intermediate  cluster,  the  second  phase  identifies  the  adjacent  nodes  of  the 
cluster  which  are  not  involved  in  other  clusters  by  using  the  topology  matrix  and 
chooses  one  of  them  by  using  the  relative  cost  matrix,  which  produces  a  maximum 
increase  in  cost  when  it  is  involved  in  the  cluster.  Then  the  maximum  node  is  added 
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Table  5.1.  Number  of  Times  in  100  Trials  a  Heuristic  Produced  a  Solution  with 
Maximal  Total  Cost  with  Respect  to  Other  Heuristics  When  m  =  3  and  1/2  x 
\n/m]  <  \Pi\  <  2  X  [n/m] 


Heuristics 

Random 

Centering 

Total 

n=19 

n=37 

n=61 

n=91 

n=19 

n=37 

n=61 

n=91 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

4 

17 

6 

2 

1 

17 

8 

1 

2 

54 

5 

84 

48 

47 

37 

88 

76 

56 

50 

486 

6 

84 

51 

43 

37 

86 

80 

70 

45 

496 

7 

85 

51 

47 

33 

85 

71 

70 

45 

487 

8 

70 

59 

39 

41 

76 

62 

57 

49 

453 

9 

70 

62 

41 

40 

76 

63 

60 

52 

464 

10 

70 

61 

43 

38 

75 

62 

61 

52 

462 

to  the  intermediate  cluster.  This  process  is  repeated  until  all  nodes  are  contained 
one  of  the  m  clusters. 

In  experimental  testing,  the  100  instances  of  the  relative  cost  matrix  were  ran- 
domly generated  for  each  of  several  values  of  n,  where  n  =  19, 37, 61,  and  91  for  H3, 
H4,  Hs,  and  He,  respectively.  For  each  value  of  n,  Heuristic.l  through  Heuristic. 10 
were  extensively  tested  on  the  100  cost  matrix  instances  for  each  of  the  four  cases 
which  are  the  combinations  of  the  two  methods  for  obtaining  the  initial  partition, 
random  or  centering,  and  the  number  of  clusters  in  =  3  or  4.  The  constraint  on  the 
cluster  size  is  given  by  1/2  x  [n/m]  <  |Pi|  <  2  x  p/t/m]. 

Table  5.1  and  Table  5.2  present  the  number  of  times  a  heuristic  produces  a  solution 
that  is  maximal  with  respect  to  the  solutions  produced  by  the  other  heuristics  when 
m  =  3  and  4,  respectively.  The  tables  show  that  Heuristic.5  through  Heuristic.  10 
which  are  the  combinations  of  the  techniques  of  interchanging  or  moving  boundary 
nodes  greatly  outperform  the  other  heuristics  using  only  one  of  the  techniques.  This 
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Table  5.2.  Number  of  Times  in  100  Trials  a  Heuristic  Produced  a  Solution  with 
Maximal  Total  Cost  with  Respect  to  Other  Heuristics  When  m  =  4  and  1/2  x 
[n/m]  <  I  P.  I  <  2  X  [n/m] 


Random 

Centering 

Heuristics 

n=19 

n=37 

n=61 

n=91 

n=19 

n=37 

n=61 

n=91 

Total 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

4 

8 

2 

2 

2 

7 

2 

2 

1 

26 

5 

43 

46 

41 

33 

73 

71 

58 

52 

417 

6 

43 

45 

43 

42 

73 

72 

61 

52 

431 

7 

64 

53 

56 

45 

84 

71 

63 

60 

496 

8 

68 

59 

34 

27 

80 

74 

66 

41 

449 

9 

68 

59 

36 

27 

80 

75 

71 

45 

461 

10 

68 

61 

37 

29 

80 

76 

66 

43 

460 

is  due  to  the  fact  that  although  the  cost  between  a  pair  of  adjacent  clusters  cannot 
be  improved  by  interchanging  their  boundary  nodes,  moving  their  boundary  nodes 
might  further  improve  the  cost,  which  in  turn  changes  the  boundary  between  the 
adjacent  clusters  for  possible  interchanging  at  the  next  step,  and  vice  versa.  The 
tables  also  show  that  as  the  number  of  nodes  n  increases,  the  probability  that  a 
heuristic  produces  a  maximal  solution  becomes  lower. 

It  is  interesting  that  Heuristic.  1  through  Heuristic. 3  do  not  produce  maximal 
solutions  at  all  even  though  Heuristic. 4  produces  a  small  number  of  solutions  being 
maximal.  This  is  because  interchanging  boundary  nodes  does  not  change  the  number 
of  nodes  for  each  cluster  in  the  initial  partition  and  so  less  flexible  than  moving 
boundary  nodes.  Thus,  unlike  general  graph  partitioning  problems,  the  algorithm 
which  is  only  based  on  one  of  the  techniques  of  interchanging  or  moving  boundary 
nodes  is  no  longer  useful  for  our  problem  which  additionally  considers  the  underlying 
topology. 
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Table  5.3.  Maximum  Difference  in  Total  Cost  From  Maximal  Solution  When  m  =  3 
and  1/2  x  [n/m]  <  |P,  |  <  2  x  fn/m] 


Heuristics 

Random 

Centering 

n=19 

n=37 

n=61 

n=91 

n=19 

n=37 

n=61 

n=91 

1 

44.2 

36.4 

34.6 

31.4 

42.8 

37.1 

35.0 

34.7 

2 

44.2 

36.4 

34.6 

31.2 

42.8 

37.1 

35.0 

34.7 

3 

44.2 

36.4 

33.1 

31.0 

42.8 

37.1 

35.0 

34.7 

4 

9.3 

6.3 

5.6 

4.3 

11.2 

17.7 

5.7 

23.4 

5 

8.7 

28.1 

25.9 

13.5 

27.3 

3.2 

23.5 

23.1 

6 

8.7 

28.1 

25.9 

18.9 

27.3 

3.2 

23.5 

23.1 

7 

8.7 

26.3 

25.9 

18.9 

27.3 

4.5 

2.4 

1.6 

8 

7.2 

23.3 

25.6 

23.9 

3.9 

3.7 

2.7 

2.0 

9 

7.2 

23.3 

25.6 

23.9 

3.9 

3.7 

2.7 

2.0 

10 

7.2 

23.3 

3.1 

23.9 

3.9 

3.7 

2.7 

2.0 

Table  5.4.  Maximum  Difference  in  Total  Cost  From  Maximal  Solution  When  m  —  i 
and  1/2  x  \n/m]  <  |P,|  <  2  x  \n/m] 


Heuristics 

Random 

Centering 

n=19 

n=37 

n=61 

n=91 

n=19 

n=37 

n=61 

n=91 

1 

40.2 

33.1 

29.5 

26.3 

40.6 

34.6 

31.5 

29.9 

2 

40.2 

33.1 

29.7 

26.1 

40.6 

34.6 

31.5 

29.9 

3 

40.2 

33.1 

29.7 

26.0 

40.6 

34.6 

31.5 

29.9 

4 

16.0 

8.2 

5.5 

5.0 

23.3 

17.9 

9.4 

5.4 

5 

11.2 

5.2 

10.3 

4.5 

14.0 

17.0 

12.8 

3.7 

6 

11.2 

5.2 

10.3 

3.5 

14.0 

17.0 

12.8 

3.7 

7 

9.8 

5.2 

2.7 

4.1 

6.8 

4.5 

3.3 

2.1 

8 

13.1 

5.4 

3.9 

4.3 

6.8 

4.7 

3.6 

3.9 

9 

13.1 

5.4 

3.9 

4.3 

6.8 

4.7 

3.6 

3.9 

10 

13.1 

5.4 

3.9 

4.3 

6.8 

4.7 

3.6 

3.9 
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Table  5.5.  Average  Difference  in  Total  Cost  From  Maximal  Solution  When  m  =  3 
and  1/2  x  |"n/m]  <  |P,|  <  2  x  [n/m] 


Heuristics 

Random 

Centering 

n=19 

n=37 

n=61 

n=91 

n=19 

n=37 

n=61 

n=91 

1 

38.2 

33.6 

31.0 

29.2 

35.1 

32.3 

31.7 

31.9 

2 

38.2 

33.6 

31.0 

29.2 

35.1 

32.2 

31.7 

31.9 

3 

38.2 

33.6 

31.0 

29.1 

35.1 

32.2 

31.7 

31.9 

4 

3.0 

2.3 

2.3 

1.7 

3.1 

2.4 

1.6 

1.5 

5 

0.4 

1.6 

1.0 

0.8 

0.7 

0.2 

0.7 

0.5 

6 

0.4 

1.3 

1.4 

0.8 

0.7 

0.2 

0.6 

0.5 

7 

0.4 

1.0 

1.3 

1.1 

0.6 

0.4 

0.2 

0.3 

8 

0.9 

0.7 

0.9 

0.9 

0.4 

0.4 

0.3 

0.3 

9 

0.9 

0.6 

0.9 

0.9 

0.4 

0.4 

0.3 

0.3 

10 

0.9 

0.6 

0.6 

0.7 

0.4 

0.4 

0.3 

0.3 

Table  5.6.  Average  Difference  in  Total  Cost  From  Maximal  Solution  When  m  —  A 
and  1/2  x  [n/m]  <  |P.|  <  2  x  [n/m] 


Heuristics 

Random 

Centering 

n=19 

n=37 

n=61 

n=91 

n=19 

n=37 

n=61 

n=91 

1 

31.3 

29.2 

27.5 

24.2 

31.8 

29.3 

26.1 

25.4 

2 

31.3 

29.2 

27.5 

24.2 

31.8 

29.3 

26.1 

25.5 

3 

31.3 

29.2 

27.4 

24.2 

31.8 

29.3 

26.1 

25.4 

4 

5.9 

3.3 

2.8 

2.1 

5.2 

3.2 

2.2 

2.0 

5 

2.3 

0.7 

0.7 

0.7 

1.1 

0.7 

0.5 

0.4 

6 

2.3 

0.6 

0.7 

0.5 

1.1 

0.7 

0.5 

0.4 

7 

0.9 

0.5 

0.4 

0.4 

0.5 

0.5 

0.3 

0.3 

8 

0.9 

0.5 

0.7 

0.6 

0.6 

0.5 

0.2 

0.5 

9 

0.9 

0.5 

0.8 

0.6 

0.6 

0.4 

0.2 

0.5 

10 

0.9 

0.5 

0.7 

0.6 

0.6 

0.4 

0.2 

0.5 
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Table  5.3  and  Table  5.4  present  the  percent  of  maximum  difference  in  cost  from  the 
solution  with  the  maximal  cost,  while  Table  5.5  and  Table  5.6  present  the  percent  of 
average  difference  in  cost  from  the  solution  with  the  maximal  cost.  The  tables  confirm 
the  superiority  of  Heuristic. 5  through  Heuristic.  10  with  respect  to  the  other  hueristics. 
In  general,  the  maximum  differences  for  the  centering  partition  of  Heuristic. 8  through 
Heuristic.  10  are  minimal  with  respect  to  the  random  partition  of  Heuristic.8  through 
Heuristic.  10  and  both  the  random  and  centering  partitions  of  Heuristic. 5  through 
Heuristic. 7  with  some  exceptional  test  cases.  The  average  differences  for  the  centering 
partition  of  Heuristic.8  through  Heuristic.  10  are  also  minimal  in  general. 

5.4    Chapter  Summarv 

In  cellular  systems,  a  hexagonal  mesh  of  n  base  stations  may  generate  multiple 
types  of  traffic  among  themselves.  We  have  considered  the  problem  of  finding  a  cover 
of  disjoint  clusters  of  base  stations,  so  as  to  minimize  the  total  communication  cost  for 
the  entire  system.  The  problem  differs  from  general  graph  partitioning  problems  in 
that  it  additionally  considers  the  underlying  physical  topology  among  base  stations. 
We  have  developed  several  heuristics  based  on  the  combinations  of  the  techniques 
of  moving  or  interchanging  the  boundary  nodes  between  adjacent  clusters.  These 
heuristics  produce  optimal  partitions  with  respect  to  the  initial  partition  obtained 
randomly  or  by  centering.  The  heuristics  are  compared  and  shown  to  behave  quite 
well  through  experimental  testing  and  analysis. 


CHAPTER  6 

A  PERFORMANCE  ANALYSIS  OF  THE  VIRTUAL  CELL  SYSTEM 

In  this  chapter,  we  deal  with  the  performance  analysis  of  the  virtual  cell  system. 
Once  an  optimal  partition  of  m  disjoint  clusters  is  obtained  by  the  algorithms  in  Chap- 
ter 4  and  Chapter  5,  the  virtual  cell  system  is  deployed  as  shown  in  Figure  3.3  where 
each  cluster  corresponds  to  a  virtual  cell.  Virtual  cell  i,  where  1  <  i  <  m,  is  imple- 
mented by  interconnecting  the  base  stations  of  the  ith  cluster  and  an  ARP/Location 
server  by  a  base  station  network.  These  in  virtual  cells  are  interconnected  by  the 
backbone  network.  The  virtual  cell  system  is  modeled  as  a  BCMP  open  multiple 
class  queueing  network.  Both  the  move  and  find  frequencies  among  base  stations 
and  the  topology  of  the  virtual  cell  system  are  used  to  determine  service  transition 
probabilities  in  the  queueing  network  model  as  well  as  the  arrival  rate  for  each  type 
of  messages.  There  are  three  types  of  messages  entering  or  leaving  the  virtual  cell 
system  via  base  stations:  the  handoff  message,  the  data  message,  and  the  address 
resolution  message.  The  handoff  messages  are  generated  due  to  move  operations  and 
the  data  and  address  resolution  messages  are  due  to  find  operations.  By  solving 
traffic  equations,  the  performance  measures  such  as  the  network  response  time  and 
protocol  processing  loads  at  network  components  are  obtained. 

6.1    Performance  Model 

We  adopt  a  BCMP  open  multiple  class  queueing  network  to  model  the  virtual 
cell  system,  as  depicted  in  Figure  6.1.  A  virtual  cell  is  modeled  as  a  number  of  base 
station  nodes,  an  ARP/Location  server  node,  and  a  base  station  network  node  which 
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Virtual  cell  i 

:  Base  staticm 

1  Base  station  j 

I  Base  station  : 

i  — KDO-  I 

:  Base  station  network  I 

I  I 

i  ARP/Location  server 

Backbone  network 

Figure  6.1.  The  Performance  Model  of  the  Virtual  Cell  System 

captures  traffic  characteristics  among  physical  cells  in  the  same  virtual  cell.  But,  in 
order  to  capture  traffic  characteristics  between  virtual  cells,  a  separate  service  node 
is  used  to  model  the  backbone  network.  Messages  from  mobile  hosts  enter  and  leave 
the  network  model,  only  going  through  base  station  nodes  in  the  virtual  cell  system. 

The  base  station  nodes  of  virtual  cell  i  are  sequentially  indexed  as  1,2, ...  ,n,-  in 
an  arbitrary  order,  where  n,  is  the  number  of  base  stations  in  virtual  cell  i  such  that 
Hill  '^i  =  Denote  ij  as  service  node  j  in  virtual  cell  i.  The  network  model  with 
N  =  n  +  2m  +  1  nodes  is  defined  as  follows: 

1.  A  base  station  node  is  labeled  ij,  where  1  <  z  <  m  and  1  <  j  <  n,-,  and 
modeled  as  an  FCFS  type  of  service  station  with  a  fixed  service  rate  for 
message  class  r.  The  messages  in  the  queue  are  transmitted  to  its  base  station 
network  node  or  leave  the  network  model. 
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(a)  The  labeling  scheme  in  the  optimization  algorithms 
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(b)  Tlie  corresponding  labeling  in  the  network  model 


Figure  6.2.  Labeling  the  Topology  of  the  Virtual  Cell  System 

2.  A  base  station  network  node  is  labeled  iB,  where  1  <  i  <  m,  and  modeled 
as  a  PS  type  of  service  station  with  a  fixed  service  rate  fiiBr  for  message  class 
r.  The  messages  in  the  queue  are  transmitted  to  its  base  station  nodes  or 
ARP/Location  server  node. 

3.  An  ARP/Location  server  node  is  labeled  iS,  where  1  <  i  <  m,  and  modeled 
as  an  FCFS  type  of  service  station  with  a  fixed  service  rate  ^isr  for  message 
class  r.  The  messages  in  the  queue  are  transmitted  to  its  base  station  network 
nodes  or  the  backbone  network  node. 

4.  The  backbone  network  node  is  labeled  Bl  and  modeled  as  a  PS  type  of  service 
station  with  a  fixed  service  rate  fisir  for  message  class  r.  The  messages  in  the 
queue  are  transmitted  to  ARP/Location  server  nodes  of  virtual  cells. 

The  analysis  is  based  on  the  following  assumptions: 

1.  Given  a  hexagonal  mesh  of  base  stations  of  size  5,  it  is  known  that  the  number 
of  base  stations  are  n  =  3s  —  3^  -|-  1  and  the  number  of  columns  in  each 
three  directions  is  d  =  2s  -  I.  From  left  to  right,  each  column  is  indexed  as 
0,1,2,..., J— 1.  Let  i  be  the  column  index.  Then  the  nodes  of  every  column  i 
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are  sequentially  labeled  i  x  d,i  x  d+l,i  x  d  +  2, . . from  bottom  to  top.  Denote 
nindex  and  cindex  as  the  node  and  cluster  indices,  respectively.  An  optimal 
partition  produced  by  the  algorithms  [23,  24]  can  be  represented  by  an  index 
pair  [nindex,  cindex)  for  every  base  station.  Figure  6.2(a)  depicts  an  optimal 
partition  for  a  hexagonal  mesh  of  size  2.  Note  that  the  node  index  reflects 
the  underlying  topology  of  the  hexagonal  mesh  of  base  stations.  Thus,  given 
the  node  index  of  a  base  station,  its  neighboring  base  stations  can  be  directly 
identified  by  the  labeling  scheme  of  the  algorithms. 

In  addition  to  the  node  index,  the  network  model  of  the  virtual  cell  system  also 
requires  base  stations  to  be  logically  indexed  because  a  base  station  in  a  virtual 
cell  can  change  into  another  virtual  cell  according  to  traffic  patterns.  Denote 
vindex  and  bindex  as  the  virtual  cell  and  base  station  indices,  respectively. 
The  virtual  cell  index  is  directly  mapped  to  the  cluster  index  and  the  base 
station  index  is  logically  assigned  so  that  the  base  stations  in  a  virtual  cell  is 
sequentially  labeled  1,2,...,  in  an  arbitrary  order.  Figure  6.2(b)  depicts  the 
labeling  of  the  network  model  which  corresponds  to  an  optimal  partition  in 
Figure  6.2(a).  Thus,  given  base  station  xy  in  the  network  model,  where  x  and 
y  are  the  virtual  cell  and  base  station  indices,  respectively,  its  neighboring  base 
stations  are  directly  identified  by  the  node  index.  Let  Adj{xy)  be  a  set  of  base 
stations  adjacent  to  xy.  From  Adj{xy)  we  can  identify  which  neighboring  base 
station  belongs  to  which  virtual  cell. 

2.  We  use  a  flow-based  mobility  model  [25]  which  assumes  that  mobile  hosts  are 
uniformly  distributed  in  the  area  of  a  physical  cell  and  the  travel  direction  of 
each  mobile  host  with  respect  to  the  border  is  uniformly  distributed.  If  we 
define  p:cy  to  be  the  density  of  the  mobile  hosts  per  km^  at  the  area  of  base 


91 


station  xy,  v^y  to  be  the  average  speed  in  km/ sec  of  a  mobile  host  at  the  area  of 
base  station  xy,  and  L  to  be  the  length  of  the  perimeter  of  the  area  of  a  physical 
cell,  then  the  average  number  of  mobile  hosts  per  sec  leaving  base  station  xy 
is  given  by  Mxy  =  PxyVxyL/r.  Denote  fm{xy,  x'y')  as  the  move  frequency  from 
base  station  xy  to  one  of  its  neighboring  base  station  x'y'  G  Adj{xy).  Then 
fm{xy,x'y')  =  Mxy/\Adj{xy)\,  where  \Adj{xy)\  is  the  number  of  base  stations 
adjacent  to  xy  and  set  to  6  in  the  hexagonal  arrangement  of  n  base  stations. 

3.  If  we  define  b  to  be  the  data  rate  in  bits /sec  of  a  wireless  channel,  /  to  be  the 
average  message  length  in  bits,  and  E  to  be  the  in-call  probability  in  Erlangs 
for  a  mobile  host,  then  the  non-blocking  data  arrival  rate  in  messages/ sec 
at  base  station  xy  is  given  by  Fx,j  =  ttB^ p^yEb/l,  where  R  is  the  radius  in 
km  of  a  physical  cell.  Denote  fj[xy,x'y')  as  the  find  frequency  from  base 
station  xy  to  base  station  x'y' .  Then  ff{xy,x'y')  is  a  fraction  of  Fxy,  such  that 
Er'y'  fj{xy,  x'y')  =  Fxy. 

4.  The  address  resolution  frequency  is  derived  from  the  find  frequency.  If  we 
define  Fxy  to  be  the  arrival  rate  in  messages/ sec  of  data  messages  at  base 
station  xy,  which  are  destined  to  base  stations  in  the  same  virtual  cell  x,  c  to 
be  the  average  number  of  messages  over  a  conversation  from  the  source  mobile 
host  to  the  destination  mobile  host,  and  h,,,  to  be  the  miss  ratio  of  the  address 
resolution  table  of  a  mobile  host,  then  the  arrival  rate  in  messages /sec  of  the 
address  resolution  messages  at  base  station  xy  is  given  by  Axy  =  F^yhm/c,  where 

^xy  —  ITxy' 

ff{xy,xy').  Denote  fa{xy,xy')  as  the  address  resolution  frequency 
from  base  station  xy  to  base  station  xy'.  Then  fa{xy,xy')  is  a  fraction  of  Axy, 
such  that  Zxy'  fa{xy,xy')  =  A^y. 
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It  should  be  noted  that  the  move  frequency  not  only  affects  the  arrival  rate  of  the 
handoff  messages  but  also  the  forwarding  rate  of  the  data  and  address  resolution  reply 
messages.  For  example,  consider  the  move  frequency  fm{xy,  xy').  When  a  mobile  host 
at  base  station  xy  moves  into  bcise  station  xy\  it  initiates  a  handoff  request  message 
to  xy' .  But  a  data  message  destined  to  the  mobile  host  may  be  forwarded  to  base 
station  xy'  at  base  station  xy  during  the  handoff  process.  In  the  same  way,  the 
address  resolution  reply  message  for  the  mobile  host  from  the  ARP/Location  server 
should  be  also  forwarded  to  base  station  xy'  at  base  station  xy. 

6.1.1    Multiple  Class  Traffic  Model 

There  are  three  classes  of  messages  entering  the  network:  the  handoff  request 
message,  the  data  message,  and  the  address  resolution  request  message.  As  each 
class  of  message  traverses  through  the  network,  it  not  only  requires  different  service 
requirements  and  different  routing  behavior,  but  changes  its  class.  For  example,  the 
base  station  received  a  handoff  request  message  from  a  mobile  host  sends  the  message 
to  the  previous  base  station  of  the  mobile  host.  Then  the  previous  base  station  mul- 
ticasts  the  message  in  its  virtual  cell  to  update  the  distributed  location  information 
of  the  mobile  host,  and  at  the  same  time  it  sends  a  handoff  response  message  to  the 
new  base  station  from  which  the  mobile  host  receives  a  handoff  confirmation  message. 
Thus,  as  a  message  of  the  handoff  request  class  progresses  through  the  network,  it  is 
changed  into  a  message  of  the  multicast  class  and  next  into  a  message  of  the  handoff 
response  class. 

After  the  handoff  completes,  a  data  message  destined  to  the  mobile  host  will  be 
directly  deUvered  to  the  new  base  station.  However,  during  the  handoff  the  data 
message  will  be  first  delivered  to  the  previous  base  station,  which  in  turn  forwards  it 
to  the  new  base  station.  Thus,  as  a  message  of  the  data  class  progresses  through  the 
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network,  it  might  be  changed  into  a  forwarding  message  which  may  be  involved  in 
a  virtual  cell  or  between  virtual  cells.  The  forwarding  message  involved  in  a  virtual 
cell  is  referred  to  as  the  intra-forwarding  class  message,  while  that  involved  between 
virtual  cells  is  referred  to  as  the  inter-forwarding  class  message.  In  the  same  way, 
as  a  message  of  the  address  resolution  request  class  progresses  through  the  network, 
it  is  changed  into  a  message  of  the  address  resolution  reply  class  at  ARP/Location 
servers  and  possibly  the  address  resolution  reply  forwarding  class  at  beise  stations. 

Hence,  the  set  of  message  classes  in  the  network  can  be  partitioned  into  three 
subsets  which  contain  the  message  classes  related  to  handoff,  data,  and  address  res- 
olution, respectively.  Denote  z  =  1,2,  and  3  as  the  subsets  of  the  message  classes 
related  to  handoff,  data,  and  address  resolution,  respectively.  Even  though  any  class 
r  in  a  subset  z  can  visit  the  entire  set  of  nodes  N,  we  cannot  define  a  routing  chain 
for  each  subset  z  because  a  class  r  message  in  a  subset  z  may  require  different  routing 
behavior  due  to  the  topology  of  the  virtual  cell  system. 

For  example,  consider  a  handoff  request  message  from  a  mobile  host  which  moves 
from  an  adjacent  base  station  into  base  station  12  in  the  example  of  Figure  6.2(b).  If 
the  mobile  host  is  from  base  station  11,  the  message  will  be  directly  delivered  to  its 
previous  base  station  11  via  base  station  network  IB.  However,  if  the  mobile  host  is 
from  base  station  32,  the  message  will  be  delivered  to  its  previous  base  station  32, 
going  through  base  station  network  IB,  ARP/Location  server  IS,  backbone  network 
Bl,  ARP/Location  server  36*,  and  finally  base  station  network  SB .  Thus,  in  addition 
to  the  message  class,  the  topology  of  the  virtual  cell  system  should  be  considered  to 
determine  the  routing  behavior  of  messages  in  the  network. 

In  order  to  incorporate  the  topology  of  the  virtual  cell  system  into  the  queueing 
network  model,  it  is  necessary  to  identify  which  base  station  generates  the  message 
classes  for  each  subset  z.  Thus,  we  define  a  routing  chain  for  the  message  classes  of 
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each  subset  z  generated  by  each  base  station.  Then  there  are  3n  routing  chains  in 
the  network,  denoted  cis  where  1  <  x  <  m,  1  <  y  <  Ux,  and  z  —  1,2,  and  3. 

For  each  routing  chain  Ex,y,z,  the  service  transition  probabihties  are  defined  by  the 
set  {pijr;kis},  which  describes  the  probability  that  a  class  r  message  at  node  ij  goes 
next  to  node  kl  as  a,  class  s  message,  where  r,  s  ^  z  and  ij,  kl  G  N. 

6.1.2    Arrival  Process 

We  assume  a  Poisson  state-independent  arrival  stream  for  each  routing  chain. 
Denote  \x,y,z  as  the  arrival  rate  corresponding  to  a  routing  chain  Ex,y,z-  Then  the 
arrival  rate  \x,y,z  is  determined  by  the  move,  find,  and  address  frequencies  among  n 
base  stations. 

1.  A  mobile  host  migrating  from  base  station  x'y'  E  Adj{xy)  to  base  station  xy 
initiates  a  handoff  request  message  to  base  station  xy.  Thus,  the  arrival  rate 
for  the  handoff  request  class  of  messages  at  base  station  xy  can  be  represented 

by  Xx,y,l  =  Ilx'y'&Adj(xy)  fm{x'y',Xy). 

2.  A  data  message  transmitted  by  a  mobile  host  at  base  station  xy  needs  a  find 
operation  to  locate  base  station  x'y'  to  which  the  destination  mobile  host  be- 
longs. Thus,  the  arrival  rate  for  the  data  class  of  messages  at  base  station  xy 
can  be  represented  by  Xx,y,2  =  Ex',/  ff{xy,x'y'). 

3.  The  source  mobile  host  performs  address  resolution  operations  only  when  it  re- 
sides in  its  native  virtual  cell  and  the  network  address  of  the  destination  mobile 
host  indicates  the  same  native  virtual  cell.  Thus,  the  arrival  rate  for  the  address 
resolution  request  class  of  messages  at  base  station  xy  can  be  represented  by 

^x,y,3  =  T.x'y',x<=x  fa{xy,x'y'). 
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6.1.3    Service  Transition  Matrix 

A  service  transition  matrix  Px,y,z  =  [pi]T;kis]  is  to  be  defined  for  each  routing  chain 
Ex,y,z-  To  determine  service  transition  probabilities,  it  is  necessary  to  consider  not 
only  the  move,  find,  and  address  resolution  frequencies,  but  also  the  topology  of  the 
virtual  cell  system.  For  the  handoff  message  classes,  base  station  xy  may  send  and 
receive  a  handoff  request  class  message  and  a  handoff  response  class  message  to  and 
from  an  adjacent  base  stations  st  6  Adj{xy)  on  a  different  route  in  the  network, 
depending  on  whether  base  station  st  is  located  in  the  same  virtual  cell  or  a  different 
virtual  cell.  For  the  data  message  classes,  in  the  same  way,  base  station  xy  may  deliver 
a  data  clciss  message  to  bcise  station  st  G  Adj(xy)  on  a  different  route,  depending  on 
whether  base  station  st  is  located  in  the  same  virtual  cell  or  a  different  virtual  cell. 
If  the  data  class  message  arriving  at  base  station  st  is  to  be  forwarded  to  an  adjacent 
base  station  s't'  G  Adj{st),  the  forwarding  message  may  also  take  a  different  route  in 
the  network  model  depending  on  whether  s't'  and  st  are  in  the  same  virtual  cell  or 
different  virtual  cells. 

Handoff  message  classes.  Consider  the  handoff  request,  multicast,  and  handoff 
response  message  clcisses  originated  from  base  station  xy,  i.e,  a  routing  chain  Ex,y,i. 
Define  S^t  =  fm{st,xy)/Xx,y,i  to  be  a  probability  that  a  mobile  host  at  base  station 
st  €  Adj(xy)  moves  into  base  station  xy.  Denote  xy'  as  a  base  station  adjacent  to  xy 
in  the  same  virtual  cell  x  and  x'y'  as  a  base  station  adjacent  to  xy  in  a  different  virtual 
cell  x' .  Then  we  can  compute  the  following  probabilities  from  the  move  frequency: 

•  The  probability  that  a  mobile  host  at  base  station  xy'  moves  into  base  station 
xy  is  given  by  S^y-  =  fm{xy',xy)/\^^y^i. 

•  The  probabiHty  that  a  mobile  host  in  virtual  cell  x  moves  into  base  station  xy 

is  given  by        =  Exy'eAdJ^xy)  ^xy'- 
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•  The  probability  that  a  mobile  host  at  base  station  x'y'  moves  into  bcise  station 
xy  is  given  by        =  fm{x'y' ,xy)l\^^y^i. 

•  The  probability  that  a  mobile  host  in  a  different  virtual  cell  x'  moves  into  base 
station  xy  is  given  by  S^'  =  T,x'y'€Adj{xy)  ^x'y'- 

The  service  transition  matrix  Px,y,\  foi"  the  handoff  message  classes  originated 
from  base  station  xy  is  given  as  follows: 

•  For  the  handoff  message  classes  when  a  mobile  host  moves  from  base  station 
xy'  €  Adj{xy)  into  base  station  xy,  service  transition  probabilities  are  as  fol- 
lows: 

—  Pxyi;xB\  =  1-0  and  PxB\\xy'\  —  ^xy'  fof  the  handoff  request  cleiss, 

"~  Pxy'\;xB2  —  1-0  ^.nd  PxB2;xy'2  —  ^xy'/^x  for  the  multicast  class,  and 

—  Pxy'2;xB3  =  1-0  Bud  PxB3;xy3  =  1-0  for  tlic  haudoff  responsc  class. 

•  For  the  handoff  message  classes  when  a  mobile  host  moves  from  base  station 
x'y'  €  Adj(xy)  into  base  station  xy,  service  transition  probabilities  are  as  fol- 
lows: 

—  Pxyl;xBl  =  l-0,PiBl;xSl  =  PxSl;Bll  =  1-0-,  PB11;x'S1  =  ^x' /  T,x'eAdj(xy)  ^x' , 

Px'Si;x'Bi  =  1-0,  and  Px>Bi;x'y'i  =  ^x'y'/K'  for  the  handoff  request  class, 

—  Px'v'i;x'B2  =  1-0  and  Px'B2;x'y'2  —  ^x'y'j^x'  for  the  multicast  class,  and 

—  Px'y'2;x'B3  =  1-0,  Px'B3;x'S3  =  1-0,  Px'53;Bl3  =  1-0,  PB13;xS3  =  1-0,  PxS3;xB3  = 

1.0,  and  PxB3;xy3  =  1.0  for  the  handoff  response  class. 

Data  message  classes.  Consider  the  data,  intra-forwarding,  and  inter-forwarding 
message  classes  originated  from  base  station  xy,  i.e,  a  routing  chain  £^x,y,2-  Define 
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=  ff{xy,st)/Xx,y,2  to  be  a  probability  that  a  mobile  host  at  base  station  xy  sends 
a  data  message  to  base  station  st.  Denote  xy'  as  a  base  station  in  the  same  virtual 
cell  X  and  x'y'  as  a  base  station  in  a  different  virtual  cell  x'.  Then  we  can  calculate 
the  following  probabilities  for  the  data  message  class  from  the  find  frequency: 

•  The  probability  that  a  data  message  from  base  station  xy  is  destined  to  base 
station  xy'  is  given  by  a^j^'  =  ff{xy,xy')/X.,:^y'i. 

•  The  probability  that  a  data  message  from  base  station  xy  is  destined  to  the 
base  stations  in  virtual  cell  x  is  given  by  a.,,.  =  Ylxy'  ctxy'- 

•  The  probability  that  a  data  message  from  base  station  xy  is  destined  to  a  base 
station  x'y'  is  given  by  a^'y'  =  ff{xy,  x'y')/ X.,,^y;2. 

•  The  probability  that  a  data  message  from  base  station  xy  is  destined  to  the 
base  stations  in  a  different  virtual  cell  x'  is  given  by  a^'  =  J2x'y'  ctx'y'- 

Let  us  now  consider  forwarding  message  classes  in  the  network.  Consider  the 
case  that  the  source  mobile  host  at  base  station  xy  trys  to  send  a  data  message  to 
the  destination  mobile  host  which  is  moving  from  base  station  st  into  base  station 
s't'  G  Adj{st).  If  the  message  transmission  occurs  after  the  handoff  completes,  the 
message  should  be  directly  sent  to  base  station  s't'.  If  the  message  transmission  occurs 
before  the  handoff  completes,  however,  the  message  will  be  sent  to  the  previous  base 
station  st  where  it  will  be  forwarded  to  the  new  base  station  s't'.  Define  T^t  to  be  the 
average  time  in  second  for  a  mobile  host  to  stay  at  base  station  5^  before  it  leaves. 
Then  =  H^t/Mst,  where  H^t  is  the  total  number  of  mobile  hosts  at  base  station  st 
and  Mat  is  the  average  rate  of  mobile  hosts  per  second  moving  out  of  base  station  st. 
Given  the  average  handoff  time  r/i,  the  probability  that  a  message  destined  to  base 
station  st  is  forwarded  to  base  station  s't'  e  Adj{.<it)  is  given  by  a^t  =  Th/Tat\Adj{st)\. 
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This  should  be  a  valid  estimate  if  an  MH  migrates  independent  of  when  it  receives 
data  messages. 

Define  q°gt\ni  to  be  the  forwarding  frequency  going  out  of  base  station  st  to  base 
station  s't'  G  Adj{st),  to  be  the  forwarding  frequency  going  out  of  base  station  st 
to  all  neighboring  base  stations  in  Adj{st),  and  f;°"*  to  be  the  forwarding  frequency 
going  out  of  all  base  stations  in  virtual  cell  s  to  different  virtual  cells.  Then  q^^gui  = 

ff{xy,st)a-,t,  q^t*  =  'Es't'eAdjist)  Qtt,l't'^  ^^'^        =  T^st  T,s't'eAdj(st),s'ji,  ^°styt'- 

Define  q'J},it'  forwarding  frequency  going  out  of  base  station  s't'  G  Adj{st) 

into  base  station  st,  q*^  to  be  the  forwarding  frequency  going  out  of  all  neighboring 
base  stations  in  the  same  virtual  cell  s  into  base  station  st,  q^  to  be  the  forwarding 
frequency  going  out  of  all  neighboring  base  stations  in  different  virtual  cells  into 
base  station  st,  and  q*J^  to  be  the  forwarding  frequency  going  out  of  all  neighboring 
base  stations  in  different  virtual  cells  into  all  base  stations  in  virtual  cell  5.  Then 
9lt,s't'  =  ff{^y,s't')<Ttiti,  qll  =  ^s't'eAdj{st),s'=s  qlt,s't'^  '/if  =  T!,s't>eAdj{st),a':^s  qlt,s't'i  and 

From  the  above  forwarding  frequencies,  we  can  compute  the  following  probabilities 
used  to  determine  service  transition  probabilities  for  the  forwarding  message  classes: 

•  The  probability  that  a  forwarded  message  at  base  station  st  is  routed  to  a 
different  virtual  cell  is  given  by  =  qT^  T^st  qT- 

•  The  probability  that  a  forwarded  message  at  base  station  st  is  routed  to  base 
station  st'  in  the  same  virtual  cell  5  is  given  by  /Ji"/'""  =  q'Jli/JZst'  9lt'- 

•  The  probability  that  a  forwarded  message  at  base  station  st  is  routed  to  the 
base  stations  in  a  different  virtual  cell  s'  is  given  by  7^7*^'"  =  qiVf^gi  q^V. 


99 


•  The  probability  that  a  forwarded  message  at  base  station  st  is  routed  to  base 
station  s't'  is  given  by  7,/t.  =  q^/  Zs't'  Q'^- 

The  service  transition  matrix  Px,y,2  for  the  data  message  class  generated  at  beise 
station  xy  is  given  as  follows: 

•  When  a  data  class  message  is  destined  to  base  station  xy'  in  the  same  virtual 
cell  X,  service  transition  probabilities  are  as  follows: 

~  Pxyi;xBi  =  1-0  and  PxBi;xy'i  =  Oi^y' / otx  iov  the  data  class, 

-  Pxy>l;xB2   =   (T^y'\Adj{xy')\,  PxB2;xy'2    =    (1  "  ^^*")P'xy>^'' ,   and  PxB2;xS2  = 

^mter  £qj,  ^y^^  intra-forwarding  class,  and 

-  PxS2;B13  =  1-0,  PB13;x'S3  =  7i"'^%  Px'53;x'B3  =  1-0,  and  Px'B3;xy3  =  Ix'y'  for 

the  inter-forwarding  class. 

•  When  a  data  class  message  is  destined  to  base  station  x'y'  in  a  different  virtual 
cell  x',  service  transition  probabilities  are  as  follows: 

-  Pxyl;xBl  =  1-0,  PxBl;x51  =  I  -  a^,  PxSl;BU  =  1-0,  PBU;x'S1  —  Olx> / {I  -  OCx), 

Px'si;x'Bi  =  1-0,  and  Px'Bi;xyi  =  oi^iyi /a^-i  for  the  data  class, 

-  Px'y>V,x'B2  =  (Tx'y'\Adj{x'y')\,  Px>B2;x'y>2  =  (1  " /S^^'^O/^xv"'  Px'B2;x'S2  = 

^*?**'"  for  the  intra-forwarding  class,  and 

-  Px'S2;BlZ  =  1-0,  PBl3;x'53  =  Ti"'""^ )  Px'53;x'B3  =  1-0,  and  Px'B3;x'y'3  =  Ix'y'  for 

the  inter-forwarding  class. 

Address  resolution  message  classes.  Consider  the  address  resolution  request,  ad- 
dress resolution  reply,  and  address  resolution  reply  forwarding  message  classes  origi- 
nated from  base  station  xy,  i.e.,  a  routing  chain  E^^y^z-  Denote  Kbs  as  the  hit  ratio 
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of  the  cache  table  at  a  base  station  for  a  mobile  host.  The  service  transition  matrix 
■fx,y,3  for  address  resolution  message  classes  originated  from  base  station  xy  is  given 
as  follows: 

—  Pxy\;xBi  =  1  "  ^Bs  and  PxBi;xSi  =  1-0  for  the  address  resolution  request 
class, 

~  PxS\\xB2  =  1-0  and  PxB2;xy2  =  1-0  for  the  address  resolution  reply  class, 
and 

~  Pxy2;xB3  —  <^xy 

\Adj{xy)\  and  PxB3;xy'3  —  1  / 1 (^J/)  I  f^r  thc  addrcss  resolu- 
tion reply  forwarding  class. 

6.1.4    Traffic  Equations 

Suppose  that  e.jr  is  the  average  throughput  of  class  r  messages  through  node  ij 
in  a  routing  chain  Ex,y,z-  Then  for  node  ij  and  class  r,  {ij,r)  €  Ex,y,z,  we  can  write 
the  following  traffic  equations: 

^  ]         ^klsPkls;ijr  "i"  PO;ijr  —  (^ijrt 

where  poiov  is  the  probability  that  an  external  arrival  is  for  node  ij  and  class  r  and 
determined  by  the  rate  of  external  arrivals  of  class  r  messages  to  node  ij.  Since 
Po;ijr  >  0  for  some  {ij,r)  e  Ex,y,z,  Ex,y,z  is  open  and  so  the  traffic  equations  have  a 
unique  solution  for  {e,j>}. 

6.2    Performance  Measures 

The  BCMP  theorem  states  that  multiple  class  queueing  networks  with  the  FCFS, 
PS,  IS,  and  LCFS-PR  types  of  nodes  have  a  product  form  solution  for  the  steady 
state  joint  probability  distribution  of  the  node  states  [26].  Since  in  an  open  network 
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a  message  sees  the  network  in  the  steady  state  when  it  leaves  or  enters  the  network, 
Little's  result  can  be  applied  to  any  given  class  of  messages  at  a  node. 

Node  ij  has  a  fixed  service  rate  of  jj,  J  J  y.  3,11  d  relative  throughput  of  e,jv  for  class 
r  messages.  The  service  time  of  FCFS  nodes  is  independent  of  the  class  and  so 
1/Hijr  =  ^/fJ'ij  for  all  classes  r.  The  relative  throughput  e,jr  is  the  average  number 
of  times  that  a  class  r  message  visits  node  ij  before  leaving  the  network.  Then  the 
total  service  demand  of  a  class  r  message  at  node  ij  is  Dijr  =  Cijr/fiijT- 

Let  rjr  be  the  external  arrival  rate  of  each  class  /•.  The  performance  measures  we 
can  obtain  from  the  described  queueing  network  model  include: 

•  The  throughput  of  class  r  messages  in  the  steady  state  is  Tr  —  T]r- 

•  The  utilization  of  node  ij  by  class  r  messages  is  f/,jv  =  rjrDijr  by  Little's  result. 
Thus,  the  utilization  of  node  ij  is  given  by  =  J2reR''lrDijr,  where  R  is  the 
set  of  all  message  classes. 

•  The  mean  waiting  time  of  a  class  r  message  at  node  ij  for  each  visit  is  given 
by  Wi,,  =  (l//io>)/(l  -  Ui,). 

•  The  mean  time  that  a  class  r  message  spends  at  node  ij  during  its  stay  in 
the  network,  i.e.,  the  mean  residence  time  at  node  ij  for  class  r  is  given  by 

•  The  mean  time  that  a  class  r  message  spends  in  the  network,  i.e.,  the  system 
response  time  for  class  r  is  given  by  Q,.  =  YlijeN  Qijr- 

•  The  mean  number  of  class  r  messages  at  node  ij  is  Lijr  =  rjrQijr  by  Little's 
result. 
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6.3    Performance  Evaluations 

This  section  gives  the  results  of  computations  from  the  performance  measures  of 
the  queueing  network  model.  By  changing  the  density,  velocity  and  in-call  probability 
of  mobile  hosts  among  mobility  and  traffic  parameters,  we  evaluate  the  utilization  of 
various  network  components  and  the  system  response  time  due  to  various  messages. 
With  the  same  model  parameters  we  also  compare  two  different  virtual  cell  systems: 
one  deployed  according  to  an  initial  partition  and  the  other  deployed  according  to  an 
optimal  partition  with  respect  to  the  initial  partition.  The  initial  model  parameters 
are  listed  as  follows: 

•  System  parameters 


the  number  of  base  stations(n)  19 

the  number  of  virtual  cells(m)  3 

base  station  network  service  rate(/ijB)  45  Mbps 

backbone  network  service  rate(/iBi)  45  Mbps 

base  station  service  rate(/i,j)  4.5  Mbps 

ARP/Location  server  service  rate(/f,5)  •        4.5  Mbps 

•  Mobility  and  traffic  parameters 

average  density  of  mobile  hosts  at  physical  cell  xy{p^y)  500  mh/km^ 
average  velocity  of  mobile  hosts  at  physical  cell  xy{v^y)     15  km/hr 

cell  radius(i?)  1  km 

data  rate  of  a  wireless  channel(6)  2  Kbps 

average  message  length  (/)  4  Kbits 

Erlangs  per  mobile  host(^)  0.04 

average  number  of  messages  over  a  conversation(c)  240 

address  resolution  miss  ratio  at  a  mobile  host(/i„)  0.5 

address  resolution  hit  ratio  at  a  base  station( /1B5)  0.5 

probability  that  a  mobile  host  is  powered  on  0.5 

average  handoff  time(T;i)  1  sec 


The  analysis  assumes  that  the  initial  density  of  mobile  hosts  in  each  of  n  physical 
cells  is  randomly  varied  between  325  mh/km^  and  650  mh/km^  such  that  the  average 
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Figure  6.3.  The  Utilization  of  the  Network  Components  in  the  Virtual  Cell  System 
When  an  Initial  Partition  Is  Used 

density  becomes  500  mh/km^.  The  initial  velocity  of  mobile  hosts  in  each  of  n 
physical  cells  is  also  randomly  varied  between  5  km./hr  and  25  km/hr  such  that 
the  average  velocity  becomes  15  km/hr.  The  initial  density  of  each  physical  cell  is 
increased  with  the  same  ratio  at  each  step  such  that  the  average  density  is  raised 
from  500  mh/km'^  to  2000  mh/km'^.  After  the  average  density  increased  up  to  2000 
mh/km^,  the  average  velocity  is  then  increased  from  15  km/hr  up  to  30  km/hr  in 
the  same  way  and  then  the  in-call  probability  is  raised  from  0.5  up  to  1.0.  The  base 
station  network  service  rate  ms  and  the  backbone  network  service  rate  fiBi  represent 
the  service  rate  of  the  same  physical  transport  network. 

Figure  6.3  and  Figure  6.4  depict  the  utilization  of  the  network  components  in 
the  virtual  cell  system  which  is  deployed  according  to  an  initial  partition  and  an 
optimal  partition,  respectively.  Since  the  same  physical  transport  network  is  used 
for  both  base  station  networks  and  the  backbone  network  to  construct  a  virtual  cell 
system,  the  network  utilization  of  Figure  6.3  and  Figure  6.4  represents  the  sum  of  all 
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Figure  6.4.  The  Utilization  of  the  Network  Components  in  the  Virtual  Cell  System 
When  an  Optimal  Partition  Is  Used 

utilizations  for  three  base  station  networks  and  one  backbone  network.  The  utilization 
of  the  ARP/Location  server  and  base  stations  represents  the  average  utilization  of 
an  ARP/Location  server  or  a  base  station. 

As  depicted  in  Figure  6.3,  the  high  utilization  of  the  ARP/Location  server  for 
the  initial  partition  implies  that  there  is  a  large  volume  of  mobility  and  data  traffic 
between  clusters  before  the  optimization  process.  It  is  shown  in  Figure  6.4  that  the 
inter-cluster  traffic  is  significantly  reduced  after  the  optimization  process.  The  less 
inter-cluster  traffic  means  that  the  utilizations  of  the  backbone  network  and  base 
station  networks  are  also  reduced.  This  is  due  to  the  fact  that  for  handoff  or  data 
transfer  operations,  intra-cluster  traffic  involves  operations  at  only  one  base  station 
network  while  inter-cluster  traffic  involves  operations  at  two  base  station  networks 
and  one  backbone  network. 
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It  is  interesting  to  observe  that  the  average  utilization  of  a  beise  station  for  the  op- 
timal partition  is  slightly  higher  than  that  for  the  initial  partition  and  it  is  especially 
sensitive  to  the  velocity  of  mobile  hosts.  This  is  due  to  the  fact  that  the  optimiza- 
tion process  produces  as  large  clusters  as  possible  within  the  cluster  size  constraint 
in  order  to  reduce  the  total  communication  cost  for  the  entire  system.  Consider  the 
impact  of  handoff  multicast  operations  in  terms  of  cluster  size.  The  number  of  mobile 
hosts  leaving  a  base  station  is  determined  by  the  density  and  velocity  of  mobile  hosts 
in  its  physical  cell,  not  the  cluster  size.  For  each  leaving  mobile  host,  the  base  station 
invokes  a  multicast  operation  to  all  other  base  stations  in  the  same  cluster  in  order 
to  maintain  the  consistency  of  the  distributed  location  information  of  the  mobile 
host.  Thus,  a  handoff  multicast  operation  within  a  large  cluster  involves  more  base 
stations  than  that  within  a  small  cluster.  This  increases  the  total  utilization  of  base 
stations.  When  considering  only  the  handoff  multicast  operation,  the  best  partition 
would  consist  of  equal-weighted  clusters  where  the  weight  of  a  cluster  is  the  product 
of  the  number  of  base  stations  in  the  cluster  and  the  number  of  mobile  hosts  moving 
out  of  physical  cells  in  the  cluster.  As  the  velocity  of  mobile  hosts  increases,  more 
handoff  multicast  operations  will  be  needed,  which  in  turn  results  in  more  utilization 
of  base  stations. 

The  utilization  of  the  ARP/Location  server  for  the  optimal  partition  is  less  likely 
saturated  compared  to  that  for  the  initial  partition.  However,  the  average  utiHzation 
of  the  ARP/Location  server  for  the  optimal  partition  approaches  to  approximately 
0.72  at  the  saturation  point  of  the  in-call  probability  0.8.  This  means  that  the 
ARP/Location  server  for  the  largest  cluster  of  the  optimal  partition  is  saturated  at 
the  in-call  probability  0.8.  The  bottleneck  in  the  ARP/Location  server  comes  from 
the  fact  that  the  destination  of  data  messages  are  randomly  distributed  over  the 
entire  system  and  the  locality  of  data  traffic  patterns  is  not  considered.  Thus  the 
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Figure  6.5.  The  Utilization  of  the  ARP/Location  Server  for  Each  Type  of  Message 
When  an  Initial  Partition  Is  Used 

ARP/Location  server  may  use  a  faster  processor  to  resolve  it.  Adjusting  the  cluster 
size  constraint  in  the  optimization  process  could  also  alleviate  the  effect  of  the  largest 
cluster  on  the  utilization  of  the  ARP/Location  server.  In  this  analysis,  the  cluster 
sizes  of  the  initial  partition  are  ni  =  6,  ^2  =  6,  and  =  7,  and  the  cluster  sizes  of 
the  optimal  partition  obtained  by  using  the  cluster  size  constraint,  4  <  <  12,  are 
Til  =  4,  n2  =  4,  and  ns  =  11. 

From  Figure  6.5  and  Figure  6.6,  it  should  be  noted  that  the  significant  differ- 
ence in  the  utilization  of  the  ARP/Location  server  between  the  initial  partition  cind 
the  optimal  partition  comes  from  much  less  inter-cluster  data  traffic  in  the  optimal 
partition.  Thus,  the  sensitivity  of  data  messages  to  the  density,  velocity,  and  in-call 
probability  variations  dominates  the  utilization  of  a  network  component  in  the  virtual 
cell  system. 

For  data  messages,  the  increases  in  the  density  and  in-call  probability  of  mobile 
hosts  directly  affect  the  arrival  rate  of  data  messages.  However,  the  increase  in  the 
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Figure  6.6.  The  Utilization  of  the  ARP/Location  Server  for  Each  Type  of  Message 
When  an  Optimal  Partition  Is  Used 

velocity  of  mobile  hosts  does  not  raise  the  arrival  rate  but  the  forwarding  rate  of  data 
messages.  Thus,  as  the  velocity  of  mobile  hosts  increases,  the  utilization  of  data  mes- 
sages at  the  ARP/Location  server  is  slightly  increased  by  the  inter-cluster  forwarding 
data  messages,  but  it  could  be  neghgible  compared  to  that  of  data  transfer  class  of 
messages.  The  utilization  of  the  ARP/Location  by  handoff  messages  is  more  sensi- 
tive to  the  velocity  variation  than  the  density  variation.  The  in-call  probability  does 
not  affect  handoff  messages  because  it  only  raises  the  arrival  rate  of  data  messages. 
Compared  to  data  and  handoff  messages,  the  utilization  of  the  ARP/Location  server 
by  address  resolution  messages  could  be  negligible. 

The  system  response  time  Qr  for  each  type  of  messages  are  shown  in  Figure  6.7 
and  Figure  6.8.  In  the  optimal  partition,  the  system  response  time  for  data  and 
address  resolution  messages  is  reduced  at  the  cost  of  slightly  higher  response  time  for 
handoff  messages.  This  is  due  to  the  impact  of  handoff  multicast  operations  on  the 
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Figure  6.7.  The  System  Response  Time  for  Eacli  Type  of  Message  When  an  Initial 
Partition  Is  Used 
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Figure  6.8.  The  System  Response  Time  for  Each  Type  of  Message  When  an  Optimal 
Partition  Is  Used 
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relatively  larger  cluster  size  with  respect  to  the  cluster  size  of  the  initial  partition,  as 
described  previously.  When  considering  high  mobility  and  data  parameters,  i.e.,  2000 
mh/km^,  30  km/hr,  and  0.8  in-call  probability  of  mobile  hosts,  the  mean  handofF 
response  time  80  msec  is  quite  acceptable.  Given  some  mobility  and  data  traffic 
parameters,  the  analysis  shows  a  trade-off  between  the  data  and  handofF  response 
times.  The  much  larger  volume  of  data  traffic  over  handofF  traffic  in  mobile  data 
communications  reveals  that  even  the  small  difference  in  the  data  response  time 
could  improve  the  overall  system  performance  significantly.  To  reduce  the  handofF 
response  time,  more  strict  constraint  on  the  cluster  size  is  needed  in  the  optimization 
process. 

6.4    Chapter  Summary 

To  analyze  the  performance  of  the  virtual  cell  system  we  adopt  a  BCMP  open 
multiple  class  queueing  network.  Since  the  same  type  of  message  may  have  a  different 
routing  behavior  depending  on  which  base  station  belongs  to  which  cluster,  a  routing 
chain  is  defined  for  each  type  of  message  generated  by  each  base  station.  Both 
mobility  and  data  traffic  patterns  among  base  stations  and  the  topology  of  the  virtual 
cell  system  are  used  to  determine  service  transition  probabilities  in  the  queueing 
network  model.  With  various  performance  measures  such  as  the  utilization  of  network 
components  in  the  virtual  cell  system  and  the  system  response  time  for  various  types 
of  messages,  we  have  conducted  sensitivity  analyses  of  those  performance  measures 
as  mobility  and  data  traffic  parameters  vary.  We  also  compared  the  performance 
measures  of  two  different  virtual  cell  systems  which  are  one  deployed  according  to  an 
initial  partition  used  for  the  optimization  process  and  the  other  deployed  an  optimal 
partition  with  respect  to  the  initial  partition,  respectively. 


CHAPTER  7 
CONCLUSIONS 


We  have  presented  the  network  infrastructure  of  the  virtual  cell  system  for  the 
transmission  of  IP  datagrams  in  mobile  computer  communications.  With  the  virtual 
cell  concept,  physical  cells  are  grouped  into  larger  logical  cells  where  host  mobility 
is  shielded  from  the  IP  layer.  In  order  to  achieve  this  concept,  we  have  designed  the 
virtual  cell  protocol  which  consists  of  handoff,  address  resolution,  and  data  transfer 
modules,  based  on  the  distributed  hierarchical  location  information  of  mobile  hosts. 
It  eliminates  the  need  of  IP-level  protocols  for  host  mobility  and  provides  the  flexible 
coverage  area  of  a  virtual  cell  that  can  be  properly  engineered  according  to  mobility 
and  communication  patterns  among  physical  cells. 

Mobility  and  communication  patterns  among  physical  cells  can  be  represented 
by  the  move  and  find  frequencies  among  base  stations  because  base  stations  are  re- 
garded as  traffic  sources  and  destinations  from  the  prospective  of  the  virtual  cell 
system.  Given  the  move  and  find  frequencies  among  base  stations,  the  optimal  de- 
ployment of  the  virtual  cell  system  is  equivalent  to  the  problem  of  finding  an  optimal 
partition  of  disjoint  clusters  of  base  stations.  For  each  type  of  operation,  inter-cluster 
communication  is  more  expensive  than  intra-cluster  communication.  The  objective  is 
to  minimize  the  total  communication  cost  for  a  sequence  of  move  and  find  operations 
in  the  entire  system.  The  optimization  problem  differs  from  general  graph  partition- 
ing problems  in  that  it  additionally  considers  the  underlying  topology  of  base  stations 
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such  as  the  linear  arrangement  of  base  stations  in  highway  cellular  systems  and  the 
hexagonal  mesh  arrangement  of  base  stations  in  cellular  systems. 

For  highway  cellular  systems,  we  have  showed  that  the  optimization  problem  is 
solvable  in  polynomial  time  of  (9(mn^)  by  dynamic  programming  for  an  arbitrary 
number  of  clusters  and  size  of  a  cluster,  where  n  is  the  number  of  base  stations  and 
m  is  the  number  of  clusters  in  a  partition.  The  algorithm  also  finds  all  valid  partitions 
in  the  same  polynomial  time  if  given  a  constraint  on  the  size  of  a  cluster  and  the 
total  allowable  communication  cost  for  the  entire  system. 

For  hexagonal  cellular  systems,  several  heuristics  are  considered  using  the  tech- 
niques of  interchanging  or  moving  boundary  nodes  between  adjacent  clusters.  Exper- 
imental testing  and  analysis  show  that  unlike  general  graph  partitioning  problems, 
some  heuristics  relying  on  only  one  of  the  techniques  are  no  longer  useful  for  the  op- 
timization problem.  We  have  developed  several  heuristics  ba^ed  on  the  combinations 
of  the  techniques  of  moving  or  interchanging  the  boundary  nodes  between  adjacent 
clusters.  These  heuristics  produce  optimal  partitions  with  respect  to  the  initial  par- 
tition obtained  randomly  or  by  centering.  The  heuristics  are  compared  and  shown  to 
behave  quite  well  through  experimental  testing  and  analysis. 

Finally  we  have  analyzed  the  performance  of  the  virtual  cell  system  under  the 
assumption  that  it  is  deployed  according  to  the  topology  of  the  optimal  partition  such 
that  each  cluster  corresponds  to  a  virtual  cell.  The  virtual  cell  system  is  modeled 
as  a  BCMP  open  multiple  class  queueing  network.  In  addition  to  the  move  and 
find  frequencies  among  base  stations,  the  topology  of  the  optimal  partition  is  used 
to  determine  service  transition  probabilities  and  the  arrival  rates  of  each  type  of 
message.  By  solving  traffic  equations  of  the  network  model,  we  obtain  the  interesting 
performance  measures  such  as  the  network  response  time  for  each  type  of  message 
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and  the  utilization  of  the  base  station  networks  and  the  backbone  network  of  the 
virtual  cell  system. 
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